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On the irrational root sieve sequence 
Xiaowei Pan? and Baoli Liu ¢ + 


{+ Department of Mathematics, Northwest University, Xi’an, Shaanxi, P.R.China 
t Xi’an Aeronautical Polytechnic Institution, Xi’an, Shaanxi, P.R.China 
Abstract Let a be any fixed square free number, the irrational root sieve sequence is defined 
as taking off all k-powers (k > 2), of all square free numbers a from the natural numbers 
except 0 and 1. The main purpose of this paper is using the elementary method to study the 
properties of the irrational root sieve sequence, and give an interesting asymptotic formula 
for it. 


Keywords Irrational root sieve, square free numbers, asymptotic formula. 


§1. Introduction and Results 


For any positive integer n > 1, we call n be a square free number if for any prime p with 
p|n, then p? + n. Now we define the irrational root sieve sequence as follows: 

Definition. From the set of all natural numbers (except 0 and 1), 

-take off all powers of 2*, k > 2 (i.e. 4, 8, 16, 32, 64, ---); 

-take off all powers of 3°, k > 2; 

-take off all powers of 5", k > 2; 

-take off all powers of 6", k > 2; 

-take off all powers of 7", k > 2; 

-take off all powers of 10", k > 2; 

- and so on (take off all k-powers (k > 2), of all square free numbers). 


Now, we can get all irrational root sieve sequence: 

2, 3, 5, 6, 7, 10, 11, 12, 138, 14, 15, 17, 18, 19, 20, 21, 22, 23, 24, 26, 28, 29, 30, 31, 33, 34, 
35, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 
63, 65, 66, 67, 68, 69, 71, 72, 73,---. 

In reference [1], Professor F.Smarandache asked us to study the properties of the irrational 
root sieve sequence. About this problem, I do not know whether someone had studied it. At 
least I have not seen any related papers before. In this paper, we use the elementary method 
to study the this problem, and give an interesting asymptotic formula for it. That is, we shall 
prove the following conclusion: 

Theorem. Let A denotes the set of all elements of the irrational root sieve sequence. 


Then for any real number x > 1, we have the asymptotic formula 


Sris2-Sve- Set +0 (et Ima). 
nee T T 
acA 


2 Xiaowei Pan and Baoli Liu No. 4 


§2. Proof of the theorem 


In this section, we shall use the elementary method to complete the proof of the theorem. 
First we need the following simple lemma. 
Lemma. For any real number t > 1, we have the asymptotic formula 
6 
Y la(n)| = t+ 0 (vi), 
nxt 
where ji(7) denotes the Mobius function. 
Proof. For any real number t > 1 and positive integer n, from the properties of the Mobius 
function ju(n) ( See reference [3] ): 
e(n)| = D5 u(a) 
d?|n 


and note that 


we have 


dole = SEM e@= YF ua) 


nt n<t d2|n md?<t 
= SS nla) 1 
d<vt my 
t 
= > ud) = +0(1)) 
d<vt 
fd L(d 
= (49 O) 401 Xia) 
d=1 d>vVt d<vt 


This proves Lemma. 
Now we use this Lemma to complete the proof of our theorem. From the definition of A 
and the above Lemma we have 


Yi se Visz- OY DY vol 


a<e@ ak<a 2<k< PS 2<a< Ye 
acA |(a)|40 

k>2 

a>2 


lI 
8 
| 
= 
SS 
| 
aM 
= 
= 
+ 
1) 
a 
8 
IK 
a 
8 
NY 


1 
2<as Ve 2<a<ax3 
6 6 4 1 
= @ Vx x3 4 O (« Inz) 
1 m2 


This completes the proof of Theorem. 
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§3. Some notes 


If we use the deeply result for S- |(n)|, then we can get a sharper asymptotic formula for 
nt 
the mean value » 1. For example, assume RH, then (See Theorem 3 of reference [4]) 


a<ax 
acA 


6 sate 
2 la) = —t+0 (a+), 


where e€ be any fixed positive number. 


Using this result we may immediately deduce that assume RH, then we have the asymptotic 


formula 6 6 6 6 6 
1 1 1 1 _ 
) 1=2% svi 5? ae au de +0 (x2), 
7 T T T T 
aca 
acA 


where e€ be any fixed positive number. 
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Abstract The definition of sub-self-conformal set is given, and we obtain the formula for 
their Hausdorff and box dimension. 


Keywords Sub-self-conformal set, symbolic space, Hausdorff dimension, box dimension. 


81. Introduction 


Self-similar sets presented by Hutchinson [1] have been extensively studied. See for example 
[2] [3] [5] [8]. Let’s recall that. Let X C R” be a nonempty compact set, and there exists 
0<C <1 such that 

|w(x) — w(y)| < Cla — yl, Va,y eX. 

Then we say that w: X — X is a contractive map. If each w,;(1 < 7 < m) is a contractive 
map from X to X , then we call (X, {w;}%,) the contractive iterated function systems(IFS). 
It is proved by Hutchinson that if (X, {w;}™,) is a contractive IF'S, then there exists a unique 
nonempty compact set E C R”, such that 


E=|)wi(£). (1) 
j=1 
Set F is called a invariant set of IFS {w;}, . If each w, is a contractive self-similar map, then 
we call (X, {w;}%%,) the contractive self-similar IFS. Set E in (1) is a invariant set of w; , and 
it is called a self-similar set. If each w; is a contractive self-conformal map, then (X, {w;}/2,) 
is called the contractive self-conformal IFS, and set E decided by (1) is called a self-conformal 
set..[4] [6] 
Recently, Falconer introduced sub-self-similar sets. Let w; is self-similar IFS, and let F' be 
a nonempty compact subset of R” such that 


This set F' is called sub-self-similar set for w;. Easy to see that, self-similar sets is a class of 
special sub-self-similar sets. At the same time, Falconer obtained the formula for the Hausdorff 
and box dimension of the sub-self-similar sets, when w; satisfies the open set condition (OSC). 
In this paper, we will give the definition of sub-self-conformal sets similarly, and obtain the 


formula for their Hausdorff and box dimension. 
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At first, we give some definitions and remarks. By a symbolic space we mean the in- 
finite product space © = {1,2,---,m}N, ie. © = {(i1,é2,---) 1 <i; < m Let Q = 
J 41,2,--- ,m}”. For I = (i1,%2,--+ ,im) € E and k € Nt, we write I, = (i1,%2,--- , ix), 
n>0 
TP = (ingastegas*+), Ue = {1|, : I © DS}. The shift transformation on © is defined by 
6(1) = 1|'. For I = (ii, i2,---), J = (j1,J2,-++) € U , define a metric d on © by 


d(I, J) =e*,k = max{l: I|; = Jj}. 


Obviously, if J = J, then d(I, J) = 0; and if 4 Ji, then d(J, J) = 1. 

Let {w;}%, define as above, and B be any closed ball in R” large enough to ensure that 
wi(B) C B,i = (1,2,--- ,m). Since each w; is contracting, it is easy to see that the sequence 
of balls wy), (B) := wi, ° Wiz 0+ ++ 0 W;,(B) decreases with k and has intersection a single point. 
Then we let 7: © — R” be defined by 


(1) = () Wi, 0 Win --- Wi, (B). 
k=1 


Or 


m(E) = lim wi, © Wig ---wa,(z), V2 ER”. 


Definition 2.1. Suppose (X, {w;}!,) is a contractive self-conformal IFS, and F' be a 
nonempty compact subset of X satisfies the condition 


Then we call F' a sub-self-conformal set for {w;}?7). 


§2. Symbolic space 


For studying the Hausdorff and box dimension of the sub-self-conformal set, we set up 
a Closed relationship between a sub-self-conformal set and a compact subset of the symbolic 
space. 

Proposition 3.1.[4] Suppose 7 is defined as above, then the map 7 : & — R” is continuous. 
Similarly as the case of sub-self-similar set, we have this conclusion. 

Proposition 3.2. F is a sub-self-conformal set for the IFS {w;}1”, if and only if there is 
a set A, which is a shift invariant closed subset of 4, such that F = 7(K). 

m 


Proof. Let F' be a sub-self-conformal set for {w;}”™,, then F’ C U w;(F). We define Kk 
i=1 
as follows 


K = { (in, %2,---) : (ik, tkga, ++) € F,WK € Zt}. 
Obviously, for any I € K, we have 6(I) € K. Then we will prove that F = 7(K). By the 


definition of K, we have m(K) C F. Since F C U w;(F), then for any vo € F, we can find 
i=1 
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some 21 € F' such that xp = w;, (21), for some i; : 1 < i; < m. Similarly we have x1 = w;, (x2) 
for some xg € F and 1 < ig <‘m, and so on. We get a sequence {x,} C F andi,:1 <i, <m 
satisfy the condition 


co 
rp-1 = lim wy, 0 Wi,,, 077° Wi, (En) = ih) Wi, o++-0 w;,(B) Er Vke TD. 
I= 


So, to € F. 
Conversely, let K be a compact subset of ©, and @([) € K for any I € K. Suppose 
F =7(4K), then we will prove that FC U w;(F). In fact, for any « € 7(K), there exist some 


i=1 
I = (t1,%2,--+) € K& such that 


so m(EX) C U wi(m(K)), ie. F is a sub-self-conformal set for {w;}}"",. 
i=1 


From {8][9], we can find a similar proof. 
Next, we will define a measure M*(A) on the symbolic space. By using this measure we 
can prove that: 7(s) = jim ( S- jw (x)|*)* =1, then 
ee IEK, 


dimy F < dim, F < dimgF <s. 


Moreover, if M*(K) < co, then H*(F’) < oo, where 1(K) = F. Finally we apply this theorem 
to attain the main conclusion of this paper: Let F be a sub-self-conformal set for {w;}), which 
satisfy the open set condition, then s = dimy F = dimg F, here s is the number satisfy 


§3. Hausdorff and box dimension of sub-self-conformal sets 


To make the definition of sub-self-conformal sets plain or comprehensible, we give some 
examples. 

e.g.1. Suppose {w;}%, is a self-conformal IFS and a nonempty compact set F' satisfies 

m 


P= U w;(F), here m > 2, i.e. F is an invariable set for {w;}%,, then F is a sub-self-conformal 


i=1 
set for {wi }?4. 
e.g.2. Let F' be defined as e.g.1 and let OF be the boundary of F’, then OF is a sub-self- 


conformal set for {w;}/”). 
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Since for each x € OF, we have x € F. So there exist some i: 1 < i < m, such that 
x €w;(F). At the same time, every neighborhood 6(2) of « contains points outside w;(F’). ie. 


x € Ow; (F) = w;,(O(F), 


so a 
OF Cwi(OF) C (J wi(OF). 
w=1 
e.g.3. Let F be a sub-self-conformal set for {wi,we,---,Wm}, and E be a sub-self - 
conformal set for {51, 82,--- , Sn}, then FU E is a sub-self-conformal set for 
{w1, W2,°+* Wm, $1, $2,°°* Sn}. 
e.g.4. Let F be a sub-self-conformal set for {w1,w2,--: ,Wm}, and E is a compact subset 


of F, and w is a contract self-conformal map, then F' U w(£) is a self-self-conformal set for 
{w1, W2,°°+ ,Wm, wh. 
Lemma 1./6] Suppose X and {w;}%%, be defined as above, and log |w/(x)| satisfy Dini 
condition, let ry = inf |w?(x)|, Rr = sup |w{(x)|. 
TEX eX 


(i) There exist C, such that 
Rr<Cirr, forany le. 
(ii) There existC2 > 0 such that for any x,y,z € X,|a — y| < 0, 


w7(z) — wry) 
oy 


Cz *|w7(z)| < <Cylw}(z), for any 1 EQ. 


(iii) There exist C3 > 0,such that for any x,y € X, 
|wrs(2)| < C3|w7(x)||w5(x)| for any I, JE Q. 
Next, we consider the dimension of the sub-self-conformal set. Let 
K = {(i1,i2,+++) : wigs iegi,:++) € Fy Wk E ZT} CY, 


and Ay, as above. For some x € F,, we write ry = |w{(x)|, here I = (%1,%2,:-- 2%). By lemma 1 
(iii), we know that for any I, J € Q there exists C such that r7y < Cryryz , hence there exists 


C" such that 
Mews DD eC OL: 


LEK R41 TEK, &JEK, TeK, IEk;, 
Lemma 2.[3] Let {b,}22, be a positive real number sequence. Suppose for any positive 
br 
integer k, m we have beim < be + bm, then jim Ss exist. 
00 


Lemma 3. The limit 7(s) := lim ( S- r3)'/* is exist. 
k— oo IK, 
Proof. By lemma 1(iii), we have 


Carrz < (C3rr)(Csrz). (2) 
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We choose logarithm to (2), and let b, = log S- r7 + slogC3 . Through lemma 2, we can 


IeK, 
b 
prove that log( y reyk = — is convergent. Hence we know that lim ( y rz)+/* exist. 
k k—-00 
IeKy IEK 
Then we define a measure on the symbolic system. By using this measure we obtain the 


dimension of the sub-self-conformal set. We define 


M§ (A) = inf ie Ir? : ACen, UI = 7 ' 
I I 
here ACY, s >0,k =1,2,--- , ando; = {1J: J © S}, for I € Ky. Simility to the hausdorff 
measure we define 
M*(A) = jim M; (A). 


Propose 4.[10][11][12] Let A be a Borel set of 5. If 0 < M*(A) < oo, then there exists a 
compact subset Ap C A and a b > 0, such that 0 < M*(Ao) < oo, it follows that 


M*(AgNor) <br7, VI EQ. 


Propose 5. Let K,r; be defined as above, then we have some conclusion 
(i) There exists s; which satisfies ($1) = 1 ; 
(ii) so := inf {s > 0: M*(K) =0} =sup{s >0: M*(K) =o}; 


CO 


(iii) mint {s20. es ri <ce} = wp 20: s > i=} 


k=1IEK k=11€K 
At the same time we have s; = s2 = $3, and 


M*(K) > 1. 


Proof. At first we will prove that s; is existent. Note thatr = min {r;}, R= max {rj}. 
1l<i<m 1<i<m 


Then for h > 0, we have 


I 
kh IeKr kh 
aS 7 SR 
ry 
IEK,, 
Let k — oo , it follows that 

h TLE) phy, 

T(s1) 


5, >0, h>0. 
Hence 7 is continuous and strictly decreasing. Moreover, 7(0) > 1 and 7(s1) — 0 if 


8; — oo. So there is a unique s; > 0 satisfies r(s;) = 1. By the standard feature of measure 
and progression [13][14], we know that 


inf {s >0: M°(K) =0} =sup{s >0: M*(K) =o} 


and 


it {1> 0:5. Yo vi <aoh amp{a2 0:30 Siac} 


k=1IcCK, 
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We write them s2 and s3 apartly. Next, we will prove that s; < sg. Suppose M*(K) <1, 

for some s > 0. Here we have s > sg. Then there is a covering U o; of K such that 5° re<il. 
IeEQ IeQ 

Hence there exists t: 0 < t < s, such that S- ry <1. Since K is compact, we can find a finite 


TEQ 
Q. Write g = max {|J| : J € Q}. Let 


Qn ={hlh---In: 1; EQ, |IyI2+-+Ip-1| <k< |I,I2---Ip 


}, 


here k > q, and |I| denote the length of I. Let K be a positive integer, we have K No; C 
U (K Norz,7) , for I € Ky. Thus if I € kK, , then IJ € Q;, for some J € Q. Hence, we have 


JEQ 
that for each k 


) rt <r ot ) re <rt, 
TeK, V'EQ: 


Thus, if M*(ic) < 1, then for some 0 < t < s we have 


So t(s) < T(t) < 1, ie. s > 51. It follows that s1 < sg. sg < s3 is obviously. We cover Kk 
CoO Co 


by cylinders o7, here I € U Ky. Then if S- S- rt <0o, we have Mp(K) < o S- r+ for 
k=0 k=11€Kp j=1IEK; 


k — co. Thus, if s > 53, we have s > so. 
[oe) 


To see that s; = s3, we observe that S- S- r7° converges if T(s1) < 1 and diverges if 
k=11€K 
T(s1) <i. 


In conclusion, we have s1, 52,53 exist and s; = s9 = 83. 

On the side, we can know that if M*!(i) < 1, then r(s1) < 1. From the proof of s1 < se, 
we have that if 7(s1) =1 then M*(K) > 1. 

Propose 6. Let F' be a sub-self-conformal set, and 7(s) = 1, then 


dimy F < dim, F < dimgF < Ss. 


Hence, if M*(Ix) < oo, then H*(F’) < oo. 
Proof. Suppose B is a closed ball such that w;(B) C B, for any i(1 <i <m). Let 6 satisfy 
0<6< |B). For any I = (i1,%2,---) € K, we can find a corresponding k € Zt such that 


rd < |wy,(B)| = lw’, ©) |B] <6, VEE B. 
Write 
Q5 = {1 € Kx i rd < |w7(§)||Bl = |wr(B)| < 5}. 


We have that FC U w7(B) and U wy(B) is a cover of F’. If Ns is the number of sets in 


IEQs IEQs 
this cover, then for t > s, 


N50 = SO 5 <r * BY SE lw] < Cor * |B SS SO rr SM < co. 


IEQs IEQs k=1IE€K, 
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It follows that dimg < t for allt > s , and so dimgF < t. 
Since kK C U o7, then FC U wr(B). So 
IEQ IeEQ 


H3(F) < C3|Bl’ My(K). 


Here 6 > r®. Thus if M*(K) < oo, then H°(F) < oo. 

Lemma 7.[12] Let {V;} be a collection of disjoint open subsets of R” such that each V; 
contains a ball of radius a,r and is contained in a ball of radius agr. Then any set U of diameter 
at most r intersects at most by := (1 + 2a2)"a;,” of the closures {Vi}. 

Theorem 8. Let F be a sub-self-conformal set for {w;}/", which satisfy the open set 
condition, and let 7(s) = 1, then H*(F’) > 0 and 


s=dimy F = dim,F = dimgF. 


Proof. We can know that M*(K) > 0 from propose 5. Let Ao be a compact subset of K 
which confirm by propose 4. The Borel measure js which supported by K is defined by 


ww: w(P) = M*(Aon P), 

where P C &. By the propose4 we have 
Wor) <borz, VIE Ky, kEN. 

Let 

((U) = w{ J: (J) €U}. 
For U C R”. Such jf is a Borel measure supports by F’. Let V be an open set satisfying OSC, 
and let U C R” satisfy 0 < |U| < |V]. Let 

Q = { (it, a tk) TH Tig Tin |V | < |U| s THTig Tinea \V|} . 


For any fix « € V. Since rC4|U| < |wr(V)| < C4|U| and using Lemma 7, we have that there 


are at most b; element in the family 
Qo := {lEQ:UNwr(V) #¢}. 


Thus, if w(J) € U, then there exists k € N such that J|, € Q, and so J € a; for some I € Qo. 
Hence 
j(U) = So wf Teor} <b Yo v7 < db, |V|-* OI’. 
IEQo IEQo 

Since ji is supported by F’, the mass distribution principle implies that H*(F’) > 0. So s = 
dimy F = dim, F = dimpF, which s is the number satisfying T(s) = 1. 
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Abstract In this paper we shall study several kinds of mean values of the multiplicative 


arithmetic function D(n) by the convolution method. 


Keywords Arithmetic function, mean value, convolution method. 


§1. Introduction and Results 


We define the arithmetic function D(n) by 


ie ifn=1, 
D(n) = II ap*', ifn>1. 
pe|[n 


In her doctoral thesis, Wang Xiaoying [3] proved the asymptotic formula 


where 


_ 6 p 1 1 
om all Osea Ose) 


In this short paper we first show that the asymptotic formula (1) can be slightly improved 
by the following 
Theorem 1. There exists an absolute constant c > 0, such that the asymptotic formula 


YS pay ott Oweter a) (2) 


n<ux 


1 
~ 3 


is true, where 6(x) := log® x(log log x) 

Remark. It is very difficult to improve the exponent 1/2 in (1) and (2), unless we have 
substantial progress in the study of the zero region of ¢(s). Therefore it is reasonable to study 
the problem in short intervals. In this case we have the following 


, This work is supported by National Natural Science Foundation of China (Grant No. 10771127) and 
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Theorem 2. If 15 *8° < y < x, then 
1 =e 
a<n<at+y 
Finally we prove the following Theorem 3, which studies mean values of the function 
log D(n). 
Theorem 3. We have 


> log D(n) = cox + O(a? **), (4) 
n<u 
where 
1 co 
CQ i= 5 (: - ~) be mi log(np"~'). 
p n=2 


Notations. Throughout this paper, ¢ > 0 denotes a small positive constant, s(n) denotes 
the Mobius function, 6(2) := log? a(log log a)~ 5. 


§2. Proof of the theorems 


ee 
We first prove Theorem 1. Let f(s) = S- (Res > 1). By the Euler product we 


24 Denne 
get 
f(s) = I] (1 a < ii sie aT a ) 
=I) bP (6 (1 see) + gs) 
= Fay » 


1 1 
where G(s) = II c | 7am (108 (1 7 = 7 ‘=)}: 
Pp ps . . 


co 
Let G(s) = S- gr) It is easy to prove that the series is absolutely convergent for 


n=1 
Res >¢> 0. So we have 


Ye | 9m) |< a. (6) 
n<u 
From (5) we see that the relation 
1 
Dw = SoH Da0 (7) 
dl=n 
holds. 
We have the well-known asymptotic formula 
6 1 e96(a 
DL wd) = yt + O(e2 eH“), (8) 


n<u 
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where cp > 0 is an absolute constant. 
Now Theorem 1 follows from (6)-(8) and Theorem 14.2 of Ivié[2]. 
In order to prove Theorem 2, we need the following result(see [1]): If 3+? < y < «, then 


we have 6 
~~ HO) = jut Oya). (9) 
u<d<a+y 
Lite 
Suppose x5 t3* < y <a and let M = x”. For 1 < M, we have (=) an ; < 2 By (7), 
we get 
Lay = LY POs 
D(n) 
a<n<at+y a<di<at+y 
= So a) y?(d) 
l<a+y Bed< ate 
= So gl) wdyt S> gl) u(d) 
1<M ® cdc B4¥ M<i<a+y # cdc +4 


=yo+>5 (10) 


say. 
Using (9), (10) and partial summation, we get 


Di 


S = dao(arre(r(7) )) 


= © Gy +0 (oy: atl!) +0 | ya-€ S- oO 
<M 


I 
Q 
=< 
S 


l>M 
=cyt O(yx“), (11) 
SS «K YS [gy +1) « yo. (12) 
2 M<iI<a+y 


Theorem 2 follows from (10)-(12). 
Finally we prove Theorem 3. Suppose 4 is a complex number with —e < RA < e. Consider 


1 
F(s,A) = S- Din) *n* (Rs > 1). By Euler product, we have 
n=1 
1 1 1 
F(s,X = 1 } pase sae 
a II ( D*(p)ps_— D(p?)p?s Dippy 
7 i ee ee 1 i 
_ Il WF ps : QD p2str ! 3rp3st2r ia nApns+(n—1A mys 
Pp 
a 1 1 1 1 
“U0-m) [0-p) +p tape ape) 
Pp Pp 


Pp 
= ¢(s)K(s, A), (13) 
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where 
1 1 1 1 
K(s,A) =|] (1- =) (1+- pun Pree pret ain -) (14) 
Pp 
k(n, A) es a 
Let K(s,A) = S- . It is easily to see that the series is absolutely convergent for 
n> 
n=1 
1 RA 
Rs > ee So we have 
. k(n, d) < eV 2+RA/2+., (15) 
n<x 


From (13)-(14) and Theorem 14.1 of Ivié[2], we get 
YD) = h(A)e + O(e/74RN748), (16) 
n<u 
where h(\) = K(1, A). 
Taking the derivative for \ from both sides of (16), and then letting 4 — 0, we get 
S © log D(n) = —h!(0)a"/? + O(a'/?**). (17) 
n<x 
It is easy to check that 
1 = —n n— 
h'(0) = -»— (1 - .) a: log (np"~*) . 
p n=2 


This completes the proof of Theorem 3. 


References 


[1] M. Filaseta and O. Trifonov, The distribution of squarefull numbers in short intervals, 
Acta Arith., 67(1994), 323-333. 

(2) A.Ivié, The Riemann zeta-function John Wiley & Sons, New York, 1985. 

[3] Wang Xiaoying, Doctoral thesis, Xi’an Jiaotong University, 2006. 


Scientia Magna 
Vol. 3 (2007), No. 4, 16-20 


On an infinite series 
related to Hexagon-numbers! 


Lingling Wang? and Yanni Liu? 
{ School of Mathematical Sciences, Shandong Normal University 
Jinan, 250014, P.R.China 


t Department of Mathematics, Northwest University 
Xi’an, Shaanxi,P.R.China 
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§1. Introduction 


For any positive integer m, if n = m(2m — 1), then we say n is a Hexagon-number(see 
reference [1]). For any positive integer n, let m be the largest positive integer which satisfies 


the inequality 


m(2m —1) <n < (m+1)(Qm+1). 


Define a(n) := m(2m — 1), and call a(n) the Hexagon-number part of n. The second-named 


1 
author[1] proved that for any real s > 1 the infinite series f(s) = S- OC) is convergent and 
a(n 
n=1 
5 
that f(2) = ral —4ln2. 
In this short paper, we shall prove the following 
ae 
Theorem. The Dirichlet series f(s) = S- Poe can be continued to the whole complex 
a(n 
n=1 


plane as a meromorphic function with simple poles 1, 1/2, —1/2, —3/2, —5/2, --- 
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§2. Proof of the theorem 


In this section, we shall complete the proof of Theorem. First let 


4m+1, ifn=m(2m-—1), 
b(n) = 
0, otherwise. 


Then, if s = 0 > 1, we have(see reference[1]) 


=< —  4m+1 b(n) 
P(s) = aS(n) = oa @m =p aS ns 


1 
For any real u, let [u] denote the greatest integer not exceeding u, bi(u) = u— [u] — 5 We 
define a series of functions by,(u) (k > 2) by the following relation: 
U 1 
bi +1(u) — be+1(0) = by (t)dt, } beti(t)dt =0,k =1,2,---. 
0 0 


Then for each k > 1 the function the b,(u) is a periodic function with period 1. Especially we 
have (see [2]) 


1 2 LI vol 
For any u > 1, define B(u) := s b(n), Uo = : jeg a Then 
— —. 4 2 16 
Bu) = S° (4m+1) 
m(2m—1)< 
= > (4m+1) 
m<uo 


= 2[uo] ({uo] + 1) + [uo] 


=2 | u = hi (ag) ‘| + ug — bi (uo) 


5 
= Qua + Ug — 6 + 4b2(uo) = 4ugb1 (uo) = by (ug) 
= By(u) + Bo(u), 


where 


B,(u) := 2u3 + up — a 
Bo(u) = 4b (uo) = 4ugb; (uo) = by (uo). 


Suppose x > 1 is a fixed real number, we write 
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Using Abel’s summation formula, we have 


ae er 


n>x 


= - = dBi (u) +f = dB,(u). (2) 


Now we evaluate the integrals | = dBy(u) ) and [ “dBa(u ), respectively. Our main 


tool is partial summation and the power series. 


ie dBiu) +4 TOF a, 


us 


are 


qi-s du 
= + lim 2 
s—l M-—oo _ us/8u +1 
gi-s [ du 
= lim 2 ai 
s—1 Moo J, u8(8u)2(1+ Ly)? 
qgi-s 7 1 M du 
— im 
s—l Moo 4/2 x wta(1+ 1)2 
a | a ee (Q2n—1)! 1.) 
> s—l Moo Dd x usta Dd ) (2n)! Su) a 
gis 1a, 1, @e—DW 7? 4 
= li n d 
s—1l hanes 5) d 3) (2n)!! : yntsts a 
y= tieay LenS 1 1 
= li n M27-s-" 5—s—n 
s—l1 viene 5 | ry (2n)!! aa —& ) 
oe 2. I, w= Ie 
— Le eee dC ) 


(Qn)! $—s—n 


1, 1,,(2n-1)!! 2278" 
aes 3)" (Qn)! $—s—n 


n—1)! 1 i 


1 = (2 

iin 1 Mi 1)” n 

ae eet geen 

= 0 5—s—n 

Recall Rs = a > 1, so 
— (Qn-1)!, 1 1 oY 

1)” n tin 1 
DI (2n)!! (sa7? To san * Xulgy) <1; 


and correpondingly 


ic (Qn—1)!! 1 1 1 
M?-* 1)" nm M?-°, 
: 2 ony ‘Baa? eal 
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Thus we get(o > 1), 


say. It is easy to see that the power series in (3) is absolutely and uniformly convergent in any 
compact domain in C \ {1/2, —1/2, —3/2, —5/2, ---}. So f,(s) represents a meromorphic 
function on C with simple poles 1, 1/2, —1/2, —3/2, —5/2, 


~ dB 
Now consider i. dBa(u) By partial integral, we have 
zx us 
~ dB B ~ &B 
i a(u) Bau) I | 2(u) 4, 
‘ us us * ystt 
Bo(x) °° b2(uo) °° ugbi (uo) °° by (uo) 
= a +4s f ett du —4s [ ueet du >| ees) du 
B 
= (x) as [ as [ sf (4) 
xs 1 2 3 
4u—-1 1 1 
We study only |. Let gi(v) = Ge yer z= 7 ; 16? and k > 1 be an integer. 
By the change of variable v = up and repeated partial integration we get 
be(u 4v — 1) 
| =[ 3 oS —v)stt Oa oeen 
=f bo(o)ai(o)ae 
ro 
= pi(opba(e)|age — f gh(vba(ude 
zo 
= =pu(croba(v0) — f of(w)aba(v 
xo 
= —g1(2o)b3(ao) + g(a) ba(to) — gi (#0)bs (x0) +--+ + (—1)¥g{*~? (wo) be-+2(0) 
+(-1)* f dnsa(v)al (wae. (5) 
nae) 


For v > 1, we have the power series expansion 


gi(v) = (4v — 1)(202)-*-1 (1 7 i) —s-1 
= (4p = 1)(202)- 17 (s+1)(s+2)-»-(s +n) F ) 


n! 2 
n=0 @ 


fore) . tants ra 2s+n+2 
= (io SoG tMED- etm (1 


U 
n=0 


which implies that g\") (v) < vu 2e-k-1_ Thus the integral / is absolulely convergent for 
1 


k 
a> <5 and represents an analytic function in this half plane. Similarly, the i: and i are 
2 3 
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analytic in the range o > —k/2, too. It is easy to check that the other terms in (5) form an 


°° dBo(u 
entire function in the whole plane. From the above we see that / ebay) has an analytic 
u 


s 
x 
dot ka, er ; * dBo(u) : a 
contiuation to 0 > = Since k is arbitrary, the integral ——— has an analytic contiuation 
to C. 


Combing(1), (2), (3), (4) and (5) completes the proof of Theorem. 


s 
x U 
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81. Introduction 


Recall that a semigroup S is called an rpp semigroup if all its principal right ideals aS, 
regarded as right S'-systems, are projective. 
According to J.B. Fountain[5], a semigroup S' is rpp if and only if, for any a € S, the set 


Ma={e € E| Sta C Se and for all x,y € S!, ax = ay > ex = ey} 


is a non-empty set, where £ is the set of all idempotents of S. An rpp semigroup S is called 
strongly rpp if for every a € S, there exists a unique idempotent e in M, such that ea=a. It is 
easy to see that regular semigroups are rpp semigroups and completely regular semigroups are 
strongly rpp semigroups. Thus, rpp semigroups are generalizations of regular semigroups. A 
strongly rpp semigroups S is said to be a right C-rpp semigroup if £* V R is a congruence on 
Sand Se C eS for all e € E(S). 

It is clear that a right C-rpp semigroup is a generalization of a right inverse semigroup (see 
[8]). Right C-rpp semigroups have been investigated by Guo and Shum-Ren in [3] and [2]. 

In this paper, we will give another construction of such semigroups by using right cross 


product of semigroups. 


§2. Right cross product structure 


We first introduce the concept of right cross products of semigroups. 

Let Y be a semilattice. Let M7 = [Y; Ma, Ya,g] be a strong semilattice of left cancellative 
monoids M, with structure homomorphism yg,g and A = Ugey Aq be a semilattice decompo- 
sition of a right regular band A into right zero bands Ag. 


'This research is supported by National Natural Science Foundation of China (Grant No:10671151) 
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Denote the direct product of a cancellative monoids Mg, and a right zero band Ag by 
So = Ma x Na- 

For any a, 8 € Y with a > (, we define a mapping 09,3: Ma x Ag — Ag by (a,j) # [a,j] 
satisfying the following conditions: 

(i) If a = @, then for any a € Ma, j € Ag, [a,j] is constant; 

(ii) If (a,2) € Sa, (b,7) € Se, k € Ay, then [ad, [a, iz]¢k][a, i7]i = [a, o[b, pk] J]: 

(iii) If (a,7) € Sa, 7 € Ag, k € Ay and [a,ij]¢ = [a, ik]2, then [la, 27] = [la, ike. 


Then, consider the set S = Usey Sq with a binary operation “o”: for any (a,i) € Sa, 
(0,3) € Sp, 


(a, i) : (b, j) _ (ab, [a, ijl), 


where ab and 77 are the usual semigroup products of a,b and 7,7 respectively in the semigroups 
M and A. Using conditions (ii) and (iii), we can verify that (S,0) is a semigroup. Denote this 
semigroup by S = M ™¢ A and call it a right cross product of M and A. 

To give a new construction of a right C-rpp semigroup, we list a description for such 
semigroups. 

Lemma 2.1.[2] Let S be a rpp semigroup. Then the following statements are equivalent: 

(i) S is a right C-rpp semigroup; 

(ii) S is astrongly rpp-semigroup, £* VR is a congruence on § and Se C eS for all e € E(S); 

(iii) S is a semilattice of direct products My x Aq, where M, is a left cancellative monoid 
and A, is a right zero band for every a € Y; 

(iv) S is a strongly rpp semigroup such that D™ is a semilattice congruence on S, and 
D |regs= R |Regs, where Reg S is the set of all regular elements of S; 

(v) S is semilattice of D-simple strongly rpp semigroups, and D | RegS= FR |Regs- 

Now we are ready to establish a construction theorem for right C-rpp semigroups. 

Theorem 2.2. Let Y be a semilattice. Let M = [Y; Mg, Ya,,] is a strong semilattice of left 
cancellative monoids M, with structure homomorphism yg and A=Uacy Aa be a semilattice 
decomposition of right regular band A into right zero bands Ay. Then, a right cross product 
M xo A of M and A is a right C-rpp semigroup. 

Conversely, any right C-rpp semigroups can be constructed in this way. 

Proof. To prove the direct part of Theorem 2.2, we first consider the set of idempotents 
of the semigroup M x9 A. Let EF = Uaey{(la,t) |i € Aa, la is the identity of M, }. Clearly, 
for any (la,?), (la,#) ° (la, 2) = (La, [la, tJi) = (a, 2). 

Conversely, if (a, i)? = (a,i) € Ma x Aq for a € Y, then it is trivial to check that a = 1a, 
where 1g is the identity of M,. This shows that F is the set of all idempotents of M x9 A. 

It is easy to see that M x9 A is a semilattice of direct products My x Ag. 

To show that M x, A is a right C-rpp semigroup, by Lemma 2.1, we only need to prove 
that the semigroup M x9 A is an rpp semigroup. 

Suppose that « € M x, A. Then there exists a € Y such that x = (a,i) € Sg = My x Ag. 
Hence, for any u = (b,j) € Sg’, v = (c,k) € Sy’, if zou = roy, that is, (a,7) 0 (b,j) = 
(a,i) o(c,k), then ab = ac clearly, and [a,ij]i = [a,ik]i. Since M is a strongly semilattice of 
left cancellative monoids M,, we have that b=. 
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Let e = (lq,7) € E. By the condition (iii), we have [ly,ij]i = [lo,ik]i and eou=eov. 
Again since (a, i) © (1q,7)=(a, [a, ii) = (a,), that is, r = xo e, then we have S'x C Ste. Thus 
M, 4 2. By the definition of an rpp semigroup, « € M 9 A is an rpp semigroup and so it is 
a right C-rpp semigroup. 

We now begin to show how a right C-rpp semigroup S$ become a right cross product of 
and A, where / is a strongly semilattice of left cancellative monoids and A is a right regular 
band. Our proof is divided the following steps: 

(I) Assume that S is a right C-rpp semigroup. By Lemma 2.1, S can be written as a 
semilattice Ugey(Ma x Aa) of Ma x Aa, where M,, is a left cancellative monoid and Aq is a 
right zero band. Denote Sg = Ma x Aa. 

Now we let M = Ugey M, and A = Ugey Ag. 

Firstly, we will show that M is a strongly semilattice of M,. For this purpose, let E = 
Uaey{(la,t)|7€ Aa, la is the identity of 1}. Then it can be easily checked that F is the 
set of idempotents of S. 

Now, for any a, € Y with a > £, we let(a,i) € Sa, (1s,j) € Sg NE. Clearly, 


(a,7)(1,, 9) = (b, k) € Sa, 


for some b € Mg and k € Ag. But we have that (a,i)(1g,j)(18,3) = (0,k)(1¢,7) € Sg and 
thereby 


(a,7)(1,, 7) = (0, j) € Sa. 


This implies that the choice of b is independent on j. Moreover, since F is a right regular 
band, we see that for any 1 € Ag, (4,!)(1¢,7) = (4,2)(1a, (16,9) = (a, 7) [a Ng, 7)]As,7) = 
(a,t)(1g,7). This also shows that the choice of 6 is independent on 7. 

Consequently, if we define a mapping Yo,g : Ma — Mg by a+ ava, if and only if 


(a,4)(1a,3) = (b,) = (aPa,6,5): 


Then the mapping Yq,g is well define. By routine checking, all the mappings ¢..,g are indeed 
the structure homomorphisms of a strong semilattice Y of monoids. Thus, M = [Y; Ma, ¥a,g] 
is a strong semilattice of My. 

Next, we claim that A is a right regular band. For this purpose, it suffice to show that A is 
isomorphic to E(S) which is the set of all idempotents of S. Since S is a right C-rpp semigroup, 
the mapping w: E(S') — A = UaeyAa given by (1a,7) + 7 for i € Aq is clearly bijective. We 
now define ij = k for i € Ag,j € Ag if and only if 


(la, 4)(1, 9) = (lag, k). 


Then, we can easily see that the set A = Uaey Aa, under the above multiplication, is isomorphic 
to E(S) which is a right regular band as mentioned above. 

(II) Now we consider a construction for structure mapping 64,3 of a right cross product 
M xo A. 

For any a, 8 € Y, with a > #, let (1g,7) € Sg ML, and (a,i°) € Sq, where 2° is a fixed 
element. Then by Lemma 2.1, we have 
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(1e,9)(a,0°) = (a’,k) € Sp 


and 
(12, 7)(a, )(1,, k) = (12, 7)(a, )(1g, k°)(1,, k) 
>= (1e,3)(@¥a,8,k°) (1g, k) 
= (ava,e,k). 
Consequently, 


(12, 3)(a,%°) = (ayva,g, k). 
Thus, we may define a mapping 64,3: Ma x Ag — Ag by (a,j) + [a,j] € Ag if and only if 
(1g, 7) (a, 7°) = (ava,s, [a,7]). By using the above result, we have 


(1e,9)(a,2) (1a, 7) (a, 0°) (a, 4) 
= (a%o,8; [4 J])(1g, [a,5]) (as *) 
= (4 Pa,¢;[4, 5]) (Ae; la, dé) 
= (4Pa,8; [4, J]*). 


It is easy to see that, if a = G for any a € Ma,7 € la, [a,j] = i°. This shows that the 
condition (i) is satisfied by mapping 04,2. 

(III) We will verify that the conditions (ii) and (iii) in the right cross product of M and A 
are satisfied by the mapping 0,,,. 

It follows by associativity of the semigroup product that the mapping 09, satisfies the 
condition (ii). 

Next, we will show that the mapping 04,g satisfies the condition (iii). 

Since a right C-rpp semigroup is an rpp semigroup, we know that for any x = (a,i) € 
S, M, ={e€ E | S'a C S'e and for all u,v € S1,2u = zv > eu = ev} is a non-empty set. 
There exists an idempotent e = (1q,i) € M, such that for any u = (b,7),v = (c,k) € S', if 
xu = xv, then eu = ev. From this, we can easily verify that the condition (iii) is satisfied by 
the mapping 04,2. 

(IV) It remains to show that the multiplication on the right C-rpp semigroup S coincides 
with the multiplication on the right cross product M x9 A and so S = M x9 A. 

For any (a,i) € Sa, (b,7) € Sg,a,28 € Y, by using the above result, we have 


(6,3)(a,4) = (6,9)(1, 7) (a, #)(@, 4) 
= (6,5)(6a,5t)(a, ¢) 
= (6,5) (1a, [4, Ji]i)(@Pa,8a, [@, Jé]4) 
= (696,80; (4, Ji]t)(aPa,ga, [@, Jé4) 
=  (b96,6a%Pa,pa; (4, Jt)t) 
= (ba, (a, jit) 
= (b,j) 0 (4,4). 


This completed our proof. 
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Abstract In this paper we recal the definition of non closed fuzzy algebraic structures namely 
fuzzy sub-halfalgebra of a BCl-algebra. Here we give the necessary and sufficient condition 
for the two fuzzy sub-halfalgebra to have open fuzzy extension.Also we recall yet another 
new notion of the concept called antifuzzy sub-halfalgebra and antifuzzy extension and obtain 
some characterization theorems about them. 
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§1. Introduction 


The notion of BCK/BCT-algebras were interoduced by Iseki [3], and were extensively in- 
vestigated by many researchers. They are two important classes of logical algebras. The notion 
of fuzzy sets was introduced by Zadeh [5] was applied to BCI-algwbras by Xi [4]. In this paper, 
we discuss the notion of fuzzy sub-halfalgebra and we also recall the notion of normalized fuzzy 


extension and maximal fuzzy extension and obtain a relation between these two concepts. 


§2. Fuzzy extension 


By a BCl-algebra we mean an algebra (X,*,0) of type (2,0), satisfying the following 
conditions : 

(I) ((w«y) «(a 2)) *(2¥y) =0, 

(II) (v (w#y)) «y = 0, 

(III) «*xz=0, 

(IV) «xy=Oand yxx=0 imply c=y, 
for all x, y,z,E X. 

We can define a partial ordering “<<” on X by x < y if and only if «*y = 0. Throught 
this paper X will always mean a BCI-algebra unless otherwise satisfied. 

Definition 2.1. Let X be a BCl-algebra. A map yp : X — [0,1] is said to be a fuzzy 


sub-halfalgebra of X if a * y < z implies u(z) > min {p(x), w(y)} for all x,y,z € X. 
1 1 it 1 
20°91" 92° "on 


Example 2.2. Let X = { , under usual Subtraction. 
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Define uw : X — [0,1] by u(x) = 1-2 for every x € X. Now it is easy to verify that pis a 
fuzzy sub-halfalgebra of X. 

Theorem 2.3. Let jz be a fuzzy sub-halfalgebra of X and x,y,z € X. Then p(z) > 
min {pi(x), u(y)} dose not in general imply x *y < z in X. 


1 1 
Proof. In the above example we see p(z) > min {uu(x), u(y)}, but y A z for x = 507 ¥ = or 


and z= —. 
Definition 2.4. Let w and v be any two fuzzy sub-halfalgebras of X such that 
(i) p(a) < v(x) for every x € X. 
(ii) ifax*xy < z and p(z) = min {p(x), u(y)} then v(z) = min {v(x),v(y)}. Then we say 
that v is a fuzzy extension of p. 


1 11 1 
Example 2.5. Let X = { 


90791? 927 7 On 


\ be a BCl-algebra with respect to usual 


Subtraction. 

Define w: X — [0,1] by u(a) = 1-—« for every x € X and v(x) = 1— 5 for every x © X 
and it is easy to verify that v is a fuzzy extension of p. 

Theorem 2.6. Let Let 4 and v be two fuzzy sub-halfalgebra of X. If v is a fuzzy extension 
of wand «xy < zin X with v(z) = min{v(x),v(y)} then it need not in general imply 
u(z) = min {u(x), u(y) }- 


Proof. The proof is by an example. Take X and yz as in Example 2.5. Define : X — [0,1] 
1 1 


by v(x) = 1 for every x € X. Now if we take 7 = = and y = x then x *y = a Pe It 
is easy to verify that v(z) =1, min {v(x),v(y)} =1. u(z) = 2 and min {y(x), u(y)} = 0. This 
prove that v(z) = min {y(x),v(y)} and p(z) ~ min {py (2), ay This completes the proof. 

Definition 2.7. Let and v be any two fuzzy sub-halfalgebra of X. Then the fuzzy 
extension v (of) is said to be an open fuzzy extention of w if x*y < z in X with p(z) > 
min {(x), (y)} implies v(z) > min {v(x), v(y)}, v is said to be closed fuzzy extension of y if 
it is not an open fuzzy extension of ju. 

Theorem 2.8. Every open fuzzy extension is a fuzzy extension but not conversely. 

Proof. Every open fuzzy extension is a fuzzy extension, directly follows from the definitions 
of open fuzzy extension and fuzzy extension. Converse is not true.This is explicit by the 
following example. 

Example 2.9. Let X be any BCI-algebra and ys be any fuzzy sub-halfalgebra of X. With 
|Im(,)| > 2. Then it is easy verify that vy : X — [0,1] given by v(x) = 1 for every x € X isa 
fuzzy extension of pu. It is easily seen v is not an open fuzzy extension of p. 

Definition 2.10. Let X be a BCL-algebra. A countable collection of fuzzy sub-halfalgebras 
of X, denoted by {y;|¢ = 0,1,2,--- } or {pu;} lim®2, is called fuzzy extension chain for X if wi+1 
is a fuzy extension of yu; for i = 0,1,2,---. A fuzzy extension chain {j:;}limj2, is said to 


generate X if UJ uw; = lx (Here 1x denotes the map yw: X — [0,1] such that u(x) = 1 for 
i=0 
every x € X). A fuzzy extension chain {y;}lim72, is said to be a fuzzy extension chain for X 
if f4i41 is an open fuzzy extension of yu; for i =0,1,2,---. 
Example 2.11. Let X be any BCl-algebra and py, : X — [0,1] be defined by: puz(x) = 


pa for every k = 0,1,2,---. Then we get a countable collection of fuzzy sub-halfalgebra 
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{ux} limp, of X. If oxy < z and pegi(z) = min {pp41(%), Ux+i(y)} then we have puz(z) = 
min {14%(x), px(y)}. Hence x41 is a fuzzy extension of juz, for all k = 0,1,2,---. This establishes 
that {u,} lim? is a fuzzy extension chain for X. Now we will prove {ju, } lim72.) of X generates 
the X. Let x be an arbitrary element of the X then we have: 


(Um) = sup {ux(2)}= sup {par} =t= 1x) 


Pe k=0,1,2,- k=0,1,2,-~ 


That is (Ulim?? 9 ux)(%) = 1 x(a) for every « © X. Thus we have proved that (LU lim? 9 ux) (x) = 
lx. This implies that {u,}lim7f2> generates the X. 

Theorem 2.12. Let js and v be any two fuzzy sub-halfalgebra of X. If the fuzzy ex- 
tension v(of w ) is an fuzzy extension of yw then pu(z) = min {p(ux), u(y)} if and only if 
y(2) = min {v(2), v(y)}. 

Proof. We prove the converse of the contra positive method. Let x *y < z in X and 
p(z) A min {u(x), u(y)} to prove v(z) £ min {v(x), v(y)}, since ps is a fuzzy sub-halfalgebra of 
X, then we have p(z) > min {u(x), u(y)}. Since v is an open fuzzy extension of 4. This implies 
that v(z) > min {v(x),v(y)}. That is v(z) ¢ min {v(x), v(y)}. 

Hence we have proved that if «*y < z in X and p(z) A min {p(x), u(y)} then v(z) F 
min {v(x),v(y)}. This proves that if the fuzzy extension v (4) is an open fuzzy extension ju 
then u(z) = min {p1(x), (y)} if and only if v(z) = min {v (x), v(y)} for every xx y < zin X. 


§3. Antifuzzy extension 


Now we proceed on to give the definition of antifuzzy extension. 

Definition 3.1 A fuzzy subset ys of X is said to be antifuzzy sub-halfalgebra of X if 
x*xy <zin X implies p(z) < max {p(z), w(y)}. 

Definition 3.2. Let yw and v be two antifuzzy sub-halfalgebras of X. Then v is said to be 
an antifuzzy extension of yu if the following two condition hold: 

(i)u(x) > v(x) for every x € X, 

(ii) if a xy < z and p(z) = max {p(x), u(y)} then v(z) = max{v(x),v(y)}. As in case 
of fuzzy extension we can prove that if v is an antifuzzy extension of w and x *y < z and 
v(z) = max {v(x),v(y)} then it need not in general imply that p(z) = max {u(x), u(y) }. 

Definition 3.3. Let y and v be two antifuzzy sub-halfalgebras of X. Then the antifuzzy 
extention v ( of yu ) is said to be an open antifuzzy extension of pu if: «* y < z and p(z) < 
max {pi(x), u(y)} then v(z) < max {v(x),v(y)}. If vy is not an open antifuzzy extension of iu 
then we say that v is a closed antifuzzy extension of pu. 

Definition 3.4. Let X bea BCI-algebra. A countable collection of antifuzzy sub-halfalgebra 
of X , denoted by denoted by {ju;|¢ = 0,1,2,--- } or {j0;} limP2y is called an antifuzzy extension 
chain for X if w;41 is an antifuzzy extension of yw; for 1 =0,1,2,---. 


1 1 1 
Example 3.5. Let P = { \ take usual subtraction in X. Define pp: X —- 


(0, 1] by u(x) = x for every x € X. It is clear that yu is an antifuzzy sub-halfalgebra of X. 
Define p , uv: X — [0,1] by: w(x) = & for all x € X v(x) = 5 for all ¢ € X. Then it is 


easy to verify that v is an antifuzzy extension of pu. 
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Theorem 3.6 Let X be a BCI-algebra. Then yu is a fuzzy sub-halfalgebra of X if and only 
if w° is an antifuzzy sub-halfalgebra of X. 

Proof. Let s be a fuzzy sub-halfalgebra of X and x * y < z in X. Since p(x) = 1 — p(x) 
for all x € X, we nave pu(z) > min {p(x), u(y)} if and only if w°(z) < max {p°(x), u°(y)}. This 
prove that py is a fuzzy sub-halfalgebra of X if and only if u° is an antifuzzy sub-halfalgebra of 
Xx. 

Theorem 3.7. Let yz and v be two fuzzy sub-halfalgebra of X. Then: 

(i) v is a fuzzy extension of yz if and only if v° is an antifuzzy extention of p°. 

(ii) v an open fuzzy extension of yz if and only if v° is an open antifuzzy extension of °. 

Proof. Let v be a fuzzy extension of j: to prove Vv“ is an antifuzzy extension of u°. Since v is 
a fuzzy extension of jz, we have p(x) < v(x) for all x € X that is u°(x) > v°(«) for alla € X. Let 
vey < zin X and p°(z) = max {y°(a), w°(y)} that is 1—p(z) = max {1 — p(x), 1 — p(y) }. Then 
1— p(z) = 1— min {pu (x), u(y)}. Hence u(z) = min {p(x), w(y)}. Since v is a fuzzy extension 
of pp and a *y < z in X, we have v(z) = min {v(x),v(y)}. Using the properties of Min and 
Max function we have 1 —v(z) = max {1 — v(x), 1 — v(y)}. That is v°(z) = max {v“(z), v“(y)}. 
Thus v° is an antifuzzy extension of j°. 

(ii) Let v be an open fuzzy extension of yw. To prove v~ is an open antifuzzy extension 
of w°. Since v is an open fuzzy extension of y we have p(x) < v(x) for all « € X that 
is u°(x) > v(x) for alla € X. Let xxy < zin X and p°(z) < max{p°(x), u°(y)} then 
1 — p(z) < max {1 — p(x),1— p(y)}. This implies 1 — p(z) < 1 — min {p(x), u(y)}. That is 
p(z) > min {p(x), u(y)}. As v is an open fuzzy extension of uw and « *y < z. We have v(z) > 
min {v(x),v(y)} then 1 — u(z) < max {1 — p(x), 1 — u(y)} that is v°(z) < max {v*(x),v“(y)}. 
Thus v° is an open antifuzzy extention of p°. Conversely let v° be an open extension of pu° then 
we prove that v is an open fuzzy extension of jz using the fact 1—max {a,b} = min {1 — a,1 — b} 


for all a,b € [0,1]. 


T«d 1 
Example 3.8. Choose p and v the example3.5 where P = {a Tee xh , p(x) = 1-2 


and v(a#) = 1— : for all  € X. It is easy to verify that v is a fuzzy extension of u. Now we 


calculate ° and Vv as p°(x) = & and v(x) = as for every c € X. Clearly v is an antifuzzy 
extension of u°. This shows that v is a fuzzy extention of w if and only if v° is an antifuzzy 


extension of p°. 


84. Maximal fuzzy extension 


Now we proceed on to define the concept of maximal fuzzy extension. 

Definition 4.1 Let ys and v be any two fuzzy sub-halfalgebra of X. A fuzzy extension 
v(of 1) is said to be maximal fuzzy extention of y if there dos not exist any fuzzy sub-halfalgebra 
f of X such that v C f. 

Example 4.2 Let P = {1,2,3,--- ,n} be a BCl-algebra with respect to usual subtraction. 
Define p(a) = é and v(x) = 1 for all x € X. It is easy to verify that both y and v are fuzzy 
sub-halfalgebra of X. Now we cannot find any fuzzy sub-halfalgebra f of X such that v C f. 


Hence v is the maximal fuzzy extension of w. 
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Definition 4.3 Let yz and v be any two fuzzy sub-halfalgebra of X. A fuzzy extension 
v(of 1) is said to be a normalized fuzzy extension of ps if v(x) = 1 for some x € X. In the above 
example the fuzzy sub-halfalgebrs v is a normalized fuzzy extension of ju. 

Definition 4.4. Let ys be a fuzzy subset of X and a € [0,1 — sup {u(x) : « € X}]. A map 
ut (ax) : X — [0,1] is called a fuzzy translation of w if w2(x) = u(x) + a for alla € X. 

Example 4.5. Let X = {2,4,--- ,2n}. Define : X — [0,1] by u(x) = . for alla € X. 


1 1 1 
Then for a = — we have po (a) = — + Z for all z € X. It is easy to verify that uZ is a fuzzy 
x 


translation of yp. 

Theorem 4.6. Let wu be a fuzzy sub-halfalgebra of X and a € [0,1 — sup {p(x): a € X}| 
then every fuzzy translation 2 of ys is a fuzzy sub-halfalgebra of X. 

Proof. Let u be a fuzzy sub-halfalgebra of X and a € [0,1 — sup {p(x):a€ X}]. For 
“xy <zin X , we have u(z) > min {py(x), p(uy)} that is a+ u(z) > a+ min {p(x), u(y)} that 
is a+ u(z) > min {a+ pu(x),a + p(y)} by the definition of fuzzy translation 
wi (z) > min {u2 (x), u2(y)}. That is uZ is a fuzzy sub-halphalgebra of X. 

Theorem 4.7. Let js be a fuzzy sub-halfalgebra of X and a € [0,1 — sup {u(a) : « © X}). 
Then every fuzzy translation 2 (of) is a fuzzy extension of pu. 

Proof. Let jz be a fuzzy sub-halfalgebra of X and a € [0,1—sup {u(a) : « € X}], we know 
by Theorem4.6 ju2 is a fuzzy sub-halfalgebra of X. Clearly w2(x) > (zx) for all « € X. Further 
ifaxy <zin X , and p(uz) = min {p(z), w(y)} then w2(z) = min {u2 (x), uo (y)}. Hence yt 


is a fuzzy extension of pL. 


1 1 1 
3073n xh be a BCLalgebra with respect to usual sub- 
traction. Define uw: X — [0,1] by w(a#) = 1 —« for all e © X. Now it can be easily checked 


Example 4.8. Let X = 


that ys is a fuzzy sub-halfalgebra of X. For a = 3 we have wi(r) =1—a2+ : for all x € X. It 
is easy to verify that pZ is a fuzzy extension of p. 

Theorem 4.9. Let js be a fuzzy sub-halfalgebra of X and a € [0,1 — sup {u(a) : « © X}). 
If u2 is fuzzy translation of pz then: 

(i) wt is a maximal fuzzy extension of ju if and only if p is a constant map. 

(ii) If has sup property then 2 is normalized fuzzy extension of ys but not conversly. 

Proof. (i) Let 4 be a fuzzy sub-halfalgebra of X and a € [0,1—sup{pu(a) : 7 € X}]. Then 
by Theorem 4.7 y7 is a fuzzy extension of yp. Let w“T, be a maximal fuzzy extension of y then 
w(x) = 1 for all x € X , that is u(uur) = sup {p(y)|y € X} for every x € X. Hence pisa 
constant map. Conversly let 2 be a constant map, that is u(a#) = @ for all x € X ( where @ is 
a fixed element of [0, 1]). 

Now consider the fuzzy translation yw? of pw, u2 (x) = w(ur) +a = B+1-—6 =1 for all 
x € X. Hence pl is a maximal fuzzy extension of ys. (ii) Let ps have sup property, that is for 
every subset T of X there exists to € T such that (to) = sup {u(t)|t €¢ T}. If possible let as 
assume that p2 is not a normalized fuzzy extension of ys then 2 (x) <1 for all e € X. That 
is p(x) < sup {u(y)|y € X} for all « € X. Since yz has supproperty for every subset T of X 
there exists t) € T such that u(to) = sup {u(t)|t € T}. So we get u(x) < p(to) which is a 
contradiction. 


Hence pZ is a normalized fuzzy extension of py. If 2 is a normalized fuzzy extension of ps 
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then it does not in general imply that y has sup property. This can be seen from the example 
given below. 

Consider X = [0,1] under usual subtraction. Define «: X — [0,1] by u(x) = 1-2 for all 
x € X. It is easy to see that wis a fuzzy sub-halfalgebra of X. Let a € [0, 1—sup{u(ax) : 7 € X}] 
then pi is a normalized fuzzy extension of jz for there exists 0 € X such that 7(0) = 1. Take 
T = [0,1] then T C [0,1] and sup {p:(t)|t € T} = 1. Now it is easy to verify that there is no 
to € T such that (to) = 1. Hence ys dosenot have sup property. 

Definition 4.10. A fuzzy subset jz of a set X has the weak sup property if there exist 
xo € X such that p(ao) = sup {u(a)|a2 € X}. 

Theorem 4.11 Let u be a fuzzy subset of X. If has sup property then yz has weak sup 
property but not conversly. 

Proof. Let w be a fuzzy set of X. If has sup property then by the definition of sup 
property for every subset T of X there exist to € T such that p(to) = sup {u(t)|t € T}. The 
above condition is true for every subset T of X. If we take T as X then there exists tp € T 
such that pu(to) = sup {p(x)|a € X}. 

Hence yp has weak sup property.However the converse is not true this can be seen by the 
following example. Consider the set [a, 3] where a and (@ are any two arbitary fixed numberes 
in interval [0,1] with a < @. Define p: [a, 8] — [0,1] by w(x) = @ for alla ec X. 

Then pu(3) = sup {u(x)|2 € X} but for T = [a, 3) we have sup {u(x)|z € T} = @ and there 
is no tp) € T such that p(to) = G. That is w has weak sup property but ~ dose not have sup 
property. 

Theorem 4.12. Let ys be a fuzzy sub-halfalgebra of X and a € [0,1 —sup {u(x) : 7 © X}). 
Then pu is a normalized fuzzy extension of yw if and only if 4 has weak supproperty. 

Proof. Let ys be a fuzzy sub-halfalgebra of X and a € [0,1 —sup{pu(x): a2 € X}]. Let 
ut be a normalized fuzzy extension of jz hen by he definition of normalized fuzzy extension we 
have 2 (xo) = 1 for some ro € X. That is (xo) = sup {u(x)|x € X} for some xp € X. Hence 
pt has weak sup property. Conversly let js have weak supproperty. 

Then there exists zo € X such that u(x) = sup {u(x)|a € X}. Now consider the fuzzy 
translation w2 of wp. wl(x) = u(x) +1—sup{pu(2)|z7 eX}. Take x = xo, then we have 
uz (ao) = w(xo) +1 — u(x) = 1. Hence pZ is normalized fuzzy extension of u. This completes 
the proof of the theorem. 

Definition 4.13. Let jz be a fuzzy subset of X and 6 € [0,1]. A map pg: X — [0,1] is 
called a fuzzy multiplication of pw if ue(x) = 6.u(x) for all x € X. 

Theorem 4.14. Let 4s be a fuzzy sub-halfalgebra of X and a € [0,1 —sup {pu(x) : a € X}| 
then every fuzzy translation 7 (of) is a fuzzy extension of fuzzy multiplication pg(of 1). 

Proof. Let jz be a fuzzy sub-halfalgebra of X. If a € [0,1 —sup {u(x) : x € X}] then fuzzy 

translation w2 and the fuzzy multiplication jg are fuzzy sub-halfalgebras of X. We see that 
for all x € X, ug = B.u(x) < u(x) < at p(x) = pT (2). 
Also if cxy < z and ug(z) = min {4¢(z), we(y)} then this implies G.u(z) = @. min {4g(2), a(y)} 
that is u(z) = min {y(x), u(y)} that is a+ u(z) = min{a+p(x),a+ u(y)} that is w2(z) = 
min {2 (x), u2(y)}. Hence pZ is a fuzzy extension of jug for all a € [0,1—sup {u(x) : « € X}). 
Hence the theorem. It is to be noted that for G = 0 the result is not true. 
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Theorem 4.15. Let ys be a fuzzy sub-halfalgebra of X and a € [0,1—sup {u(x) : « © X}). 
hen {uf} limg2.g where (S = {y:0 <7 < a}, and ¥ rational) is an open fuzzy extension chain 
for X. Further { uz } limgeg generates the BCI-subhalf algebra X if and only if pz is a constant 
map. 

Proof. It is easy to see from the proof of the Theorem 4.14 k € S and 7 is a fuzzy 
sub-halfalgebra of X. Choose ki, ko € S such that 0 < ky < kg <k then be (a) < Mio (a) for 
all x € X. Thus for kj,kj € S(t < j), ue. (2) < Mh, (x) for all 2 € X without loos of generality 
we have 0 < ky < ky < +++ Shy < hig Sk. 

Then by the above construction uj, (7) < Hi, (2) Ss Sf (2) < Mh, (2) <s1+ < pT (a) 


for all a € X. Now ifaxy < z and pi (z) = min {pg (x), we (y)}, then it is easy to see 


that Mie (2) = min {ub.,.(0),uk,,()} for i = 0,1,2,---. Further if a *y < z and pz (z) > 
min { Lg, (x), wR, (y)}, then we can pins, (2) > min {uh 4 (2), Heyy (y) }for i =0,1,2,---. Hence 


{ un} limgeg is an open fuzzy extension chain for X. 

Now we will prove that { tn } limgeg generate the BClI-algebra X if and only if p is a 
constant map. Let {uz} lim?2 generate the BCI-algebra X. Then we have Ulim@e.5 wy = 1x 
that is sup {ug (x)|k € S} = 1 for every x € X, that is sup {u(x) + klk € S} =1 for alla e X, 
that is u(x) +1 —sup {u(y)|y € X} = 1 for all « € X. That is u(x) = sup {u(y)|y € X} for all 
x € X. This proves that ys is a constant map. 

Conversely let 4. be a constant map that is (a) = y for all « € X (where 7 is fixed elements 
in (0, 1]). Consider (Ulimf2g wz) (x) for an arbitrary x € X. 


(lim wf)(2) = sup {uf (x)|k € S} = sup {n(z) + blk € 5} 
= p(x) +1—sup{u(ylye X}=y+1l-y=l=l1xa) 


since x is an arbitary element of X we have (Ulimg2g ug, )(2) = 1x 2) for all x € X. Hence 
{uz } lime g generates the BCLalgebra X. 
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Abstract Let D # 1 be a positive non-square integer. Let a be a positive integer. In this 


paper is considered the Pell’s equation «? — Dy? = +a and some recurrence relations are 
obtained. 
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§1. Introduction 


Let D #1 be any positive non-square integer and a be any fixed integer. The equation 


x? — Dy? = +a 


is known as Pell’s equation. It is named mistakenly after John Pell (1611-1685) who was a 
mathematician who in fact did not contribute for solving it (see [8]). 


For a = 1, the Pell’s equation x? — Dy? = +1 is known as classical Pell’s equation and it 
has infinitely many solutions (#,,y,) for n > 1. There are different methods for finding the 
first non-trivial (21, y1) solution called the fundamental solution from which all others are easily 
computed (see [3], and [9]). 

There are many papers in which different types of Pell’s equation are considered (see [1], 
[2], [4], [5], [6], [7]). 


In his paper [1] A. Tekcan considered the equation x? — Dy? = +4 and obtained some 


formulas for that equation. 


In this paper we will consider a more general equation, indeed x? — Dy? = +a, with a 
positive integer. Some similar results to those in [1] and [2] will be obtained. The ideas used in 
[1] and [2] are going to be used again in order to prove the main results. 


§2. The pell’s equation 2? — Dy? =a 


Theorem 1.1. Let (x1, y1) be the fundamental solution of the Pell’s equation x? — Dy? = a 
and let 


n 


Un xy Dy 1 


= (1) 
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for n > 1. Then the integer solutions of the Pell’s equation x? — Dy? = +a are given by 


Cor) (4 —-) (2) 


Proof. We prove the theorem using the method of mathematical induction on n. 


Clearly it is true for n = 1. We assume that x? — Dy? = ais satisfied for (ap_1, Yn—1), ie. 


2 2 2 2 
Un—1 Dvr_4 uy, — Dun 
an a Dyn, = i 5 =— A= w=. (3) 
(Ja “)? (va *)? a 


From (3) we have x and y should have index n (See the original) u2_,Dv2_, =a"~1. 


Then we show that x2 — Dy? = +a. 
First we note that 


n 


Un Xy1 Dy 1 
Un Y1 Ly 0 
n-1 
7 ry Dy zy Dy 1 
Yr 2 Yr @y 0 
_ x, Dy Un-1 
Y1 Ty Un-1 


tUn—1 + DyUn-1 


YiUn—1 1 L1Un-1 
So 


Un = ©1Un—1 + Dy vn—-1 


Un = Y1Un—-1 + £1Un-1 
Then 


Dy, 
u Dy 
n—1 
va) 


( 
(@1Un—1 + Dy¥n—1)* — D(yrun—1 + 21Un—1)? 
qr-1 


x 


yore 22,2 2,2 2,,2 
Ujus 1 + 241Un—-1Dyitn—-1 + D*yjus_, — D(yfus_1 + 2y1un—-1%1Un—1 + L{VZ_1) 


qn-l 
_ (et — Dyf)un_s — Dun_s(@i — Dyi) 
qn-1 
_ (xt = Dyj)(un_a — Dut) 
_ qr-l 
aa} 
= qn-l 
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completing the proof. 

From the theorem 1.1, the following corollary can be proved. 

Corollary 1.2. If (x1, y,) is the fundamental solution of the Pell’s equation 2? — Dy? = a, 
then 


2, = Tent + Dyiyn-1 | — Yitn—1 + L1Yn—-1 (5) 
ni ’ nm 
Ja Ja 

and 

In Ln— 

7 a = —Jay, 

Un Yn-1 

for n > 2. 


Proof. By (4) Un = @1Un—-1 + Dyivn-1;. Un = Y1Un—-1 + L1Un-1. 
By (2) un = va ‘Ini Va. = Ja” “Yn, 80 then un_1 = Va”? Typos Uncy = 
Ja" *Un-1- 
Hence we get: 
va" * ‘In = 21° Ja" ‘Ln + Dy Ja" * -Yn—-1 
= Va" (e1tn-19n + Dy yn-1)- 
L12n-1 + DY Yn-1 


Ja 
In a similar way, one obtains: 
-1 -2 —n—2 
Va yn = WJVO Spt ea" + Yn—1 
-2 
= Ja” (Yi2@n-1 + L1Yn—1)- 


Yr@n—1 + L1Yn-1 
It follows that = , 
Yn va 


It follows that x, = 


Finally, 


Tn En-1 
= @LnYn-1 — Yntn-1 
Yn Yn-1 
L1Xn—-1 + Dyyn-1 Yitn—1 + L1Yn-1 : 
Va Yn-1 Va n-1 
ee Zz Dy;-1) 


Ja 


completing the proof. 
Remark. For a = 4, (2) and (5) represent the results of [1]. 
Now following the same argument of proof as in [2] we can prove the following theorem: 
Theorem 1.3. If (x1, y1) is the fundamental solution of the Pell’s equation x? — Dy? =a 
then (Zn, Yn) satisfy the following recurrence relations 


Ln = (41 —1)(an 1T&y 2) — Ln 3 


(6) 


Yn = (41 — 1)(Yn-1 + Yn—2) — Yn—3 
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for n > 4. 

Proof. It is a matter of tedious calculations to show that relations (6) are true for 
n = 4.(See [2]) 

Next, we assume that (6) holds for n and we show that it holds for n+ 1. 

Indeed, by (5) and by hypothesis we have: 


£1((1 — 1)(n—1 + Ln—2) — En—3) + Dyr((21 — 1)(Yn-1 + Yn—2) — Yn—3) 
— Va 
Di (Cap Fy 9) DY (Yet + Uno) L1%n—-3 + Dy yn-3 
(a=) ie a 
(x 1) T1Ln—-1 + DyiYn-1 A T1Ln—-2+ DyiYn—2 e 
i Va Va = 


(x1 a 1)(2n a ti) — In-2, 


y = yi (21 _ 1) (an 1 t In 2) In 3) t xi((44 ~~ 1)(Yn 1 + Un 2) — Un 3) 
n+1 Va 
G 1) Y1(Gn—1 + Fn—2) + L1(Yn—1 + Yn—2) Y1Ln—3 + L1Yn-3 
: Ja Ja 
(2, —1) YrTn—-1 + 21Yn-1 , YiTn-2 + L1Yn—-2 
1 Ja T Va Yn-2 


= (1 —1)(Yn + Yn—1) — Yn-2; 


completing the proof. 


§3. The negative Pell’s equation x? — Dy” = —a 


Theorem 2.1. If (x1,y1) is the fundamental solutions of the negative Pell’s equation 


x? — Dy? = —a, then the other solution are (%2n41, Yon+1) Where 
U2n+1 V2n+1 
(Tan41; Yont1) = (4, “nit) ’ (7) 
for n> 0. 
Proof. Clearly it is true for n = 0. 
We assume that x? — Dy? = —a holds for (a2n—1, yen—1) 80 
2 2 
Wan—1 ~ DvIn-1 
Woe _ DySn—1 _ = qen-2 = = —a. (8) 
From (8), we have u3,,_; — Dv3,_, = —a?"—1. 


Then we show that 73,4, — Dy3,41 = —a. 
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and 


First we note that 


So 


2n+1 
U2n+41 _ ry Dy 1 
Vont+1 Yo Ly 0 
2 
- x, Dy x, Dy 
Y1 Beal Y1 xy 
2 
_ xz, Dy U2n—1 
Y1 vy V2n-1 


2n-1 


Uon—1(27 + Dy?) + 2Daiyiven—1 


20141 Uan—1 + (27 + Dy?)ven-1 


U2n4+1 = Uon—1(22 + Dy?) + 2D21Y1V2n— 


Vandi = 2 yi Uan—1 + (a7 + Dy?)ven—-1 


Using (9) we obtain: 


2 2 
Lint — DYan41 


2 2 
Vang — Dvon41 
q2” 


1 


(uan—1(24 + Dy?) + 2D21y1V2n—1)? - D(221y1U2n-1 + (x? + Dy?)ven—1)" 


q2” 


Udn—1 (21 + Dyt)* — Dv3n—1(@i + Dyj)? + 4D? xi y{V3n—1 — ADUZn TT YT 


q2n 


(a7 + Dy{)? (U3n—1 7 Dv5n—1) — AD xii (Udn—1 — Dv3n_1) 


q2n 
(U5n—1 = D3n—1) (24 7 Dyj)? 
a2” 
—q2"-1 a2 
q2n 
—a. 


From the theorem 2.1, the following corollary can be proved. 


Corollary 2.2. If (x1, y,) is the fundamental solution of the Pell’s equation x? — Dy? = —a 
then: 


(2? + Dy?)zen—1 + 2D2r1y1Yen—1 


(10) 


L2n+1 = , 
a 
2 2 
221Y1L2n-1 + (aq + Dy?)yon-1 
Y2n41 = 
a 
TIn4+1 L2n-1 
= 2211 
Yant+i Y2n-1 
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forn> 1. 
Proof. By (9) 
U2n+1 = Urn —1(X4 + Dy?) + 2D2141Van-1, 
Vont1 = 2r1y1Uan—1 + (27 + Dy?) ven-1. 
By (7) 
U2n+1 = a” Lon+1, V2n4+1 = O”'Yon+1; 

so then 

ten-1 = 0" “@en-1, Van—1 =a" “Yon-1- 


Hence we get 
a” fon41 => alan me Coy + Dy?) + 2Dayy,a"—1yon—1 
= a" (@an—1(a} + Dy) + 2Daryiyen—i). 


It follows that 
(27 + Dy?) @on-1 + 2Daiy1yen-1 


L2n+1 = i 
In a similar way we obtain 
= 2 2) n—-1 
a"Yont1 = 2a’ ton-1%i1y1 + (27 + Dy{ a” yon-1 


= a™*(2ton—121y1 + (xt + Dy?) y2n—1)- 


It follows that ; r 
221 Y1Lan—1 + (27 + Dy{)yon-1 


Y2n+1 = 
a 
Finally, 

T2n4+1 V2n-1 

Y2nt+l Y2n-1 
= £9n41Y2n-1 — Y2n4+1%2n-1 
(at + Dy?)?tan-1 + 2aryry2n-1 20141 on—1 + (a? + Dy?) yan-1 
= a Y2n-1 re 2n-1 
_ 2De1YiY9n—1 ~ 2L1Y1T3n—1 

a 
= —201y1(@3n-1 — DySn—1) 
a 

= 20191; 


completing the proof. 
Remark. (7) and (10) represent a general case of the results of [1]. 
Now following the same argument of proof as in [2] we can prove the following theorem: 
Theorem 1.6. If (1, y1) is the fundamental solution of the Pell’s equation x? — Dy? = —a, 
then (Z2n+41, Yan+41) satisfy the following recurrence relations 


Tanti = (27 + 1)(Lan—1 + Lan—3) — Lon—5 


(11) 


Yan+1 = (x? a 1)(Yyan—1 le Yon—3) — Y2n-5 
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for n> 3. 

Proof. Tedious calculations, yield that relations (11) are true for n = 3,(see [2]). 

Now we assume that (%n+41,Yan+1) satisfies(11) and we prove that also (%an43, Yan+3) 
satisfies(11). 

Indeed by (10) and hypothesis we obtain. 


L2n+3 
(at + Dy?) (a7 + 1) (@en-1 + Fan—3) — Cons) 
a 
re 2Daxyyi((z7 + 1)(yan—1 + Yon—3) — Y2n—s) 
a 
_ 2 (a? + Dy?)(tan—1 + Zan—3) + 2Dx1y1(Yon—1 + Yon—3) 
(xy + 1) 
17 
a 
(a? + Dy?) ®on—5 + 2Dr1y1Yon—5 
Ga 
= @ aij (af + Dy{)ton—1 + 2Dr1y1Y2n-1 de (aj + Dyj)ton—3 + 2Dariy1yon—3 
1 ; 7 
— L2n-3 
= (a7 +1)(en41 + Zan—1) — Len—3, 
and 
Y2n4+3 
_ Qaryi((27 + 1)(ton—1 + Pon—3) — Lan—s) + (ej + DyZ)((@7 + 1)(Yon—1 + Yon—s) 
a 
= (241) (Sunes + ®on—3) + (27 + Dy?) (Yon-1 + ea} 
- 1 
a 
20141 Fon—5 + (27 + Dy? )yon—s 
a 
= (#241) 221y1%an—1 + (a7 + Dy?)yen-1 _ 281Y1L2n—3 + (27 + Dy?) yon—3 
‘ll a T rm 
— Y2n-3 
= (x} ar 1)(Yanti = Yon-1) — Y2n-3; 


completing the proof. 
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Abstract For any positive integer n, the F.Smarandache LCM function S'L(n) is defined as 
the smallest positive integer k such that n | [1,2,--- ,k], where [1,2,--- ,k] denotes the least 
common multiple of 1,2,--- ,&, and let n = pf{'p5? ---p$* be the factorization of n into prime 
powers, then Q(n) = arp: + azp2 +-:-+asps. The main purpose of this paper is using the 
elementary methods to study the mean value properties of Q(n)SL(n), and give a sharper 


asymptotic formula for it. 
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§1. Introduction and Results 


For any positive integer n, the famous F.Smarandache LCM function SZ(n) defined as 
the smallest positive integer k such that n | [1, 2, ---, k], where [1, 2, ---, k] denotes 
the least common multiple of 1, 2, ---, k. For example, the first few values of SL(n) are 
SL(1) = 1, SL(2) = 2, SL(3) = 3, SL(4) = 4, SL(5) = 5, SL(6) = 3, SL(7) = 7, SL(8) = 4, 
SL(9) = 6, SL(10) =5, SL(11) = 11, SL(12) = 4, SL(13) = 13, SL(14) =7, SL(15) =5, ---. 
About the elementary properties of SL(n), some authors had studied it, and obtained some 
interesting results, see reference [2] and [3]. For example, Lv Zhongtian [4] studied the mean 
value properties of SL(n), and proved that for any fixed positive integer & and any real number 


xz > 1, we have the asymptotic formula 


mw E bj - x? x? 
SL(n) = —-— O | ——— 
y (n) 12 Inz + d. a + (a5 -| 2 


nN<au 


where b; (¢ = 2,3,--- ,&) are computable constants. 
On the other hand, Chen Jianbin [5] studied the value distribution properties of SL(n), 


and proved that for any real number x > 1, we have the asymptotic formula 


sum) —Pemyr=2-¢(2) Z +0 ( e | 


n<ux 


where ¢(s) is the Riemann zeta-function, and P(n) denotes the largest prime divisor of n. 
Now we define a new arithmetical function Q(n) as follows: Q(1) = 0; for n > 1, let n = 


pi ps? +++ p&be the factorization of n into prime powers, then Q(n) = api + agp2 +--+ sDs. 
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Obviously, for any positive integers m and n, we have Q(mn) = Q(m) + O(n). That is, Q(n) 
is the additive function. The main purpose of this paper is using the elementary methods to 
study the mean value properties of Q(n)SL(n), and give a sharper asymptotic formula for it. 
That is, we shall prove the following conclusion: 


Theorem. For any real number zx > 1, we have the asymptotic formula 


5 B dix? 2° 
Q(n)SL(n) = “+0 ; 
2 OE = ir (=a) 


n<u 


where d; («= 1, 2, --- , &) are computable constants. 


§2. Proof of the theorem 


In this section, we shall use the elementary methods to complete the proof of the theorem. 
In fact, for any positive integer n > 1, let n = p{'ps? ---p®* be the factorization of n into 


prime powers, then from [2] we know that 


SL(n) = max{p{",p5?,-°- Sate (1) 


and we easily to know that 


Q(n) = ayp; + agp2 +--+ + asps. (2) 


Now we consider the summation 


S> Q(n)SL(n). (3) 


n<u 


We separate all integer n in the interval [1,2] into four subsets A, B, C and D as follows: 

A: p> J/nandn=™m-p; 

B: n3 <p, < po < Vn and n=™m- py: po, where p; (¢ = 1,2) are primes; 

Ci ns <p<VJnandn=m-p?; 

D: otherwise. 

It is clear that if n € A, then from (1) we know that SLZ(n) = p, and from (2) we know 
that Q(n) = Q(m) +p. Therefore, by the Abel’s summation formula (See Theorem 4.2 of [6]) 
and the Prime Theorem (See Theorem 3.2 of [7]): 


where a; (¢ = 1,2,...,&) are computable constants and a, = 1. 
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We have 
S7Q(n)sL(n) = SY Q(m)+p)p= So w+ SS (A(m)p) 
neEA mp<x mpcx mpcx 
m<p m<p m<p 
= Yd +O’) 
M<J/xzM<PS 7 
2 z. 
_ Ee ce 2 - 2 
miVz 
1 x? a byx x? 
= «~C(3 : O 4 
36( Noe 2h (<ar-). (4) 
where ¢(s) is the Riemann zeta-function, and 6; (¢ = 2,3,--- ,&) are computable constants. 
Similarly, ifn € B, then we have SL(n) = pz and Q(n) = Q(m) + pi + po. So 
S Q(n)SL(n) = > (Q(m) + pi + pz) po = yy p, +O oe Pip2 
neB mpip2<@ mpip2<x mpip2<@ 
m<pi<p2 m<pi<p2 m<pi<p2 
=) dD dL wB+OCe?) 
m<a3 m<piSr/Z P<P2s pire 
2 rat 
= S- S- |. (=) a _ m(p1)pr - 2 | , wt 
1 Pim) pym Pi 
m<x3s m<pis/ = 
+0(27) 
k 3 3 
Cyr av 
= “+0 ; 5 
re In’ a (<a -) (5) 
where c; (¢ = 1,2,--- ,&) are computable constants. 


Now we estimate the error terms in set C. Using the same method of proving (4), we have 


SL(n) = p? and Q(n) = 


S > Q(n)SL(n) 


nec 


Q(m) + 2p, so 


S- (Q(m) + 2p) p=2 S- p+ S- (Q(m)p”) 


mp? <x mp? <x mp? <x 
m<p m<p m<p 

3 

=a s s p> + O(a?) 


iL 
m<a3 M<PL/ = 


O(2”). 


(6) 


Finally, we estimate the error terms in set D. For any integer n € D, if SL(n) = p then 
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p< Jn; if SL(n) = p’, then p < n3; or SL(n) = p*,a > 3. So we have 


S72 Q(n)SL(n) «< S> Q(m)+p)p+ SY) (Q(m) + 2p)p 


ne€D mp<ax mp? <2 
pgm psm 
=o 2 
S> Q(m)ap)p* « ae (7) 
mpe <x 
ped a>3 


Combining (4), (5), (6) and (7) we may immediately obtain the asymptotic formula 


S>Q(n)SL(n) = S—A(n)SL(n) + SY Q(n)SL(n) 
n<u neA neB 

+ $2 O(n) SL(n) + S> O(n) SL(n) 

nec neD 
* dja? x 
= —+0(—__ }, 
2 In' x (<A -) 
where d; (¢ = 1,2,--- ,k) are computable constants. 


This completes the proof of Theorem. 
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Abstract For any positive integer n, the Pseudo-Smarandache-Squarefree function Zw(n) 
is defined as the smallest positive integer m such that m” is divisible by n. That is, 
Zw(n) =min{m: me N, n| m"}. In reference [2], Felice Russo proposed many problems 
and conjectures related to the Pseudo-Smarandache-Squarefree function Zw(n). The main 
purpose of this paper is using the elementary methods to study several problems in [2], and 
four of them are solved. 
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§1. Introduction and results 


For any positive integer n, the famous Pseudo-Smarandache-Squarefree function Zw(n) is 
defined as the smallest positive integer m such that m” is divisible by n. That is, 


Zw(n) =min{m: me N, n|m"}. 
It is easy to see that ifn > 1, then 


Zw(n) =| [p, 


pin 


where Il denotes the product over all different prime divisors of n (see reference [1]). From this 


pin 
formula, we can easily get the value of Zw(n). For example, Zw(1) = 1, Zw(2) = 2, Zw(3) = 


3, Zw(4) = 2, Zw(5) =5, Zw(6) = 6, Zw(7) =7, Zw(8) = 2, Zw(9) = 3, Zw(10) = 10, ---. 
In fact if m is a square-free number, then Zw(n) = n. In reference [2], Felice Russo studied the 
properties of Zw(n), and proposed the following four problems: 

Problem 1. Find all the values of n such that Zw(n) = Zw(n + 1) - Zw(n +4 2). 

Problem 2. Solve the equation Zw(n) - Zw(n+ 1) = Zw(n + 2). 

Problem 3. Solve the equation Zw(n) -Zw(n+ 1) = Zw(n+ 2) - Zw(n + 3). 

Problem 4. Find all the values of n such that S(n) = Zw(n), where S(n) is the Smaran- 
dache function. 


The main purpose of this paper is using the elementary methods to study these four 
problems, and solved them completely. That is, we shall prove the following conclusions: 
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Theorem 1. The following three equations have no positive integer solution. 


Zw(n) = Zw(n+ 1) - Zw(n + 2); (1) 
Zw(n) + Zw(n +1) = Zw(n + 2); (2) 
Zw(n) + Zw(n+ 1) = Zw(n + 2) - Zw(n + 38). (3) 


Theorem 2. There exist infinite positive integers n such that the equation S(n) = Zw/(n), 
where S(n) is the Smarandache function defined by S(n) = min{k: ke N, n| kl}. 


§2. Proof of the theorems 


In this section, we shall complete the proof of our theorems. First we prove that the 
equation Zw(n) = Zw(n+1)-Zw(n-+ 2) has no positive integer solution. It is clear that n = 1 
is not a solution of this equation. In fact if n = 1, then 


1 = Zw(1) 2-3 = Zw(2)- Zw(3). 
If n > 1, suppose that the equation (1) has one positive integer solution n = no, then 
Zw(no) = Zw(np + 1) + Zw(no + 2). 


For any prime divisor p of n = no, it is clear that p | Zw(no). 


From Zw(no) = Zw(no + 1) Zw(no + 2) 
we deduce that p| Zw(no +1) + Zw(ng + 2). 
That is, p| Zw(no +1) or p | Zw(no + 2). 


(a) If p | Zw(no + 1), then p | (no + 1), combining p | no and p | (mo + 1) we get 
p|no +1— no = 1, this is a contradiction. 
(b) Ifp| Zw(no + 2), then p | (no + 2), combining p | no and p | (no + 2) we deduce that 
p|no +2 — no = 2, then we get no = p = 2. 


From equation (1) we have 
2= Zw(2) #3-2 = Zw(3)- Zw(4), 


It is not possible. So the equation (1) has no positive integer solution. 

Using the similar method as in proving problem 1, we find that the equation (2) and 
equation (3) also have no positive integer solution. This completes the proof of Theorem 1. 

In order to prove Theorem 2, we need some important properties of the Smarandache 
function S(n), which we mentioned as the following: 


1 2 Ok; 


Lemma 1. Let n = pf{'p5’---py* is the prime powers factorization of n, then 


Proof. See reference [3]. 
Lemma 2. If p be a prime, then S(p*) < kp; If k < p, then S(p*) = kp, where k is any 


positive integer. 
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Proof. See reference [4]. 

Now we use these two simple lemmas to prove our Theorem 2. It’s clear that all prime 
p are the solutions of the equation S(n) = Zw(n). So there are infinite positive integers n 
satisfying the equation S(n) = Zw(n). 

Now we construct infinite composite numbers n satisfying the equation S(n) = Zw(n), 
let n = p - p++: Pe—1+ Py”, where p; are the different primes, and py > ag = pip2+ ++: Pk—1- 
This time, from the definition of S(n) and Zw(n) we have S(n) = pi: po-:+Pr—1* Pe and 
Zw(n) = pi: P2°+*Pr—1* Pk. So all composite numbers n = p; - p2---Ppe—1- py," (where p; are 
the different primes, and py > a, = pip2--+:pr—1 ) satisfying the equation S(n) = Zw(n). 

Note that k be any positive integer and there are infinite primes, so there are infinite 
composite numbers n satisfying the equation S(n) = Zw(n). 

This complete the proof of Theorem 2. 
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Abstract For any positive integer n, the famous pseudo Smarandache function Z(n) is 


defined as the smallest positive integer m such that n evenly divides 3 k. That is, Z(n) = 
k=1 
m(m + 1) 
2 
purpose of this paper is using the elementary and analytic methods to study the mean value 


min {m : n| , m€N >, where N denotes the set of all positive integers. The main 


properties of In Z(n), and give an interesting asymptotic formula for it. 
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§1. Introduction and results 


For any positive integer n, the famous pseudo Smarandache function Z(n) is defined as 
m 


the smallest positive integer m such that n evenly divides S- k. That is, 
k=1 
m(m +1) 


Z(n) = min fm : n| 5 


,meEN \ ; 
where N denotes the set of all positive integers. For example, the first few values of Z(n) are: 


Z(1)=1, Z(2)=3, Z(3) =2, Z(4) =7, 2(5) =4, Z(6) =3, Z(7) =6, Z(8) =15, 
Z(9) =8,Z(10) =4, Z(11) =10, Z(12) =8, Z(13) =12, Z(14) =7, Z(15) =5, 
Z(16) =31, Z(17) =16, Z(18) =8, Z(19) =18, Z(20) =15, ------ 


This function was introduced by David Gorski in reference [1], where he studied the elementary 
properties of Z(n), and obtained a series interesting results. Some of them are as follows: 

If p> 2 bea prime, then Z(p) = p—1; 

If n = 2", then Z(n) = 2**1 -1. 

Let p be an odd prime, then Z(2p) = p, if p = 3(mod4); Z(2p) = p — 1, if p = 1(mod4). 

For any odd prime p with p | n and n# p, Z(n) > p—1. 

The other contents related to the pseudo Smarandache function can also be found in 
references [2], [3] and [4]. In this paper, we consider the mean value properties of In Z(n). 
About this problem, it seems that none had studied it yet, at least we have not seen any related 
results before. The main purpose of this paper is using the elementary and analytic methods 


Vol. 3 On the pseudo Smarandache function 49 


study this problem, and give an interesting asymptotic formula for it. That is, we shall prove 
the following conclusion: 
Theorem. For any real number z > 1, we have the asymptotic formula 


S/n Z(n) =a#lnae+ O(a). 


Whether there exists an asymptotic formula for the mean value 


are two open problems. 


§2. Proof of the theorem 


In this section, we shall complete the proof of Theorem. First we need the following simple 
conclusion: 
Lemma. For all real number x > 1, we have the asymptotic formula 


> = =Inz+O(1), 


pSau 


where ‘ denotes the summation over all prime p with 2< p< za. 
pSau 
Proof. See Theorem 4.10 of reference [5]. 


Using this Lemma we can prove our Theorem easily. In fact for any positive integer n > 1, 


2n(2n — 1 
Ze). from the definition of Z(n) we know that Z(n) < 2n —1. So by the 


Euler’s summation formula we can get 


S "In Z(n) < $7 nn - 1) < eine + O(2). (1) 


n<a nsx 


note that n | 


Now let A denotes the set of all square-full numbers n (That is, if p | n, then p? | n) in the 


interval [1, 2]. Then we have 


S7inZ(n) = Sn Z(n) + $6 In Z(n). (2) 


neA ngA 


From reference [6] we know that 


> l= MDA + a +0 (x? exp (-—Clog? x(log log «)~*) ) ; 


n<ux 
neEA 


where C > 0 is a constant. By this estimate and note that In Z(n) < In(2n), we may get 


S/n Z(n) < Vz-Inz. (3) 


n<u 
neA 
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If n ¢ A, then n = 1 or there exists at least one prime p with p | n and p? + n. So from Lemma 


we have 
S > In Z(n) = S- In Z(np) > >> a In(p — 1) 
n<a npSa psx n<e 
ng€A (n, p)=1 (n, p)=1 
= S- F = sa + on) In(p — 1) 
= 2 
Pe 
1 1 
= S- Oe os Ge S- a + O(a) 
psx P psa - 
= «-n#+O(a). (4) 
Combining (1), (2), (3) and (4) we may immediately get the asymptotic formula 


S/n Z(n) =alna+ O(2). 


n<x 


This completes the proof of Theorem. 
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Abstract In this paper, I present some problems involving Mersenne primes for Scientia 
Magna. 


The primes 2? — 1 are called Mersenne primes. The set of p is 


T={y =2, =3, ¥3=5, 7, 13, 17, 19, 31, 61, 89, 107, 127,---}. 


If g = 2? — 1 is a Mersenne prime, then p is prime (by Cataldi and Fermat) and n = 2?~1g is 
an even perfect number (Euclides proved that every perfect number is of this type). There is 
an open problem about Mersenne primes: 

An open problem. Are there infinite many Mersenne primes? 

The numbers 2% — 1 (with a prime) are called Mersenne numbers. There are many open 
problems about Mersenne numbers: 

(1) Are there infinitely many composite Mersenne numbers? 

(2) Is every Mersenne number square-free? 


In this paper, I present the following sequence related with Mersenne primes 
eS Py Wop TOD? eg DS OO a OF ae tt OF OP as 


where Wy, = Wn, + 2771 for all n > 1. The prime factorizations of the first terms are: 
3, 7, 23, 3-29, 47-89, 13-31-1738, 3-7- 15803, 631 -1702177,---. 

My problems for Scientia Magna are: 

Are there infinitely many primes in the sequence w;? 

Are there infinitely many composite numbers in the sequence w;? 


Is every term of the sequence w; square-free? 
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denotes the least common multiple of 1, 2, ---, &. The main purpose of this paper is using 
the elementary methods to study the mean value properties of In SL(n), and give a sharper 


asymptotic formula for it. 
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§1. Introduction and Results 


For any positive integer n, the famous F.Smarandache LCM function SL(n) is defined as 
the smallest positive integer k such that n | [1, 2, ---, k], where [1, 2, --- , k] denotes the 
least common multiple of 1, 2, ---, k. For example, the first few values of SL(n) are SL(1) = 1, 
SL(2) = 2, SL(3) = 3, SL(4) = 4, SL(5) = 5, SL(6) = 3, SL(7) = 7, SL(8) = 8, SL(9) = 9, 
SL(10) = 5, SL(11) = 11, SL(12) = 4, SL(13) = 13, SL(14) = 7, SL(15) =5, ---. About the 
elementary properties of SZ(n), some authors had studied it, and obtained some interesting 
results, see reference [3] and [4]. For example, Murthy [4] showed that if n is a prime, then 
SL(n) = S(n), where S(n) denotes the Smarandache function, i.e., S(n) = min{m: n|m!, me 
N}. Simultaneously, Murthy [4] also proposed the following problem: 


SL(n) = S(n), S(n) An? (1) 


Le Machua [5] completely solved this problem, and proved the following conclusion: 
Every positive integer n satisfying (1) can be expressed as 


n=12 or n=pf'ps?--- perp, 


where pj, P2, °*:, Pr, p are distinct primes, and aj, a2, ---, a, are positive integers satisfying 
p> D;, v=, 2, 8 Rigo ds 
Lv Zhongtian [6] obtained the asymptotic formula: 


nT x? By x? x? 
SL(n) = —-— Ole |, 


In reference [7], Professor Zhang Wenpeng asked us to study the asymptotic properties 
of S- In SL(n). About this problem, it seems that none had studied it, at least we have not 


n<ux 
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seen related papers before. The main purpose of this paper is using the elementary methods 
to study this problem, and obtain a sharper asymptotic formula for it. That is, we shall prove 
the following: 


Theorem 1. For any real number zx > 1, we have the asymptotic formula 


S- InSL(n)=alnr+O(a). 
n<x 
Using the same method of proving Theorem 1 we can also give a similar asymptotic formula 
for the F.Smarandache function S(n). That is, we have the following: 
Theorem 2. For any real number x > 1, we have 


> InS(n) =a2nz2+O(z), 


n<ux 


where S(n) denotes the Smarandache function. 


§2. Proof of the theorems 


To complete the proof of the theorems, we need the following two important Lemmas. 

Lemma 1. For any positive integer n > 1, let n = p{'ps?---p%* denotes the factorization 
of n into prime powers, if aj > 2,a2 > 2,---a,; > 2, then we call such an integer n as a 
square-full number. Let A2(a) denotes the number of all square-full integers not exceeding z, 


then we have the asymptotic formula 


_) a 
C2) 


1 
ur? 4 


x3 4 O (x exp (—Clog? w(loglogx)~*)) , (2) 


where C > 0 is a constant. 

Proof. See reference [8]. 

Lemma 2. Let p be a prime, k be any positive integer. Then for any real number x > 1, 
we have the asymptotic formula 


S- Inp=anzr+O(za). 
pk<a 
(p, k)=1 
Proof. From the several different forms of the Prime Theorem (See reference [2], [7] and 
[9]), we know that 


yy 2? = ne +01), 


k<au 2 
Xx 
rs 
and 7 ; 
np , 
» Pp O(c) 


where D be an positive constant. 
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Applying these asymptotic formulas, we have 


S 7 Inp a 1 


M 
a 
3 
I 


pk<a psa k<e 
(p, k)=1 (p, k)=1 
= mp (Z-4+00) 
2 Pp 
PLT 
l ] 
= oy 2-2 +0(Lmp 
psu e psa P psa 
= «lmnx+O(z). 


This proves Lemma 2. 
Now, we use these Lemmas to complete the proof of our theorems. Let U(n) = S- In SL(n). 


n<ax 


First, we estimate the upper bound of U(n). In fact, from the definition of F.Smarandache LCM 
function SL(n), we know that for any positive integer n, SE(n) < n and In SL(n) < Inn, so we 


y InSL(n) < a Inn. 


n<au n<a 


have 


Then from the Euler’s summation formula (see reference [2]), we may immediately deduce 
that 


U(n) < So inn=a2lnz—2+O(Inz) =2lnz+O(z). (3) 


n<x 


Now we estimate the lower bound of U(n). For any positive integer n > 1, let n = 
pips? +: ps be the factorization of n into prime power, we divide the interval [1, n] into two 
subsets A and B. A denotes the set of all integers in the interval [1, n] such that a; > 2 
(i=1, 2, ---, s). That is to say, A denotes the set of all square-full numbers in the interval 


[1, n]; B denotes the set of all integers n with n € [1, n] but n € A. Then we have 


U(n) = S> mM SL(n) + S~ In SL(n). 


neEA neB 


From Lemma 1 and the definition of A, we have 


S¢ In SL(n) < S Inns Si nc=hns > 1=Ina: Ao(xz) «K Vring. (4) 


n<ux n<x n<u n<x 
neA neA neA neA 
Now we estimate the summation in set B. Since SL(n) = max{p{', p5?, --- , pe} (see 


reference [4]), so for any n € B, there must exist a prime p such that p|n and p? } n. Therefore, 
from the definition of SZ(n) we have SL(np) > p. Using this estimate we may immediately 
deduce that 


SS" InSL(n) = S- In SL(np) > S- Inp. (5) 


n<ux np<x np<x 
neB (n, p)=1 (n, p)=1 
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Then from Lemma 2 and (5), we have 


S° In$L(n) > amet O(a). (6) 
neB 


Combining (3) and (6), we may immediately deduce the asymptotic formula 


S/n SL(n) =alnrz+O(a). 


n<ux 


This completes the proof of Theorem 1. 


Using the same method of proving Theorem 1, we can also prove Theorem 2. 
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Abstract In this paper we give a generalization of the abc-theorem for functions of several 
variables. 
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81. Introduction 


Let F' be a fixed algebraically closed field of characteristic 0. Let f(z) be a non-constant 
polynomial with coefficients in F' and let r(f) be the number of distinct zeros of f. Then we 
have the following. 

abc-theorem.(([3]) Let a(z), b(z), c(z) be relatively prime polynomials in F' and not all 
constants such that a+ b= c. Then 


max {deg(a), deg(b), deg(c)} < r(abc) — 1. 


In [9], there is a generalization of the above theorem. The equality a+b = c is replaced by 
fo +-::>+ fn+1 = 0 and the functions f; are polynomials in several variables, relatively prime 
by pairs. In this paper, we give a similar result for the case of the functions f; have no common 
Zeros. 

Let f is a polynomial in several variables with coefficients in F' and f has a factorization 


s 
f=], 
t=1 


where the polynomials p; are irreducible, distinct, at most one of them is constant, and the 


a; > 0 are integers. Define 
No(f) = deg (IIe) : 
i=1 


The main theorem is the following: 
Theorem 1.1. Let fo, --- , fn4i be n+2 polynomials in several variables in F'la1,--- , x7] 
have no common zero such that fo,--- , f, are linearly independent. Assume also that 


FO ae Feet = Us (1) 
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Then 
n+1 

max deg fi; <n (s No(fi) - :) 
i=0 


O0<i<n+1 


The following corollary is useful for study on the problem of diophantine equations over 
function fields. 


Corollary 1.2. Let fo,--: , fn41 be n+2 polynomials in several variables in F'[a1,--- , xi], 
not all constants and have no common zero such that f9"°,--- , f7"" are linearly independent. 
Assume also that 

feet Rar =O; (2) 
where ™o,--- ,7%n+1 are positive integers. Then 
1 1 1 
feed Se 
mo Mn+1 n 


§2. Proof of the main theorem 


Let f is a rational function in several variables, we write f in the form: 


fi 
f =F 
fo 
Where fi, f2 are polynomial functions are non-zero and relatively prime on F'[a,--- , a]. The 


degree of f, denoted by deg f, is defined to be deg f; — deg fo. 
Let p is a non-constant irreducible polynomial, we write f in the form: 
fl 
f = oe 
g2 
such that p is not a divisor of gig2, where g, g2 are polynomials, then a is called the order of 
f at p and is denoted by j4. We have the following easily proved properties of p'. 
Lemma 2.1. Let f, g be two polynomials and p € F'[a1,--- , | is a non-constant irre- 
ducible polynomial, we have 
a) we, > min(u', py), 
b) Mfg = HF + Hh 
c) wh = pe — ph. 
g 
For A is a differential operator of the form 


of Ohm 
A= hie oo ea 
(p41 L ) On" Orr? 


where j1; > 0 are integers, we denote the rank of A by 


Note that in the above definition, whenever py; is zero, we admit the following convention: we 
oi : : 
7, is an identity operator. 


eliminate j; and consider the operator ay 
vy 
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Lemma 2.2. Let y be a polynomial in several variables satisfy the following Ay £ 0, p 
be a non-constant irreducible polynomial. Then 


Wig 2 PA) Fe. 
Proof. Let »2% = m, then there exists a polynomial f such that y = p™f. We have 


Of Op 
Xi 


Op _ 
da, | 5 


m—1 


). 


From this, we have 


Therefore, 
Moe 2 alsa. 
From this, we obtain 
Hay = —P(A) + He. 
Given Ao, ---, As, such that p(A;) < 7 and polynomials ho,--: ,hs in Flai,--- , xi], a 
generalized Wronskian has the form 


W{ho,--+ , hs] = det |Aihjlo<i,j<s- (3) 


A well-known result (see [7], [8]) asserts that, if the functions h; are linearly independent over 
F’, then there exists a generalized Wronskian of the (3), which does not vanish. 
Proof of Theorem 1.1. By the hypothesis, fo,--- , fn are linearly independent, we have 


there exists a generalized Wronskian W of fo,--- , fn does not vanish. We set 
P= W(fo,-°: dn) 
fo:++ fn 
for Fntt 
Q= =. 
W (fo,° = ota) 
Hence, we have 
fn4i = PQ. (4) 
We first prove that 
n+1 
degQ <n S~ No(fi)- 
i=0 


Suppose that p is a divisor of fofi--: fn41 and p is a non-constant irreducible polynomial. 
By the hypothesis, we have there exists v, 0 < vy <n+1, such that p is not a divisor of f,. By 
the hypothesis, fo+---+ fa + fn+1 = 0, we have 


Pp _ P 
Ho for-fn41 = fo:-fy-1fvdi- fn4ti 
W (fos sfn) Wforofu—as Fotis fn41) 
n+1 


Pp Pp 
= os Bg, ~ PW (fos fo-a Fotis fat1)” 
j=0 
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W(fo, ---, fv-1, fo4i, ++: fn41) is the sum of following terms 


bAo faoAt far a Anfans 


where a; € {0,---,n+1}\{v}, 6 = +1. By Lemma 2.1, 2.2 and note that p(A;) 


have 


p 
HA0 fag At fay Anfan 


ee et eee 


O<jsntl, p/ fj O<jsntl, p/ fj 


= doe, rm dl 
j=0 


O<jsnt+l, p/ fj 


P = 
lee rae x = 
0<j<n+1, D/f; 


IV 


By Lemma 2.1, we have 


n 
O<jsntl, p/ fj 


P 
Lt for fn41 <n S- 1. 


a O<j<n+, p/fj 


Pp P = 
WG ia Pacnoead) = Mie. hy, 


Hence 


By the definition of degree of a rational function, we have 
n+1 


degQ <n S_ No(fi). 


i=0 
Next, we will prove that 
deg P < —n. 
We have the determinant P is a summa of following terms 


jose Ar far -- An fe, 
Fofp.**+ Ibn 


For every term, we have 


deg ( Sole Setn Ants = deg (Ae) + deg (=) +++ dey (AZ) 
fe fa, fe, 


Feo fp. °** Fen 0 
< ‘spi Ay) —giAa) <4 py). 
Note that, n < p(Ao) +--+ p(An) < ee from this, we have 
deg P < —n. 


deg fn41 = deg P+degQ 


n+1 
n ( No(fi) — 7 
i=0 


IA 
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Similar arguments apply to the polynomials fp, fi, --:, fn, we have 


n+1 
max deg fi <n ( No(fi) — :) . 
i=0 


O0<i<n+1 


This completes the proof. 


Proof of Corollary 1.2. By Theorem 1.1, we have 


n+1 
mj;deg fj; =deg ff? < n(S> No(f;") -1) 
j=0 


n+1 


= n(S> No(fj) — 1) 
j=0 


n+1 


< n> deg f;. 
j=0 


Then, we have 


n+1 n+1 n+1 


So deg fi <n >. Yo dex fy 
j=0 j=0 J j=0 


From this, we obtain (2). This completes the proof. 
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Abstract This paper proposes an assessment method for weight of experts at interval recip- 
rocal judgment based on support vector domain description (SVDD). In this method, firstly, 
we give the point-vector decomposition to all the columns of every interval judgment ma- 
trix. Secondly, the group information is found using the support vector domain description 
(SVDD), and several concepts are introduced such as group compatibility, group core vectors 
and group information contribution ratio. Thirdly, we present a computational method for 
the contribution ratio of every expert so as to determine the assessment weight of every ex- 
pert. Finally, the feasibility and validity of the method are shown with one example. Using 
the method can decrease unilateral influence, extract key information, and evaluate weight of 
each expert more objectively. 

Keywords Group decision making, interval reciprocal judgment, compatibility, support 


vector domain description, weight. 


81. Introduction 


The key of ensuring the quality of group decision making is the determination of the 
weights of experts in the aggregation process. Therefore, the study on the weights of experts is 
a research hot point in group decision making [7]-[13]. 

Currently, the principal problem about the evaluation to weight of experts is without 
eliminating the disturbance information provided by experts in which the given weight of experts 
is doubting. We think that it is more scientific for the weight of experts that should beforehand 
remove the unilateral disturbance information, and then appraise the weight according to the 


magnitude of contribution to group opinion. 
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In addition, though there exist many evaluation methods for the weight of experts, we 
do not have yet find the evaluation method for the interval judgment. Thus, this paper will 
propose a more objective method for the weight problem of experts at the form of interval 
reciprocal judgment. 

The approach gives firstly a point-vector decomposition to each judgment matrix using the 
method in preference [4], and then extract group information based on the technique for SVDD 
[5]. Finally, determine the contribution ratio of group opinion of each expert. 


§2. Some concepts 


In this section, we briefly describe two basic definition related interval judgment. 

Definition 1.[3][4] Matrix A = (aj;)nxn is called a point judgment matrix, for all aj; 
(i, 7 =1,--- ,n), if the following conditions are satisfied: 

(1) 1/9 < ay < 9; 

(2) aig = 1/ayi; 
where the each column vector is called a point vector. We Simply call the above judgment way 
point judgment. 


Definition 2.[3][4] Matrix A = (a;;)nxn is called a interval reciprocal judgment matrix, 


for all aj; (¢,7 =1,--- ,n), if the following conditions are satisfied: 
(1) aij = [lig, wig], 4 I; 
(2) aaj = [1,1]; 


(3) lig = L/uje; 
where [l;;, wij] is an interval number, l,;, ui; represents the lower limit, upper limit of the 
significance ratio of element 7 to element 7, respectively. Here we take also the scale values 1-9. 
In following we simply call interval reciprocal judgment matrix interval judgment matrix. 
Suppose there exist n decision objects, and m experts: B®, B@), ..., E(™. Let the 
interval judgment matrix provided by E“*) (k = 1,--- ,m) be 


k k kK k 
i Tae 
k ik k k 
ye | ree) a En abn ) 
k k k k 
pa), ul) [ase ) ul) eo (1, 1] 


§3. Point-judgment decomposition 


In this section,we shall give a point-vector decomposition to interval reciprocal judgment 
matrix (1). All the elements of each column of interval judgment matrix (1) are interval numbers 
in addition to those elements that lie on main diagonal. 

If choosing randomly one of numbers among the lower limit and the upper limit of each 
interval number, which can form some n— dimensional point vectors together with the element 
1 at the diagonal. 
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Theorem 1. Interval judgment matrix (1) can generate 2”~!- n point vectors by decom- 
position. 

Proof. Since any column elements in interval judgment matrix (1) comprise n — 1 interval 
numbers. When taking two limit values at each interval number, we can obtain 2”~! kinds of 
ways taking numbers. Integrating the element 1 lie on diagonal, then each column can construct 
2”-! point vectors. 

We note that interval judgment matrix (1) comprises n columns, therefore, interval judg- 
ment matrix (1) can produce 2”~!-n point vectors(n dimension). For example, if n = 3, 
let 


(1, 1] lio, Ui2] [lis, urs] 
A=] [lo1, wai] (1, 1] [l23, u23] | > (2) 
(031, u31]  [l32, use (1, 1] 


then the first column of matrix A can be decomposed into 4 point vectors: 


oy = [1, lor, tsi)", a2 = [1, lar, wai)”, a3 = [1, war, tsi)”, a4 = [1, war, wail’ - 


So, matrix A can generate 23-1 . 3 = 12 point vectors altogether. 

For an n x n judgment matrix, if the 2”~'-n point vectors are normalized, they will lie on 
the same cone underside [4]. In fact, the convex combination of these point vectors can form a 
feasible domain of weight vectors of the decision objects [4]. For a group decision making, each 
expert can all generate such a feasible domain. 


§4. Support vector domain description (SVDD) [6] 


In this section, we shall introduce briefly SVDD algorithm. Let a set of data be x;, x; € 
R4(i = 1,--. ,N). The algorithm for SVDD is just the description of data by finding a hyper 
sphere with minimum radius containing nearly this data. 

In order to describe better the data, usually a nonlinear mapping ¢ will be mapped onto 
a high-dimensional feature space, and a kernel function k(x, y) = $(x) - é(y) will be introduced 
to instead of inner product. 

In order to eliminate the influence of outliers on the description, the slack variables €; and 
penalty factor C are introduced. Minimization of radius of the hyper sphere is just finding the 
solution of the following convex quadratic programming problem: 


N 
min R? + ey Gi, (3) 

i=l 
s.t. (b(a4) — a)" ($(xi) — a) < R? + &, (4) 
where €& > O(i = 1,---,N), and C is a constant that punishes the samples separated by 


mistake, R is the radius of the sphere, a (vector) is the center of the hyper sphere. To solve the 
optimization problem, we set lagrange function as follows: 


N N N 
L(R, a, 03, Yi, &) = Soa [l(a ||? — 24 b(xs) + |lal|? — R? —&)| +R? 4 Cy E-S vi (5) 
i=l i=l i=l 
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where a; > 0,7; > 0. 
By KKT condition, the dual programming of above convex quadratic programming [6] is 


N N oN 
max ) > ajk(a;, v5) —SO YM aia; k(x, @5), (6) 
i=1 i=1 i=1 
N 
at.) a= 1, Vs ae OO (7) 
i=1 


We will often have that a; = 0 for most 2x;, all the x; for which a; 4 0 are called the 
support vectors. And the center of the sphere can be expressed as 


N 
a=) aga), 0< a4 <C. (8) 


i=l 


For some x; for which a; 4 0, sphere radius is 


N ON 
R? = ||6(x;) — all? = k(a;,2;) 23 ok Diy J+ S03 o auc k(xj, 2), (9) 


\|z —al|? = (z-z) )=2 ot (xj, 2 ps ajyajk(aj,2;) < R?. (10) 


Note. All the test vectors will satisfy equation (10) if C > 1. 


§5. Determination method for contribution ratio of group 


information 


From the discussion in section 2 we know that if judgment matrices are given at the form 
(1), the point-vector decomposition can all be performed. 

By Theorem 1, we can construct 2”~!-n point vectors from each judgment matrix by the 
decomposition. Each point vector represents judgment information corresponding expert. 

We hope that the weight of each expert can be evaluated in terms of the magnitude of the 
his (or her) contribution ratio to group opinion. Generally, all the point vectors of each expert 
comprise either personality or group property. 

Our concern is how to extract group information. We think that group opinion is just the 
nearer opinion among all the experts, that is the opinion indicated by the point vectors with 
smaller distance among all the point vectors of all the experts. 

By reference [4], the convex domain formed by the convex combination of all the point 
vectors of each expert can be seen as the feasible domain of the weight vector of order alterna- 
tives. It is more reasonable that the overlapping domain of feasible domains of all the experts 
is viewed as the representative domain of group information. However, to find precisely the 
overlapping domain is rather difficult. Therefore, it is needful to find an approximate method 


for convenience in computation. 
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This paper attempts to train a minimum enclosing sphere respectively by the 2"~!-n point 
vectors generated by each expert using the method for SVDD. And then we shall determine 
the evaluation weight of each expert based on the information provided by the point vectors in 
public overlapping domain of all the enclosing sphere. 

5.1. Determination of penalty parameter [5] 

The key in using SVDD to train sample is the choice of parameter. For each judgment 
matrix, 2”~!-n point vectors contain likely outliers (or noises). However, these noises are likely 
group information with significant value for whole group. Therefore, we shall take parameter 
C’ = 1 when using SVDD to describe the point vectors of each expert so as to all the point 
vectors can be contained in the trained sphere. 

5.2. Choice of kernel function [6] 

Choice of kernel function must be done in using SVDD. We only consider Gauss kernel 


function as follows 
k(x, y) = exp[-||x — y||?/s?]. 


It is well known that the enclosing sphere generated by Gauss kernel function is tighter if 
width parameter s is set a smaller value, and the domain enclosed by the sphere is smaller than 
the feasible domain formed by the convex combination of point vectors. 

By reference [5], we know the enclosing sphere generated can generate a sphere surface 
as the distribution of the point vectors when s = 1. Such an enclosing sphere is too tight 
though it can reduce redundant space in the sphere.However, when s = 5, the enclosing sphere 
generated is a ellipsoid, which is larger than the domain formed by the convex combination of 
point vectors. 

Therefore, in order to more approximate to the ideal feasible domain, we let s = 2 in this 
paper. 

5.3. Concepts and computation method 

In this section, we shall propose some concepts related to training spheres and the compu- 
tation method for group information contribution ratio. Let the least enclosing sphere be 5S; 
generated by 2”—!-n point vectors from the judgment matrix A) provided by expert E“), and 
the sphere center be az, radius be Ry(k =1,--- ,m). Suppose S;()5; represents the common 
overlapping domain of S;and S;. 

Definition 3. If there exist at least a point vector V; € S;, and a point vector V; € Sj, s.t. 
VY, € S; [| 8; and V; € S;[]S,, then call 5;(]|S; 4 O(7,9 = 1,°** ,m, 14 9). 

Definition 4. If the 2”~!-n point vectors generated by each judgment matrix are com- 
pletely same, then call the group completely compatible, or we think the group is of complete 
compatibility. 

Note. The case of complete compatibility is impossible to appear in general, usually, part 
compatibility exists. 

Definition 5. If $;()S; 4 ¢(i,j = 1,--- ,m,i 4 j),then we think that experts E and 
E®) exist compatibility. 

Definition 6. If S = 5:()S2(\--- Sm # ¢, then we say that the expert group exists 
group compatibility (m is the number of experts). 
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Definition 7. If S = Si ()\S2f)-:-QSm 4 ¢, then call point vectors in S = $1] S2f}--- 
()Sm group core vectors. 

Definition 8. If S = $;()S2()---()Sm 4 ¢, suppose S$ contains M point vectors, where 
rz point vectors belong to S,(k = 1,---,m), then call w, = r,/M group information contri- 
bution ratio of BE“), or weight. The obtained vector w = (w1,wW2,--- ,Wm)7 is called weight 
vector of the group experts. 

Note. Group opinion is determined completely by the core vectors, and those points outside 
the minimum sphere are viewed as noise points for the group, therefore, the method has the 


efficacy de-noising. 


§6. Algorithm 


The steps of evaluation method are follows: 

(1) Perform a point-vector decomposition to each judgment matrix, find all the point 
vectors of each judgment matrix; 

(2) Carry out a training using SVDD, find the least enclosing sphere of each expert; 

(3) Judge whether the group is of compatibility; 

(4) If the group compatibility exists, then calculate the group information contribution 
ratio of each expert, i.e. weight. 


§7. Example 


Suppose interval reciprocal judgment matrices provided by 3 experts for 4 evaluation al- 
ternatives are as follows: 


[1,1] [2,3] 4,6 7,8 
ao —|} (eal [et 23) B4 | 
aia [Fl 1,1 2,3 
a7 lava) tegl Tt 
[1,1] [2,4] [4,5] [6,7 
AQ = +3] (11) [2,3] [2,3 | 
pal [3a] [1] 2.3 
+4) (3.3) [3.3] [1 
[1, 1] [3,4] 3,5 6,8 
5 eae I 
4@—| wal 1 23) B38 
ell [ae 1,1 3,4 
sicl dagal Laval el 


By the steps in section 5, we can get 32 point vectors from each column of each judgment 


matrix by the decomposition, and 96 point vectors can generated from 3 judgment matrix 
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altogether. Here Gauss kernel (s = 2) is chosen. After the training of SVDD, we find S ¥ 4, 
which means that the group exists compatibility. The training result tells us that S contains 
64 group core vectors, i.e. M = 64, and r; = 30, rg = 18, r3 = 16. Then w, = 30/64, wo = 
18/64, w3 = 16/64. 

Therefore, the obtained weight vector is 


w = (0.469, 0.281, 0.250)”. 


Remark. 

(1) The core vectors in S can gradually decrease as the increase of the number of experts, 
and non compatibility can likely occurs, which accords with fact. Because the more expert, the 
more difficult keeping the compatibility. In this case, we can set 0 to the weight value of the 
experts with weaker compatibility. 

(2) When the number of judgment object is larger, the number of point vectors is more, 
In this case, we can draw randomly a sample in each point vector set, and then use the method 


again. 


88. Conclusion 


The paper is mainly applying SVDD to adverse judgment problem in group decision mak- 
ing, and gave a more objective evaluation method for the evaluation quality of expert at interval 
reciprocal judgment. The contribution in this paper is mainly two aspects as follows: 

(1) Gave an evaluation method for expert weight aiming at interval reciprocal judgment 
by decomposing the interval judgment into point judgment. 

(2) The technology for SVDD was successfully applied to the group decision making at 
interval reciprocal judgment, which can eliminate disturbance information and find group in- 
formation, then the evaluation quality of each expert is more fairly judged. 

Farther work should consider the application of SVDD to the other group decision making. 
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Abstract For any positive integer n, the famous F.Smarandache LCM function SL(n) is 
defined as the smallest positive integer k such that n | [1,2,--- ,&], where [1,2,--- ,&] denotes 
the least common multiple of 1,2,--- ,4. The main purpose of this paper is using the elemen- 
tary methods to study the mean value distribution property of (P(n) — p(n))SL(n), and give 


an interesting asymptotic formula for it. 


Keywords SL(n) function, mean value distribution, asymptotic formula. 


§1. Introduction and Result 


For any positive integer n, the famous F.Smarandache LCM function S'L(n) defined as the 
smallest positive integer k such that n | [1, 2, --- , &], where [1, 2, --- , k] denotes the least 
common multiple of 1, 2, ---, &k. For example, the first few values of SL(n) are SL(1) = 1, 
SL(2) = 2, SL(3) = 3, SL(4) = 4, SL(5) = 5, SL(6) = 3, SL(7) = 7, SL(8) = 8, SL(9) = 9, 
SL(10) = 5, SL(11) = 11, SL(12) = 4, SL(13) = 13, SL(14) = 7, SL(15) = 5, ---. From the 
definition of SZ(n) we can easily deduce that if n = p? pS? ---p®" be the factorization of n into 


primes powers, then 
SL(n) = max{p{"," os ae | eee (1) 


About the elementary properties of SZ(n), many people had studied it, and obtained some 
interesting results, see references [1], [2] and [3]. For example, Murthy [1] porved that if n be 
a prime, then SL(n) = S(n), where S(n) be the F.Smarandache function. That is, S(n) = 
min{m: n|m!, me N}. Simultaneously, Murthy [1] also proposed the following problem: 


SL(n) = S(n), S(n) An? (2) 


Le Machua [2] solved this problem completely, and proved the following conclusion: 


Every positive integer n satisfying (1) can be expressed as 


1 2 


n=12 or n=pi'ps?-:- pep, 


where pj, p2, °*-, Pr, p are distinct primes and aj, a2, ---, @, are positive integers satisfying 


p>pe,i=1,2,---,r. 
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Zhongtian Lv [3] studied the mean value properties of SZ(n), and proved that for any fixed 


positive integer k and any real number zx > 1, we have the asymptotic formula 


nT x? By 2? x 
She t | o(=ar), 


a 
n<u i=2 Ine 


where c; (¢ = 2,3,--- ,&) are computable constants. 
Jianbin Chen [4] studied the value distribution properties of SZ(n), and proved that for 
any real number x > 1, we have the asymptotic formula 


Tsu) —Pemy=2-¢(2)- +0 ( | 


n<x 


where ¢(s) is the Riemann zeta-function, and P(n) denotes the largest prime divisor of n. 
Xiaoyan Li [5] studied the mean value properties of P(n)SL(n) and p(n)SL(n), and give 

two sharper asymptotic formulas for them, where p(n) denotes the smallest prime divisor of n. 
Yanrong Xue [6] defined another new function SZ*(n) as follows: SL*(1) = 1, and if 


n = pt'ps?---pe be the factorization of n into primes powers, then 


SL* (n) = min{py", p>’, wee 3 ob (3) 


where pi < po < +++ <p, are primes. 

It is clear that function SL*(n) is the dual function of SL(n). So it has close relationship 
with SL(n). About its elementary property of the function SZ*(n), Yanrong Xue [6] proved 
the following conclusion: 


For any positive integer n, there is no any positive integer n > 1 such that 
s 
+ 
an SL*(d) 


is an positive integer, where S- denotes the summation over all positive divisors of n. 
d\n 
In this paper, we shall study the value distribution properties of (P(n) — p(n))SL(n), and 


give a sharper asymptotic formula for it. That is, we shall prove the following: 


Theorem. For any real number x > 1 and any positive integer k, we have the asymptotic 


formula 
Ed as 
P(n) — p(n))SL(n) = ¢(3) - 28 - 2 +0(—ar), 
LPO) ~ Plo) SH) = 68) 2° +0 (Sears 
1 
where ¢(s) is the Riemann zeta-function, 6; = 3? b; (i = 2,3,--+ ,k) are computable constants. 


§2. Proof of the theorem 


In this section, we shall complete the proof of the theorem directly. For any positive integer 
n > 1, we consider the following cases: 


A:n=ny-p,m <p,and SL(n)=p; 
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B: n=ng-p,n2>p,and SL(n) =p; 
C:n=m- p*, a> 2, and SL(n) = p*; 
Now, for any positive integer n > 1, we consider the summation: 
So (P(e) = p(n))SL(n). 
n<u 
It is clear that if n € A, then from (1) we know that SL(n) = p. Therefore, by the Abel’s 
summation formula (See Theorem 4.2 of [7]) and the Prime Theorem (See Theorem 3.2 of [8]): 


—_ In 
where a; (t= 1,2,...,&) are computable constants and a; = 1. 
We have 
Yo (P(r) = p(n) SL(n) = So (P(r) = v(n)) SL(n) 
n<a nx 
neEA n=n1-p,ni<p 
SL(n)=p 
= SS SS Plru-p)- lr -p))p 
Ni SV MSPS Fy 
= SS D>S &=pP(r))p 
1 SVE NSPS yy 
= © dX P- dX} Dd orlmdp, (4) 
ny SVJe M1 SPS ny SVe MSPS FP 
while 
5 x? x my 3 
p= —a(—]— f  2yn(y)dy + O (nj) 
a ny ny ny 
ni <Ve MSPS Fy nisVJr 
3 Kk 3 
x b; By 
=~ Sate +o( es) | 
nicVJe 1 j=l ni 1 
k 
b; x? ) 
= 3)-2°- —+0O ; 5 
G8) Dieta (sa (5) 
: : ; 1 
where ¢(s) is the Riemann zeta-function, b; = 3? b; (i = 2,3,--+ ,k) are computable constants. 


Note that p(ni) <1, we have 


SSS plump = SYS pm) SS pv 


n1 SVE MSPS nyi<VJx ni SPSay 


ll 
= 
3 
~~~ 9 

| 8 
y 
a™~ 
Fla 
SS 
| 
s— 
a 
= 
Q 
ed 
+ 
1) 
oe 
So 
— 


«< SY pm) ——«K YY — =002?). (6) 
nixve 
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From (4), (5)and (6) we have 


k 
(Pen) = ple) Sin) = 6(3) 2° +0 (=), (7) 


n<u i=l 
neA 


where b; = =, }; (¢ = 2,3,--- ,&) are computable constants. 
Ifn € B, SL(n) = p, then by the Abel’s summation formula and the Prime Theorem, we 
can deduce the following: 


do (P(n) — p(n))SL(n) = do (P(@) = p(n) SL(n) 


na n<a 
neB N=N2-p,Nya>p 
SL(n)=p 


= SO (p—p(ne))p 


N2-PSx 
n2>p 


Dap Xa” 


N2-pSx p<Ja P<ng< 
n2>p 


If n € C, then SL(n) = p*, a > 2. Therefore, using the Abel’s summation formula and 
the Prime Theorem, we can obtain: 


Yo (P(r) = (n))SL(n) 


I 
= 
S 

| 
. 
So 
uw 
cs 
S 


n<a nx 
nec n=mM-p* ,a>2 
SL(n)=p% 
= So (Pim: p*) = p(m- p*))p® 
mp <x 
a>2 
« Pome yd De 
mp <x pO Sa MS fe 
a>2 a>2 
x 
<> = ee 
po Sa po <a 
a>2 a>2 


I 
8 
aM 
Ss, 

I 
8 
ae 
Ss, 
nN 
Be 
S 
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Now, combining (7), (8) and (9) we may immediately obtain the fowllowing asymptotic 
formula: 


k 
bj a? 
P(n) — p(n))SL(n) = €(3) - 2? - W+0(—ar), 
2PM) SHO =O 8") ag +O pea, 
where P(n) and p(n) denote the largest and smallest prime divisor of n respectively, ¢(s) is the 
Riemann zeta-function, b; = —, b; (¢ = 2,3,--- ,&) are computable constants. 

This completes the proof of Theorem. 
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Abstract For any positive integer n, the famous Pseudo Smarandache function Z(n) is 


defined as the smallest integer m such that n evenly divides Sok That is, Z(n) = 
k=1 


i 
ey) , m€N >, where N denotes the set of all positive integers. The main 


: m 
ming m: n| 
purpose of this paper is using the elementary method to study the properties of the Pseudo 
Smarandache function Z(n), and solve two conjectures posed by Kenichiro Kashihara in ref- 
erence [2]. 


Keywords Pseudo Smarandache function, conjecture, unbounded. 


§1. Introduction and Results 


For any positive integer n, the famous Pseudo Smarandache function Z(n) is defined as 
m 


the smallest positive integer m such that n evenly divides 3 k. That is, 
k=1 
m(m + 1) 


Z(n) = min fm | 


,meEN \ 3 
where N denotes the set of all positive integers. For example, the first few values of Z(n) are: 


Z(1)=1, Z(2)=3, Z(3) =2, Z(4) =7, 2(5) =4, Z(6) =3, Z(7) =6, Z(8) =15, 
Z(9) =8,Z(10) =4, Z(11) =10, Z(12) =8, Z(13) =12, Z(14) =7, Z(15) =5, 
Z(16) =31, Z(17) =16, Z(18) =8, Z(19) =18, Z(20) =15, ------ ; 


This function was introduced by David Gorski in reference [3], where he studied the ele- 
mentary properties of Z(n), and obtained a series interesting results. For example, he proved 
that if p > 2 is a prime, then Z(p) = p—1; If n = 2", then Z(n) = 2**+! — 1. The other 
contents related to the Pseudo Smarandache function can also be found in references [2], [4] 
and [5]. Especially in reference [2], Kenichiro Kashihara posed two problems as follows: Are 
the following values bounded or unbounded? 

A) [2o+2)— oh 

n+ 
B) aCe 
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About these two problems, it seems that none had studied it, at least we have not seen 
related papers before. In this paper, we use the elementary method to study these two problems, 
and prove that they are unbounded. That is, we shall prove the following conclusion: 

Theorem. For any positive number M large enough, there are infinitely positive integers 
n, such that 

Z(n+1) 
Z(n) 


From this theorem, we know that |Z(n + 1) — Z(n)| and 


> M and |Z(n+1)—- Z(n)| > M. 


Z(n+ 1) 
Z(n) 


solved two problems posed by Kenichiro Kashihara in reference [2]. 


are unbounded. This 


§2. Proof of the theorem 


In order to complete the proof of the theorem, we need the following important conclusion: 

Lemma. Let & and h are any positive integers with (h, k) = 1, then there are infinitely 
many primes in the arithmetic progression nk + h, where n= 0, 1, 2, 3, ---. 

Proof. This is the famous Dirichlet’s Theorem, see reference [6]. 

Now we use this Lemma to complete the proof of our Theorem. In fact for any positive 
number M, we take positive integer m such that 2" > M. Note that (Cai 2h 1) = 1, 
so from Dirichlet’s Theorem we can easily deduce that there are infinitely many primes in the 
arithmetic progression: 


g?mtik4+o"41, where k=0, 1, 2,---. 


Therefore, there must exist a positive integer kg such that 2?"+1kyg + 2™+1=P bea 
prime. For this prime P, from the definition and properties of Z(n) we can deduce that 
TP) =P = lao hg +2", 
Z(P -1) = Z(2?™ 1 ko + 2™) = Z(2™(2"™* ky +:1)). 


Since 
“ss _ amtl eg (2™+1 ey +1) 
w=1 2 
2a™+1 ko 


and 2™(2"*1ky + 2™) evenly divides S- i, so we have 
i=l 


Z(P —1)<2™"'k. 


Thus 
Z(P) i g2mtlky +2” 


> 2” > M. 
Z(P—1) gmtIky 


P 
—— is ded. 
Z(P ~ 1) is unbounde 
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Similarly, we have 


|Z(P) — 2(P —1)| 


IV 


|Z(P)| — |2(P — 1)| 
OF he al gm = yl 
gmtte(2” —1) +2" > 2" > M. 


IV 


So |Z(P) — Z(P — 1)| is also unbounded. 
Since there are infinitely positive integers m, such that 2” > M, so there are infinitely 


Z 
positive integers n, such that |Z(n + 1) — Z(n)| and 


(n+ 1) 


~ Ba). are unbounded. 


This completes the proof of the theorem. 
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Abstract In his book “Comments and Topics on Smarandache Notions and Problems”, 
Kenichiro Kashihara asked us to find all positive integer solutions of the Diophantine equations 

a’. b*) oP cor (a’ Df. *) atbte _ d. He also said, “The question I would like to ask is 
whether there exists positive integer solutions where a ¥ b for either of the equations. If such 
non-trivial solutions exist, find a general method for generating them”. In this paper, we using 
the elementary method to study these two problems, and prove that the first equation has no 
non-trivial solutions. The second equation has infinite non-trivial solutions. This solved the 


problems proposed by Kenichiro Kashihara in his book. 


Keywords Diophantine equation, positive integer solutions, elementary method. 


§1. Introduction and results 


It is a very important content in number theory that to find all integer (or positive integer) 
solutions for some Diophantine equations. For example, the famous Fermat’s last theorem is 
that the Diophantine equation 

ze” +y” = 2” 
has no positive integer solutions for all integer n > 3. 

In his book “Comments and Topics on Smarandache Notions and Problems”, Kenichiro 

Kashihara asked us to find all positive integer solutions of the Diophantine equations 


(a?) =e (1) 
or 
(a? -b°-c*) TEE d. (2) 


He also said: “The question I would like to ask is whether there exists positive integer 
solutions where a ¥ 0 for either of the equations. If such non-trivial solutions exist, find a 
general method for generating them”. In this paper, we using the elementary method to study 
these two problems, and prove that the Diophantine equation (1) has no non-trivial positive 
integer solutions, and the Diophantine equation (2) has infinite non-trivial positive integer 


solutions. That is, we shall prove the following conclusions: 
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Theorem 1. All positive integers a, b and c satisfying the equation 


(a° : b*) ats =¢ 


if and only if a=b=s and c= s?, where s is any positive integer. 

Theorem 2. If (a, b, c) =t, a= ait, b= bit, c= cyt and (a1, 61) = (a1, C1) = (C1, 61) = 
1, then a, b, c and d does not satisfy the equation (2); 

If (a,b, d= (8, th. Py, & ae, 2), AL, 4, 4%. 26) or @, 46, ¢, 20°), then a, b,c 
and d satisfying the equation (2), where t is any positive integer. 


It is clear that our Theorem 1 and Theorem 2 solved two problems proposed by Kenichiro 
Kashihara in reference [2]. Whether there exists any other positive integer solutions for the 
equation (2) is an open problem. 


§2. Proof of the theorem 


In this section, we shall prove our Theorems directly. First we prove Theorem 1. For any 
positive integer a, b and c, if they satisfy the Diophantine equation (1), then we have 


(a? - 6%) FF =o, (3) 


Let (a, b) = d denotes the Greatest Common Divisor of a and b, and a = d-a,, b= d-b;. Then 
(a1, 61) =1, and from equation (3) we have 


ae pa 
ays 1 oq 1 ae 
ay by -d=c. (4) 


2b4 2a, 
In equation (4), since (a1, b:) = 1, c and d are positive integers, so aj**’! and b/1*"! both 
2 


must be positive integers. If a, = b,, then (4) become a, - b; -d? = a? = c. So all positive 


integers a = b and c = a? must be the solutions of the Diophantine equation (1). If a1 4 bi, 
2b 


1 
then without loss of generality we assume that a; > b;. Now we prove that a{!*"! can not be 
2by 
a positive integer. In fact for any prime divisor p of a with p® | a, and p®t! } ay, if af!" 
2ab 2ab 
be an integer, then parry also be an integer, and so 3 must be an integer. Note that 
ay 1 
2ab 
(a1, b1) =1, (a1 + b1, by) = 1 and 7 be an integer, we can deduce that 
ay 1 
ay+ by | 2a. (5) 


It is clear that p® | a1, so the formula (5) implies that 
2a > a,4+ 6, > a, = p*- aa, (6) 


where a2 be a positive integer. 

Since p be a prime and a be a positive integer, we must have p® - ag > p®% > 2% > 2a. 
So the inequality (6) is impossible. Therefore, all positive integers a, b and c satisfying the 
Diophantine equation (1) if and only if a = b and c = a?. This proves Theorem 1. 


Vol. 3 Diophantine equations and their positive integer solutions 79 


Now we prove Theorem 2. If positive integers a, b, c and d satisfying the Diophantine 
equation (2). Let (a, 6, c) =t, a=t-a1,b=t-b; and c=t-c, then (a1, bi, c1) = 1, so from 
the equation (2) we have the identity 


3by Bey Bayz 
b b b 
an 1+Fe1 > a 1F¢1 eet 1¥e1 Pp = d. (7) 


If ay = 6) = cy in the above formula, then note that (a1, b1, c1) = 1, we have ay = b) =c, = 1 
and d = t®. So all positive integers a = b = c and d = a are the solutions of the Diophantine 
equation (2). 

If a, = b1 < cy, then we can change the equation (7) into 


3(ay+e1) 3a} . 
2ay+ey Zay+ey = 
ay “Cy tf =d. (8) 
3(ay-ter) ‘ 3a, 
Note that (a1, c1) = 1, so a, “!~* and c,*!**! must be integers. Therefore, for any prime 


3a,Q@ 


2a + Cy 
Since (2a, +1, c1) = 1, we deduce that (2a, + c1) | 3a and 3a = k(2a; + c1). So 


divisor p of c, with p® | c, and p®t! + 1, must be an integer, or (2a; + c1) | 3a1a. 


3a > 2a, + cy = 2a, + pcg > 2a, + p* 


or 
3a — 2a, > p®. 


This inequality holds if and only if ay = 1, p = 2 =a. This time, a = b = t, c= 4t and d = 2t°. 
So for any positive integer t, a = b = t, c = 4t and d = 2¢° are the solution of the equation (2). 


3(a1+c¢1) 
If ay = b; > c, then for any prime divisor p of a, with p® | a, and p++ + ay, a,“ 


must be an integer, or (2a; + c;) | 3(a, + ¢1)G. Since (2a, + c1, a, +1) = 1, we deduce that 
(2a; +1) | 36 or 36 = k(2a, + ¢1). So 3G > 2a, +1, 38—1 > 2p%. This is not possible, since 
33 —1 < 2p® for all positive integer (. 

If a; > b1 > c1, first we prove that if (a1, b1) = (a1, c1) = (c1, 61) = 1, then (7) is not 
possible. In fact if (a1, b1) = (a1, c1) = (c1, 61) = 1, then 


Bby Bey Bay 


+01 + +01 + ay Foy +e, 
a;! 1 “=, byt 1 es and cy? 1Fe1 


all must be positive integers. Let (a1, 61 +1) = u, (b1, a1 +41) = 0, (C1, a2 +1) = w. 
Note that (a1, bi, ci) = 1, so it is clear that (u, v) = (v, w) = (w, u) = 1. If gq = 1 or 
3 


eal 
2, then b/!~'*** is not a positive integer. So without loss of generality, we can assume that 
3a 


1 
c, > 3. Ifu < a, then c{1™*" is not an integer. Otherwise, for any prime divisor p 


3 
with p® | c¢, and p*t! + aq, 2" must be an integer, so 3a,a = t(a, +6, +c). or 
aj+b, +c 
3a,Q a, +b, +c a, a,+b,4+¢ 
ees hs coi ) note that ( A, rate) = 1, so that “4 | ¢t. Therefore, we 
u u u u 
have 


Ba = ta (a1 +61 + e1) > ar tbr +e1 > 3e1 +3 > 3p" +3. 


Bbq 3cy 


This is not possible. Similarly, if v < /b; or w < \/q@, then aft*""**! or bf! 2*™ is not an 
integer. So we can assume that u > ,/aj, v > Vb; and w > ,/cy. 
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If Jay <w< 5, then note that w| ar +b +e, 0 | ar +b ter and w | a, + 6; + c1, so 


3a 


cj1*"* is not an integer. In fact this time, for any prime divisor p with p® | cy and p°t! f c1, 
3a,a@ . 

—_—— must be an integer, so 3a,a = t(a; + b1 + c1) or 

ajyt+b+c¢ 


a 
3a =to-— -u- wa > 2v-w > 2c, > 2p. 
u 
This is not possible. 
Bc, 

Similarly, if /c1 < w < s then bj1*"!** is not an integer. 

If w > a and u > al then note that u | a, and w | c,, so this time, we have a, = u and 
cy = w. From u | ai+6, +c, w | ar +b: +c, and (u, w) = 1 we have u-w | ai. +0, +c. 
Therefore, a, +6; +c, > u-w=ayz-cy > 3a, > a, +6; +c. This is not possible. 

Therefore, if a, > b; > cy, and (a1, b1) = (a1, 41) = (41, 61) = 1, then (7) is not possible. 


From the above we know that if (a, b, c) = t, a = ait, b = bit, c = cyt and (a1, bi) = 
(a1, C1) = (c1, by) = 1, then a, b, c and d does not satisfy the equation (2). If (a, b, c, d) = 
(t, t, t, #9), (t, t, 4t, 243), (4t, t, t, 2t3) or (¢, 4¢, t, 2¢8), then a, b, c and d satisfying the 
equation (2). 

This completes the proof of Theorem 2. 

Whether there exists any other positive integer solutions for the equation (2) is an open 
problem. 
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Abstract The main purpose of this paper is using the elementary method to study the prop- 
erties of the Pseudo-Smarandache function Z(n), and proved the following two conclusions: 
The equation Z(n) = Z(n +1) has no positive integer solutions; For any given positive integer 
M, there exists an integer s such that the absolute value of Z(s) — Z(s + 1) is greater than 
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§1. Introduction and results 


For any positive integer n, the Pseudo-Smarandache function Z(n) is defined as the smallest 
positive integer m such that [1+ 2+3-+---+m] is divisible by n. That is, 


ae 


2(n) =min fms men: n| 5 


where N denotes the set of all positive integers. For example, the first few values of Z(n) are: 
Z(1) =1, Z(2) =3, Z(3)=2, Z(4)=7, Z(5) =4, Z(6) =3, Z(7)=6, Z(8) =15, Z(9) = 
8,Z(10) = 4, Z(11) =10, Z(12)=8, Z(13)=12, Z(14) =7, Z(15) =5,---- 

In reference [1], Kenichiro Kashihara had studied the elementary properties of Z(n), and 
proved some interesting conclusions. Some of them as follows: 

For any prime p > 3, Z(p) = p—1; 

For any prime p > 3 and any k € N, Z(p*) = p* — 1; 

For any k € N, Z(2*") = 2*+1 —1; 

If n is not the form 2” for some integer k > 0, then Z(n) <n. 

On the other hand, Kenichiro Kashihara proposed some problems related to the Pseudo- 
Smarandache function Z(n), two of them as following: 

(A) Show that the equation Z(n) = Z(n + 1) has no solutions. 

(B) Show that for any given positive number r, there exists an integer s such that the 
absolute value of Z(s) — Z(s + 1) is greater than r. 

For these two problems, Kenichiro Kashihara commented that I am not able to solve them, 
but I guess they are true. I checked it for 1 < n < 60. 


82 Su Gou and Jianghua Li No. 4 


In this paper, we using the elementary method to study these two problems, and solved 
them completely. That is, we shall prove the following: 

Theorem 1. The equation Z(n) = Z(n + 1) has no positive integer solutions. 

Theorem 2. For any given positive integer MM, there exists a positive integer s such that 


|Z(s) - Z(s+1)| > M. 


§2. Proof of the theorems 


In this section, we shall prove our theorems directly. First we prove Theorem 1. If there 
exists some positive integer n such that the equation Z(n) = Z(n+1). Let Z(n) = Z(n+1) =m, 
then from the definition of Z(n) we can deduce that 


1 1 
n j eee) ) n+1 pa ) 
2 2 
Since (n, n+ 1) = 1, we also have 
1 1 1 
nt (pg EO ee) 
2 2 2 
Therefore, 
nom. (1) 


On the other hand, since one of n and n+ 1 is an odd number, if n is an odd number, then 
Z(n) =m<n—1 <n; Ifn+1 is an odd number, then Z(n +1) =m <n. In any cases, we 
have 


m<n. (2) 


Combining (1) and (2) we have n < m < n, it is not possible. This proves Theorem 1. 
Now we prove Theorem 2. For any positive integer M, we taking positive integer a@ such 
that s = 2° > M+1. This time we have 


Za) = 22") = 3°" = 1, 
Since s + 1 is an odd number, so we have 
Z(st+1)<s=2". 
Therefore, we have 
|[Z(s) — Z(s+1) | > (2°? -—1) -—-2% = 2% -1>M+1-1=M. 


So there exists a positive integer s such that the absolute value of Z(s) — Z(s + 1) is greater 
than M. This completes the proof of Theorem 2. 
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Abstract For any positive integer n, the famous F.Smarandache function S(n) is defined 
as the smallest positive integer m such that n divides m!. The main purpose of this pa- 


per is using the elementary method to study the number of all positive integer n such that 
S(2)-S(4)-S(6)- --- S(2n) 
S(1) -$(3)-S(5)- +--+ S(2n-1) 


is a positive integer. 


Keywords F.Smarandache function S(n), related problem, positive integer solution. 


§1. Introduction and Result 


For any positive integer n, the famous F.Smarandache function S(n) is defined as the 
smallest positive integer m such that n divides m!. That is, S(n) = min{m: me N, n|m}}, 
where N denotes the set of all positive integers. From the definition of S(n), it is easy to see 
that if n = p{''p$?---p?* be the factorization of n into prime powers, then we have 

olny = me OE 
It is clear that from this properties we can get the value of S(n), the first few values of Sn) 
are S(1) = 1, S(2) = 2, S(3) =3, S(4) = 4, S(5) = 5, S(6) = 3, S(7) =7, S(8) = 4, S(9) =6, 
S(10) =5,--++-: . About the arithmetical properties of Sn), some authors had studied it, and 
obtained some interesting results. For example, Farris Mark and Mitchell Patrick [1] studied 
the bound of $(n), and got the upper and lower bound estimates for S(p*). They proved that 


(p—- Del = Ste?) <(e— 1) [e+ 1 log, a +1: 


Lu Yaming [2] studied the solutions of an equation involving the F.Smarandache function S(n), 
and proved that for any positive integer k > 2, the equation 


S(my + me +--+ + mp) = S(m1) + S(me2) +--+ + S'(mz) 


has infinite positive integer solutions (m1, m2, --- , mx). 
Jozsef Sandor [3] proved that for any positive integer k > 2, there exist infinite group 
positive integers (m ,, m2, +: , Mx) satisfying the inequality: 


S(mi + Meo +--++ mx) > S(mz) + S(mz2) + +++ + S(mg). 
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Also, there exist infinite group positive integers (m1, m2, ---, mx) such that 
S(my + me +-+-+mg) < S(m1) + S(me2) +--+ + S(me). 


In [4], Fu Jing proved a more general conclusion. That is, if the positive integer k and m 
satisfying one of the following conditions: 

(a) k > 2 and m > 1 are odd numbers. 

(b) k > 5 is odd, m > 2 is even. 

(c) Any even number k > 4 and any positive integer m; 

then the equation 


m-S(my +m +--+ + mpg) = S(m1) + S(m2) +--+ + S(me) 


have infinite group positive integer solutions (m1, ma, ---, Mx). 

Xu Zhefeng [5] studied the value distribution properties of S(n), and obtained a deeply 
result. That is, he proved the following Theorem: 

Let P(n) denotes the largest prime factor of n. Then for any real number x > 1, we have 
the asymptotic formula 


2¢ (3) «2 v2 
= P(q))* = ee aoa 
(Stn) ~ Peay)? = “SE +0 | ee). 
where ¢(s) is the Riemann zeta-function. 

On the other hand, in the manuscript “Problems lists for collective book on Smarandache 
notions”, Kenichiro Kashihara proposed the following problem: Find all positive integer n € N 
such that 

S(2)-S(4)-S(6)- --- S(2n) (1) 
S(1)-$(3)-S(5)- --- S(2n—1) 


is an integer. About this problem, it seems that none had studied it yet, at least we have 
not seen related papers before. In this paper, we using the elementary method to study this 
problem, and solved it completely. We shall prove the following conclusion: 


Theorem. For any positive integer n, the formula 


S(2) -S(4) -S(6)- --- $(2n) 
SG) -5(3)- 5()- --- SQn—1) 


is an integer if and only ifn = 1. 


§2. Proof of the theorem 


In this section, we shall use the elementary method to complete the proof of the theorem. 
First we need the following two simple lemmas. 

Lemma 1. For any positive integer n > 5, there exists at least one prime P such that 
Peé(n, 2n-1). 

Proof. See Theorem 5.7.1 of reference [6]. 
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Lemma 2. Let p be a prime. Then for any positive integer k, we have the estimate 
S(p*) < kp. If k < p, then S(p*) = kp. 

Proof. See reference [1]. 

Now we use these two lemmas to complete the proof of our theorem. It is clear that if 


n= 1, then 
S(2)-S(4)-S(6)- ++. S(2n) 


5(1)- 5) -S(5)- -- SQn—-1) 


is an integer. In fact, this time we have 


$(2)- $(4)-$(6)- +++ S@n)__ S(2)_, 
Si) S6)-SO).-+Sn—1 say 


If n = 2, we have S(1) = 1, S(2) = 2, S(3) = 3, S(4) = 4, obviously formula (1) is not an 
integer. 

Similarly, if n = 3 and 4, we have S(1) = 1, $(2) = 2, S(3) = 3, $(4) = 4, S(5) = 5, 
S(6) = 3, S(7) =7, S(8) =4, obviously (1) is not also an integer. 

Now we assume that n > 4. From Lemma 1 we know that there exists at least one prime 
P €(n, 2n—1] such that S(P) = P. For this prime P, we have P | S(1)-S(3)-S(5)---- S(2n—1). 
But we can prove that 


P + S(2)-S(4)-8(6)» --- S(2n). 


Otherwise, there exists an integer k with 1 < k <n such that P | S(2k). Let S(2k) = aP. If 
a = 1, then S(2k) = P and P | 2k. So P| k. Therefore, k > P >n+1. This is a contradiction 
with 1<k<n. If a> 2, then from Lemma 2 we have 2k > $(2k) = aP > 2P > 2(n +1), or 
k>P>n-+1, contradict withl1<k<n. 

This completes the proof of our theorem. 
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Abstract The concept of generalized open sets in generalized topological spaces was intro- 
duced by A. Csaszar [7,8]. In this paper, we introduce a class of topological spaces called 
y-s-closed spaces by utilizing the y*-semi-closure. In [16] Maio and Noiri showed that cd- 
compactness due to Carnahan [3], weak RS-compactness due to Hong [11] and s-closedness 
are all equivalent [16]. It is shown that the concept of 7-s-closed spaces generalized s-closed 
spaces [17]. It is interesting to note that y-s-closedness is the generalization of yo-compactness 
defined and investigated in [2]. 

Keywords 4-closed (open), y-closure, y*-semi-open (closed), 7*-semi-closure, y*-semi-interior, 
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81. Introduction 


A. Csaszar [7,8] defined generalized open sets in generalized topological spaces. In 1979, 
S. Kasahara [12] defined an operation a on topological spaces. In 1992 (1993), B. Ahmad 
and F.U. Rehman [1,19] introduced the notions of y-interior, y-boundary and y-exterior points 
in topological spaces. They also studied properties and characterizations of (7, 3)-continuous 
mappings introduced by H. Ogata [18]. 

In 2005, A. Guldurdek and O.B. Ozbakir [10] defined and discussed 7-semi-open sets using 
y-open sets in topological spaces which are different from the notions of y-open sets introduced 
and studied by H.Ogata [17] in 1991. In 2006, S.Hussain, B. Ahmad and T. Noiri [6] defined 
and discussed y*-semi-open sets, y*-semi-closed sets, y*-semi-closure sets, 7*-semi-interior sets 
in a space X in the sense of H. Ogata [18]. 

In 1987, G. Di Maio and T. Noiri [16] introduced the notion of s-closed spaces. It was 
shown that cd-compactness due to Carnahan [3], weak RS-compactness due to Hong [11] and 
s-closedness are all equivalent [16]. 

In this paper, we introduce a class of topological spaces called 7-s-closed spaces by utilizing 
the y*-semi-closure. It is shown that the concept of y-s-closed spaces generalized s-closed spaces 
[17]. It is interesting to note that 7-s-closedness is the generalization of yo-compactness defined 
and investigated in [2]. 


88 S.Hussain and B.Ahmad No. 4 


Hereafter, we shall write spaces in place of topological spaces on which no separation 
axiom is assumed explicitly. We recall some definitions and results used in this paper to make 
it self-contained. 


§2. Preliminaries 


Definition 2.1.[12] Let (X,7) be a space. An operation y: 7 — P(X) is a function from 
T to the power set of X such that V C V7, for each V € 7, where V7 denotes the value of y at 
V. The operations defined by 7(G) = G, 7(G) = cl(G) and 7(G) = int cl(G) are examples of 
operation ¥. 

Definition 2.2.[18] Let A be a subset of a space X. A point x € A is said to be a 7-interior 
point of A if there exists an open nbd N of x such that NY C A and we denote the set of all 
such points by int,(A). Thus 


int,(A)={xeA:xeNer and N?CA}CA. 


Note that A is y-open [18] iff A = int,(A). 

A set A is called y- closed [18] iff X-A is y-open. 

Definition 2.3.[12] A point x € X is called a y-closure point of A C X, if UNA F 4, 
for each open nbd U of x. The set of all y-closure points of A is called y-closure of A and 
is denoted by cl,(A). A subset A of X is called y-closed, if cl,(A) C A. Note that cl,(A) is 
contained in every y-closed superset of A. 

Definition 2.4.[19] An operation y on 7 is said to be regular, if for any open nbds U, V 
of x € X, there exists an open nbd W of x such that UY 1 V7 5 W7. 

Definition 2.5.[18] An operation + on 7 is said to be open, if for every nbd U of a € X, 
there exists a y-open set B such that « € B and U7 D B. 

Definition 2.6.[6] A subset A of a space (X,7) is said to be a y*-semi-open set, if there 
exists a y-open set O such that O C A C cl,(O). The set of all y*-semi-open sets is denoted 
by SO,»(X). A is 7*-semi-closed iff X — A is y*-semi-open in X. Note that A is 7*-semi-closed 
if and only if int,(cl,(A)) C A. 

We denote ['(X), the set of all monotone operators. 

Definition 2.7.[6] Let A be a subset of a space X and y € ['(X). The intersection of 
all y*-semi-closed sets containing A is called 7*-semi-closure of A and is denoted by scl. (A). 
Note that A is y*-semi-closed if and only if scl,.«(A) = A. 

Definition 2.8.[6] Let A be a subset of a space X and y € [(X). The union of y*-semi- 
open subsets of A is called y*-semi-interior of A and is denoted by sint,«(A). 


§3. 7*-Semi-Regular Sets 


Definition 3.1. A subset A of a space X is called y-regular-open, if A = int,(cl,(A)). 
The set of y-regular open sets is denoted by RO,(X,7T). 
Note that RO,(X,7) C 1, Cr. 
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Definition 3.2. A set A is y-regular closed, denoted by RC,(X,7), if one of the following 
conditions holds: 
(i)A = cl,(int,(A)). 


(ii)X — A€ RO,(X,7). 


Clearly A is y-regular open if and only if X — A is 7-regular closed. 

Definition 3.3./6] A subset A of a space X is said to be y*-semi-regular, if it is both y*- 
semi-open and y*-semi-closed.The class of all y*-semi-regular sets of X is denoted by SR,+(A). 

Note that if 7 is regular operation, then the union of y*-semi-regular sets is y*-semi-regular. 

Definition 3.4. A subset A of X is said to be y*-semi-regular-open if A = sint.,«(scl,+(A)). 

Lemma 3.5. If A is a subset of a space X, then int+(cl,(A)) C scl,-(A). 

Proof. Let x € int,(cl,(A)). Let G be a 7*-semi-open set of Xcontaining x. Then 
U CG Ccl,(U), for some y-open set U in X. 

Since x € G Ccl,(U) and a € int,(cl,(A)), we have 


 # int, (cl,(A)) NU C el,(A) NU Cel,(ANU). (byLemma2(3)[19]) 


Therefore, we have ANU # ¢ and hence ANG ¥ ¢. This shows that x € scl+(A). 

Hence the proof. 

We use Lemma 3.5 and prove: 

Proposition 3.6. For any subset A of a space X, the following are equivalent: 

(1)A € SRy+(X). 

(2).A = sint.,« (scl,«(A)). 

(3) There exists a y-semi-regular-open set U of X such that U C A C cl,(U), where 7 is 
an open operation. 

Proof. (1)=(2). If Ac SR,-(X), then sint.,~(scly«(q)) = siny«(A) = A. 

(2)=(3). Suppose that A = sint,-(scly«(A)). By Lemma 3.5, for any subset A of X, 
inty(cl,(A)) C scl,«(A) implies 


int,,(cl,(A)) C sint,+(scl,«(A)) = A. 
Since A € SO,«(X), we have A C cl,,(int,(A)). Therefore, we obtain 
int+(cly(A)) € A € cl,(int,(A)) € cl,(int,(cl,(A))), 


where int.(cl,(A)) is y-regular open, since int,(cl,(int,(cl,(A)))) = int,(cl,(A)). 

(3)=(1). It is obvious that A € SO,+(X). Clearly from (3) and Theorem 3.6[18], we have 
int.,,(cl,(A)) = int,(cl,(U)) = U C A and hence A is y*-semi-closed. 

Thus, we obtain A € SR,«(X). This completes the proof. 

Proposition 3.7. If A € SO,«(X), then scl,-(A) € SR,-(X). 

Proof. Since scl,-(A) is y*-semi-closed, we show that scl+(A) € SO,(X). Since A € 
SO,«(X), then for y-open set U of X, U C A Cl,U. 

Therefore we have 


U C scly«(U) C sely«(A) C sely+(cl,(U)) = ely(U) 
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or 
U C scly«(A) C ely(U) 


and hence scl,«(A) € SO,(X). Hence the proof. 

Definition 3.8. A point x € X is said to be a y-semi-0-adherent point of a subset A of 
X if scly«(U) NAF ¢, for every U € SO,+(X). The set of all y-semi-6-adherent points of A is 
called the y-semi-6-closure of A and is denoted by s+clg(A). 

A subset A is called y-semi-6-closed if sclg(A) = A. 

Proposition 3.9. Let A be a subset of a space X. Then we have 

(1) If A € SO,+(X), then scl,+(A) = s,clg(A). 

(2) If Ae SR,-(X), then A is 7-semi-6-closed. 

Proof. (1) Clearly scl,«(A) C s,clg(A). Suppose that x ¢ scly«(A). Then, for some 
y*-semi-open set U, ANU = ¢ and hence AN sel,«(U) = ¢, since A € SO,«(X). This shows 
that x ¢ syclg(A). Therefore scl»«(A) = syclo(A). 

(2) This follows from (1). Hence the proof. 

Lemma 3.10. If B is 7-open set in a space X, then scl,«(B) = int,(cl,(B)). 

Proof. From Lemma 3.5, we have int,(cl,(B)) C scly«(B), for any subset B of X. 

Now we show that scl,«(B) C int,(cl,(B)). For this let x ¢ int,(cl,(B)). Then 


x € cl,(int,(X — B)) € 504+ (X). 


Since B is y-open, we have B C int,(cl,(B)) and BN cl,(inty(X — B)) = ¢. This shows that 
x ¢ scl,«(B). Therefore, we obtain scl,+(B) = int,(cl,(B)). 

Hence the proof. 

Definition 3.11. A space X is y-extremally disconnected space, if cl,(U) is y-open set, 
for every y-open set U in X. 

Lemma 3.12. A subset A of a space X is 7*-semi-open if and only if cl,(A) = cl,(int,(A)), 
where y is an open operation. 

Proof. Suppose A is y*-semi-open. Then we have A C cl,(int,(A)) and since y is an 
open operation, so cl,(A) C cl,(int,(A))[19]. 

On the other hand, we have int,(A) C A and hence cl,(int,(A)) C cl,(A). Consequently, 
we have cl,(A) = cl,(int,(A)). 

Conversely, by the supposition, we have 


inty(A) C A C cl,(A) = cl,(int,(A)). 


Hence A is y*-semi-open. 

Hence the proof. 

Proposition 3.13. A space X is y-extremally disconnected if and only if cl,(U) = 
scl,«(U), for every U C SO,«(X), where y is a regular and open operation. 

Proof. Necessity. In general scl,«(U) € cl,(U), for every subset U of X. We show that 
cl,(U) C scly«(U), for each U € SO,«(X). Let 6 4 scl,+(U) € SOy+(X) and let x ¢ scl,«(U), 
then there exists V € SO,«(X) such that x € V and VNU = ¢ and hence int,(V) Nint,(U) = 
a{19]. 
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Since X is y-extremally disconnected, we have 
cl, (int,(V)) N cl, (int,(U)) = ¢. 


Therefore, by Lemma 3.12, we have « ¢ cl,(int,(U)) = cl,(U). 
Sufficiency. For every y-open set U, U € SO,y«(X) and by Lemma 3.10 , 


cl,(U) = scl,«(U) = int,(cl,(U)). 


This shows that cl,(U) is y-open [18], for every U € Ty. 


§4. y-s-closed spaces 


Definition 4.1. A space X is said to be y-s-closed if for every cover {V,, : a € I} of X by 


*-semi-open sets of X, there exists a finite subset Jp of J such that X = U sclyx (Va). 
a€lo 
Definition 4.2. A filter base IT on X is said to y-SR-converges to x € X if for each 


V € SR,-(X), there exists F €T such that F CV. 

A filter base T is said to be 7-SR-accumulate at « € X if VNF ¥ ¢, for every V € SRy+(X) 
and every FET. 

Proposition 4.3. For any space X, the following are equivalent: 

(1) X is y-s-closed. 

(2) Every cover of X by y*-semi-regular sets has a finite sub cover. 

(3) For every family {V, : a € I} of y*-semi-regular sets such that (\){Va : a € I} = 4, 
there exists a finite subset Jp of I such that (\{Va : a € Ip} = ¢. 

(4) Every filter base 7-SR-accumulates at some point of X. 

(5) Every maximal filter base y-SR-converges to some point of X. 

Proof. (1)=(2). Follows directly from the definition of y-s-closed space. 

(2)=(5). Let T be the maximal filter base on X. Suppose that I’ does not y-SR-converge 
to any point of X. Then T does not y-SR-accumulate at any point of X. 

Therefore for each x € X, there exist F, € T and V, € SR,+(X) and such that V;N Fy, = ¢. 
The family {V,,: « € X} is a cover of X by y*-semi-regular sets of X. By (2), there exists finite 
number of points 71, %2,%3,-++ ,% such that X = U{Vz, :4=1,2,--- ,n}. 

Since [ is filter base on X, there exists Fo € I’ such that Fo C (\{Vz, : 7 = 1,2,--- ,n}. 
Therefore, we have Fo = ¢. 

This contradiction proves the required. 

(5)=>(4). Let T be a filter base on X and Ip a maximal filter base such that [ CT. By 
(5), Io y-SR-converges to some « € X. For every F €T and every V € SR,-(X), there exists 
Fo € To such that Fo C V. Therefore, we obtain VN F D> Fo NF 4 ¢. This shows that [ 
y-SR-accumulates at x. 

(4)=(8). Let {Va : a € I} be a family of y*-semi regular sets such that (\{Vq :a€ I} = ¢. 
Let A(J) denot the family of all finite subsets of J. Assume that (\{V. : a € Q} ¥ 4, for every 
Qe A(Z). 
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Then the family [ = () Va: Q€T(L) ? isa filter base on X. By (4), [ y-SR-accumulates 
aeQ 
at some x € X. Since {X —Vy: a € I} isa cover of X, x € X—V,,, for some ag € I. Therefore, 


we have X — Vy, € SR,-(X) and V,, € IT. This is a contradiction. 

(3)}(1). Let {V, € I} be a cover of X by y*-semi-open sets of X. By Proposition 3.7, 
{scl«(Vq) : a € I} is a y*-semi-regular cover of X. 

Thus {X — scly-«(V.) : a € I} is a family of 7*-semi-regular sets of X having the empty 
intersection. 

By (8), there exists a finite subset I, of J such that 


(\{x — scly«(Va): a€Ibs=¢ 


and hence X = U{scl,«(Va € Io)}- 

This shows that X is y-s-closed. This completes the proof. 

Corollary 4.4. For any space X, the following are equivalent: 

(1) X is 7-s-closed. 

(2) Every y*-semi-open and 7*-semi-closed cover of X has a finite sub cover. 

(3) Every y*-regular-semi-open cover of X has a finite sub cover, where 7 is an open 
operation. 


Proof. This is an immediate consequence of Proposition 3.6 and Proposition 4.3. 
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Abstract In this paper, we obtain a successive linear programming algorithm for solving the 
semidefinite programming (SDP) relaxation of the binary quadratic programming. Combining 
with a randomized method of Goemans and Williamson, it provides an efficient approximation 
for the binary quadratic programming. Furthermore, its convergence result is given. At last, 
We report some numerical examples to compare our method with the interior-point method 
on Maxcut problem. 

Keywords Binary quadratic programming, successive quadratic programming algorithm, 


semidefinite programming, randomized method. 


81. Introduction 
In this paper, we consider the following binary quadratic programming 


min «2? Qr+2r?x (1) 
Site, ee Stored = Lae oy 
where Q is ann Xn real symmetric matrices, r is a real n-dimensional column vectors. Without 
loss of generality, we assume that Q is a positive definite matrix, because of the equivalence 


n 
between minge{1,-1390' Qrt 2r'a, and minge{1,-1}9 27 Qrt2r? rt > yi(z? -1) for ally € R. 
i=1 
In mathematical term, it means that Q > 0. 
The binary quadratic programming is a fundamental problem in optimization theory and 
practice. VLSI design, statistical physics, combinatorial optimization, the optimal multiuser 


detection and the design of FIR filters with discrete coefficients are all sources of the quadratic 
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binary programming [1]—[5]. These problems are know to be NP hard [1]. One typical approach 
to solve these problems is to construct lower bounds for approximating the optimal value. The 
classic technique to obtain bounds is either via a continuous relaxation or via the dual problem, 
which is usually followed by branch and bound type algorithms for refining it. e.g.the paper [6] 
and reference therein. Recently, the semidefinite programming relaxation approach had been 
studied and proven to be quite powerful for finding approximate optimal solutions. See, e.g., 
[7] and references therein. But the semidefinite programming approach is limited to problems 
of moderate size, which can’t solve the problems of large scale efficiently [1]. Recently, the 
nonlinear programming methods have been proposed to solve the large scale problem, which 
is more efficient. For instance, Helmberg and Rend1 [7] introduce the spectral bundle method 
which builds on the framework of the proximal method of Kiwiel; The approach by Homer and 
Peinado [8] for using the change of variables X = VV7, V Cc R"*", where X is the primal 
matrix variable of the maxcut SDP relaxation, is to transform the maxcut SDP relaxation into 
a constrained nonlinear programming problem in the new variable V. More recently, Samuel 
Burer and Renato D.C.Monteiro [9] propose a variant of Homer and Peinado’s method based on 
the constrained nonlinear programming reformulation of the maxcut SDP relaxation obtained 
by using the change of variable X = LL7, where L is a lower triangular matrix. 


In this paper, we consider the binary quadratic programming and its corresponding refor- 
mulation of the SDP relaxation directly. A successive linear programming algorithm for solving 
SDP relaxation of the binary quadratic programming is provided by using the SDP relaxation 
and the change of variables X = VV", V Cc R"*". Furthermore, its convergence result is 
given. The step-size in our algorithm is obtained by solving n easy quadratic equations without 
using the linear search technique. The computational experience with our method indicates 
that it is substantially faster than the interior-point method. 


The paper is organized as follows. In Section 2, we present the binary quadratic program- 
ming problem and its relaxations. In Section 3, the successive linear programming algorithm of 
the relaxation problem is obtained and its convergence result is given. Some numerical examples 
are offered in the last section. 


Notation and Terminology. 


In this paper, #, R",and R"*”" denote the space of real numbers, real n-dimensional 
column vectors, and real n x n matrices, respectively. By S” we denote the space of real n x n 
symmetric matrices, and we definite St and St, to be the subsets of S” consisting of the 
positive semidefinite and positive definite matrices respectively. We write A = 0 and A > 0 to 
indicate that A € Sand A € S", respectively. We let tr(A) denote the trace of a matrices 
Ae R"*", we defined Ae B = (A, B) = tr(A’B), and the Frobenius norm of A € R"*” is 
defined to be || A ||p= (Ae A)!/?. We adopt the convention of denoting matrices by capital 
letters and matrices entries by lowercase letters with double subscripts. For example, a matrix 
A € ®”*” has entries aj; for 7,7 = 1,--- ,n, In addition, we denote the rows of a matrix by 
lowercase letters with single subscripts. For example, A € #”*%” has rows a; for i = 1,--- ,n. 
In this paper, we will often find it necessary to compute the dot product of two row vectors 
a; and b; which arise as rows of the matrices A and B. Instead of denoting this dot product 


as ab}, we will denote it as (a;,b;). 
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§2. Binary quadratic programming and its relaxations 


In this section, we consider the following binary quadratic programming and describe some 
of its relaxations. 
max 2? Cx (1) 
Sy a =I for = 12m: 
Without loss of generality, we assume that C is a positive definite matrix, because of the 
equivalence between Maxze(1,-1}» 2" Cr and Maxgef1,-1}>27 Cr+ > yi(x? —1) for ally € &. 
i=l 


In mathematical term, it means that C' > 0. 
Let X = xx7, the above problem equals the following problem, 


max CexX 

st. (eef)eX =1, fori=1,---,n (2) 
rank(X) = 1 
X > 0, 


where e; is the unit vector whose i-th component is 1 and others are all 0. 
Because rank one constraint is nonconvex, dropping rank one constraint yields a semidefi- 


nite programming relaxation of (1) as follows, 


max CeX 
st. (eef)eX =1, fori=1,---,n (3) 
X= 0. 


Now, we consider the following problem, 


max CeX 
st. (ee? )eX <1,fori=l1,---,n (4) 
X > 0. 


Obviously, any feasible solution to the problem (3) is a feasible solution to the problem 
(4), say, the problem (4) is a relaxation of the problem (3). Conversely, we assumed that 
X is an optimal solution to the problem (4) and let Y be a matrix which is satisfied that 
yi = 1,fori =1,--- ,n and yy = vi3,1 € j,for i,j =1,---,n. By using the C > 0 and the 
fact that 7; <1 fori =1,---,n, we have Ce X < CeY. Along with that Y is a feasible 
solution of the problem (3), the optimal values of the problem (3) and (4) coincide. This shows 
the problem (3) and (4) are equivalent. 

We now present the nonlinear programming reformulation of the problem (4) which is the 
basis of our algorithm for finding an approximate solution of the binary quadratic programming. 
For every X € S", there exists a matrix V € R”*” such that X = VV‘. Thus the problem 
can be stated as the following one, 


max Ce(VV7) 
st. (ee?)e(VV7) <1,fori=1,--- ,n (5) 
Yer, 
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Notice that we have replaced the requirement X > 0 with X = VV7,V € R"*”. So the 
objective function of the problem (5) is nonconvex, but the feasible set of the problem (5) is 


convex. 


§3. The successive linear programming algorithm to the 


relaxation problem 


In this section, we develop and discuss the successive linear programming algorithm to solve 
the problem (5). Before giving the basic steps of the algorithm, we definite some functions as 
follows. 

f RGR, f(V) =Ce(VV") 


gi RX" 1 RK, gi(V) = ese? @ (VVT) —1,i=1,...,n 


Obviously, the gradient of function f(V) at a point V is G = 2CV, the gradient of function 
gi(V) at a point V is H; = 2(ee7)V,i =1,--- ,n. 

Given a matrix V* feasible for the problem (5), the feasible ascent direction D* of the 
function f(V) at a point V* will be obtained by solving the following quadratic programming. 


max ¢ 
st. —GreD+t<0 (6) 
(e;e7 ) e (ve(V™)*) —1+ 2(eve? )V* eD+t < O, for t= 1, nr yn, 


where (t,D) € Rx R"*”. 

Proposition 3.1. Given (t*, D*) € R x R"*” is a optimal solution for problem (6), and if 
D* #0, then D* is the feasible ascent direction D* of the function f(V) at a point V*. 

Proof. Since t = 0,D = 0 is a feasible solution for problem (6), (t*, D*) is a optimal 
solution for problem (6), and D* 4 0, so we have t* — % DF e D’ > 0, that is to say t* > 
%D* e D* > 0. Furthermore,since —G* e D* + t* < 0, we have G¥ e D* > t* > 0. From the 
assumption D* 4 0, we obtain that D* is the ascent direction of function f(V) at a point V*. 
Next, we will prove D” is the feasible direction. 

Let 6; > 0,1 =1,--- ,n, and we have 

gi(V* + 6,D*) = (ee7’) @ (V* + 6,D*)(V* + 6,D*)? -1 


= (ere?) e V*(V*)? — 1 + 26;(e;e7 )V* @ D* + 6?(e;e7) @ D¥(D*)?. 
If (e;e7 )V* @ D* > 0, let 6; < 1, and satisfy 6?(e;e7) e D¥(D*)? < t*, We obtain 
gi(V" + 5:D*) < (exe?) @ (VE(V*)) —14+2(ee7)V¥ oe D+t <0. 
If (exe? )V* e D¥ < 0, let 5; satisfy 26;(e;e2 )V* e D* + 5?(e;e2) e D¥(D*)* < 0, we obtain 


—2(e,e7)V* e DP 


2. ‘ 
O< OES Cee) DDI 
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We have gi(V* + 6;D*) < 0. Based on the analysis, we select 6 = min {6;,i=1,--- , n}, 
which must satisfy 
gi(V* + 6D*) <0. 


Furthermore, since the feasible region of problem (5)is convex, so D* is the feasible ascent 
direction D* of the function f(V) at a point V*. 

By proposition 3.1, we obtain the feasible ascent direction D* of the function f(V) at a 
point V*. Now, we will give a simple method to obtain an appropriate step-size. 

From proposition 3.1, we know that if D* 4 0, there is a step-size 5 > 0, which satisfy 
gi(V* + 5D*) <0, for i=1,---,n +1. 

Given 6; > 0, by solving the equations g;(V* + 6D*) = 0, for i=1,--- ,n +1. We obtain 


(ese VE «DE + Vee WE DEF (eel 0 VIF — (ese! J ® DUD) 
(e;e)) e Dk( Dk)? : 


5, = 


We choose the step-size 
6 = min {6;,7 = 1,--- , n}. (7) 


Proposition 3.2. Given D* ¥ 0, the step-size from (8) is an appropriate step-size. 
Proof. Given scalar h > 0, We define y: RH RK by yv(h) = f(V* +hD*), and we have 


y(h) = f(V* + AD*) = f(V*) + Atr(G*D*) + h?C e D*(D*)?. 
y'(h) = tr(G* D*) + 2hC e G*(G*). 


Based on the proposition 3.1 and D* 4 0, we have tr(G* D*) > 0. Furthermore, since C' > 0, so 
CeG*(G*) > 0, Thus, when h > 0,y'(h) > 0, then y(h) is a monotone increasing function. 
Directly following from (7), the step-size is an appropriate search step-size. 

We are now ready to write the successive linear programming algorithm. 

The successive linear programming algorithm. 

Let V° be a feasible solution of the problem (5), which is satisfied that (e;e7)e(V°(V°)7) = 
1, fori =1,---,n. And a prespecified constant ¢ > 0. Let u = 0.5/n. 

1. Compute the gradient G* = 2CV* for the function f at the point V*. 

2. Compute the feasible ascent direction D* by solving problem (6). By the proposition 
3.1, we obtain D*. 

3. Compute||D*|| p. 

4. If ||D*||r < e, then stop; otherwise,go to 5. 

5. compute 6 by formula (7). Let V#t! = V*¥ + 6D*,k =k +1, go to 1. 

Now,we will prove the convergence of the above algorithm. 

Proposition 3.3. Assume that (t*, D*) € Rx R™*™ is the optimal solution of the problem 
(7), If t* = 0, there are some multipliers *, which satisfy 


207" =2 > pk (2e,e, )V* 
i=1 
pk (ese? ) e V*(V*)? — 1) = 0, and p* > 0 for i=1,--- ,n. 


That is to say the matrix V* is the KKT point of the problem (5). 
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Proof. The dual programming of the linear programming (6) is [1° 


max Se (1 (er 1€7_,) e VE(V*)*) 


(DLP) i 4 (8) 
\,G* + > A;(2e;_1e7_,)V* =0 


1=2 


M>O=1,-+,m41. 


Because t* = 0, due to the dual theory of the linear programming, we have that there is 
an solution A* € R™*+ for problem (8), which satisfy 


m+1 
De Mies ex_1e;_,) @ V*(V*)") =0, A, > 0, i=1,--- m+], (9) 
m+1 
MG* + S$" MF (2e:_1e7_1)V* =0, (10) 
4=2 


If \F £0, let 


By the equation (10), we have 


Ge + S- pk (Qe;e7 )V* = 20V* + S- pk (2e,e7 )V* = 0. 
i=1 i=1 
Based on the (9), we obtained that 
pi (exe, e VE(V*)T — 1) = 0, pe 2 0,1 = 1,-°: ;m+ 1. 


If \¥ = 0, following from (10), we have 


m+1 


S> AF (Qe-1eF ,)V* =0. 


Thus, we have 


m+1 
S> AF (2e_1€7-1)V* e V* =0. 
Following from (9), we have 
m+1 m+1 m+1 
So AMF - es eV) = N= a 


1=2 


m+1 

Because \* is the solution of problem (8), then we have )> A¥ = 140, which contradict the 
i=1 

equation above, so \¥ 4 0. We have proven the theory. 
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Based on the proposition 3.3, the algorithm will terminated at the KKT point of problem 
(5) if the termination criterion is t* < ¢. Where ¢ is a constant which is enough small. 

Not every KKT point of the problem (5) is a global solution. The following proposition 
give sufficient conditions for a KKT point of the problem to be a global solution. 

Proposition 3.4. Assume that V* is a KKT point of problem (5) and Let 


ji 0, fori =1,--- ,n;S* = Suk (ee?) — C. 
i=1 


If S* = 0, then V*(V*)? is a global solution to the problem (4), 
Proof. Since V" is a KKT point, together with Proposition 3.4, we have 


Ge = 2cv" -_ So ub (ee? V*, ut 2 0, for i= 1, voy. 


i=1 
Thus 
(CVE,V*) = () Tai (eies V",V") = Dai. 
i=1 i=1 
That is (C, V¥(V*)7) = >> wk. If S* = 0, let X be feasible for (4), then (S*, X) > 0, that is 
i=1 


'S pk > (C,X). Therefore, V*(V*)* is a global solution of the problem (4). 


Il 
un 


84. Numerical results 


In this section we present computational results by comparing our method with earlier 
method to find approximate solutions to the maxcut problem based on solving its SDP relax- 
ation. As stated in the introduction, the purpose of the results presented here are to show that 
our successive linear programming algorithm is considerably faster than interior-point method. 

Maxcut problem. 

The maxcut problem is one of the standard N P—complete problems defined on graphs!!"). 
Let G = (V, FE) denote an edge-weighted undirected graph without loops or multiple edges. We 
use V = {1,--- n}, ij for an edge with endpoints i and j, and a,,; for the weight of an edge 
ij € BE. 

For S C V the cut 6(S) is the set of edges ij € E that have one endpoint in S and the 
other endpoint in V \ S. The maxcut problem asks for the cut maximizing the sum of the 
weights of its edges. Here, we only work with the complete graph K,,. In order to model a 
graph in this setting, define a;; = 0 for ij ¢ E. A = (a;;) € S” is referred to as the weighted 
adjacency matrix of the graph. An algebraic formulation can be obtained by introducing cut 
vectors x € {—1,1}" with x; =1fori¢€ S and x; =—-1 forieV\S. 

The maxcut problem can be formulated as the integer quadratic program. 


<j 


st. 2, €{-1,1},2=1,---,n. 
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The matrix L(G) = Diag(Ae) — A is called the Laplace matrix of the graph G. Where e is the 
unit vector whose every component is 1. And Diag(Ae) is the diagonal matrix whose diagonal 
elements are Ae. Let C = ‘L, the maxcut problem may be interpreted as a special case of the 
problem (1). 


Numerical examples. 


We report the numerical example in this section. In the numerical example, we compare 
the computational results between our method and interior point method [5]. 


As stated before, the purpose of the results presented here is to show that our algorithm 
is substantially faster than interior point method. All the algorithms are run in the MATLAB 
6.1 environment on a AMD AthlonXP1600+ personal computer with 128Mb of Ram. 


In all the test problems, we choose the initial iterate L° to be n x n identity matrix, u to 
be 0.5/n. In interior-point method, we solve the SDP relaxation by using SDPpack software 
[12]. We adopt the randomized cut generation scheme of Goemans and Williamson [4]. The 
iteration number of this algorithm is n which is the scale of max-cut problem considered. Here 
all the tested problems are random graphs with two different edge density 0.7 and 0.3, which 
denote the dense random graphs and sparse random graphs respectively. 


We select problems in size from n = 50 to n = 300 to compare the suboptimal value of 
maxcut problem,and the total time of the three methods. In our algorithm, the iteration stops 
once t* < € is found. The result is shown in Table.1. 


In Table 1, we use “SDP” presents for interior point algorithm based on semidefinite 
programming, “SLA” for our successive linear programming algorithm, “time” for the total 
time of two methods, “value-f” for the suboptimal value of the maxcut problem based these 
methods, “density” for edge density of the random graphs. The Table 1 shows our method can 
generally reach solutions of the problems much faster than the interior-point method whether 
to the dense random graphs or the sparse random graphs. 


Table.1 Comparison of the two method 
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Size of graph | density | algo value-f time 
50 0.7 SLA | 381.316300729441 0.1720 
SDP | 381.222818884658 1.2970 
50 0.3 SLA | 384.343853121280 0.4220 
SDP | 384.343853121280 3.2190 
100 0.7 SLA | 1385.94015100020 1.4530 
SDP | 1388.66226246641 | 13.7190 
100 0.3 SLA | 7179.35285955398 | 10.0000 
SDP | 7184.64381700424 | 149.2200 
150 0.7 SLA | 3027.47658904186 3.5470 
SDP | 3030.60112063104 | 75.6090 
150 0.3 SLA | 1527.69134585512 3.7500 
SDP | 1528.54917463768 | 82.5470 
200 0.7 SLA | 5323.36289231888 | 18.3750 
SDP | 5318.11148873339 | 184.5470 
200 0.3 SLA | 2676.47956117645 | 15.1250 
SDP | 2672.86837035170 | 200.8750 
250 0.7 SLA | 8156.16895709477 | 45.4850 
SDP | 8162.65641905066 | 440.6560 
250 0.3 SLA | 4129.85944356911 | 14.8750 
SDP | 4135.51356364934 | 480.8120 
300 0.7 SQA | 11686.0672820264 | 100.1720 
SDP | 11697.2745314466 | 2116.0620 
300 0.3 SLA | 5884.30709295320 | 80.6100 
SDP | 5889.15827541389 | 1892.8590 


§5. Conclusion 


In this paper, we have proposed a successive linear programming algorithm for solving SDP 
relaxation of the binary quadratic programming. In the algorithm, we give a simple method for 
selecting the step-size. Using the randomized cut procedure of Goeman and Williamson, it can 
give a sub-optimal of max-cut problem. It is able to obtain a moderately accurate solution more 
quickly than interior point method. This paper has demostrated the single case of max-cut SDP 


relaxation, but we believe that the same results are apt to hold elsewhere. 


Vol. 3 A successive linear programming algorithm for SDP relaxation of binary quadratic programming 103 


References 


1] F. Barahon and M. Grotschel, An application of combinatorial optimization to statiscal 
optimization and circuit layout design, Operation Research, 36(1998), No. 3, 493-513. 

2| K. C. Chang and D. H., Eifficient algorithms for layer assignment problems, IEEE 
Trans, on Computer Aided Design, CAD-6(1987), No. 1, 67-78. 

3] R. Chen, A graph theoretic via minimization algorithm for teo layer printed circuit 
boards, IEEE Trans, on Circuit and System, CAS-30(1983), No. 5, 284-299. 

4] M. X. Goeman and D. P, Williamson, Improved approximation algorithms for maximum 


cut and satisfiably problem using semidefinite programming, Journal of ACM., 42(1995), 1115- 
1145. 
5] S. Benson, Y. Ye and X. Zhang, Solving large scale sprase semidefinite programming 
for combinatorial Optimization, SIAM. J. Optim., 10(2000), No. 2, 443-461. 

6] C. Helberg, F. Rendel, An interior point method for semidefinite programming, SIAM 
Journal On Optim., 10(1990), 842-861. 

7| C. Helberg and F. Rendl, A spectral bundle method for semidefinite programing, SIAM 
Journal on Optim., 10(2000), 673-696. 

8] M. Peinadoo and S. Homer, Design and performance of parallel and distributed approx- 
imation algorith, 9(1998), 141-160. 

9| S. Burer and R. D. C. Monteriro, A projected gradient algorithm for solving the max-cut 
relaxation, Optimization methods and Software, 15(2001), 175-200. 

10] Y. X. Yuan, W. Y. Sun, The Optimal theories and methods. Beijing, scientific pub- 
lishing company, 1995. 


11] C. Helmberg, Semidefinite Programming for Combinatorial Optimization, Konrad- 
Zuse-Zentrum fur informationstechnik Berlin, Germany, 2000. 

12] M. V. Nayakakuppam, M. L. Overton, and S. Schemita, SDPpack user’s guide-version 
0.9 Beta, Technical Report Yr 1997-737, Courtant Institute of Mathematical Science, NYU, 
New York, NY, 1997. 


Scientia Magna 
Vol. 3 (2007), No. 4, 104-104 


A note on primes of the form a? + 1 


Juan Lépez Gonzalez 


9, December 2007 
Madrid, Spain 


Abstract In this note I prove using an algebraic identity and Wilson’s Theorem that if a? +1 
is an odd prime, thus this prime must has the form 4k? + 1, then 5 { 2k — 3. 
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If n = a? +1 is prime and n F 2, then n is odd, thus a? is even and n must has the form 


4k? + 1, where k > 1 is an integer. The integers 4k? + 1 can be written as 
Ak? +1 = (2k — 3)? + 3(4k — 3) +1. (1) 


If 2k —-3 = —1 then k = 1, and 5 f —1. If 2k —3 = 1 then k = 2, 17 is a prime with 
(2-2—3,5) =1. If 2k-—3> 1 then 
Ak? +1 = 0+3(2k) +1(mod 2k — 3) 
3(2k — 3) +9 + 1(mod 2k — 3) 
10 (mod 2k — 3). 


By Wilson’s Theorem 
(4k?)! = —1(mod 4k? + 1). (2) 


Thus exists an integer c such that (4k?)! + 1 = c- (4k? + 1), since 4k? > 2k — 3 for all k, then 
2k — 3|(4k?)!, thus 
0+ 1=c-10(mod 2k — 3). (3) 


Then there are integers s and t, such that 
10s + (2k —3)t =1, (4) 


thus (5, 2k — 3) = 1, by contradiction if 5|2k — 3, then 0+ 0 = 1(mod 5). Thus I’ve proved the 
following 
Proposition. If a? + 1 is an odd prime different of 5, then (a — 3,5) = 1. 
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Abstract The natural partial order on an U-semiabundant semigroup is introduced in this 
paper and some properties of U-semiabundant semigroups are investigated by the natural 
partial order. In addition, we also discuss a special class of U-semiabundant semigroups in 


which the natural partial order is compatible with the multiplication. 
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§1. Introduction 


In generalizing regular semigroups, a generalized Green relation LY was introduced by M. 
V. Lawson [4] on a semigroup S as follows: 
Let E be the set of all idempotents of S and U be a subset of EF. For any a,b € S, define 


(a,b) € LY ifandonlyif (VeEeU) (ae=asbe=b); 


(a,b) € RY ifandonly if (VeeU) (ea=aeeb=b). 

It is clear that £C L* C LY andR CR*C RU. 

It is easy to verify that if S is an abundant semigroup and U = E(S) then L* = LY, 
R* = RY": if S is a regular semigroup and U = E(S) then L= LY, R= RY. 

Recall that a semigroup S is called U-semiabundant if each £Y-class and each RY-class 
contains an element from U. 

It is clear that regular semigroups and abundant semigroups are all U-semiabundant semi- 
groups. 

The natural partial order on a regular semigroup was first studied by Nambooripad [7] 
in 1980. Later on, M. V. Lawson [1] in 1987 first introduced the natural partial order on an 
abundant semigroup. The partial orders on various kinds of semigroups have been investigated 
by many authors, for example, H. Mitsch [5], Sussman [6], Abian [8] and Burgess[9]. 
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In this paper, we will introduce the natural partial order on U-semiabundant semigroups 
and describe the properties of such semigroups by using the natural partial order. 

We first cite some basic notions which will be used in this paper. Suppose that e, f are 
elements of E(S). The preorders w” and w! are defined as follows: 


ew" fofe=e and e/fsef=e. 


In addition, w = w" Nw!, the usual ordering on E(S). 

We use Dz to denote the relation (w” Uw!) U (w” Uw!)-!. Assume (9,U) is an U- 
semiabundant semigroup. 

It will be said that U is closed under basic products if e, f € U and (e, f) € Dg then 
ef EU. 

For terminologies and notations not given in this paper, the reader is referred to Howie [3]. 


§2. The natural partial order 


Let S(U) be an U-semiabundant semigroup and a € S. The £LY(RY)- class containing the 
element a will be denoted by LY (RY) respectively. 

We will denote an element of LY MU by a* and an element of RY NU by at. 

Recall in [4] that a right ideal J of a semigroup S is said to be an U-admissible right ideal 
if for every a € I we have RU Cl. 

For a € S, we define the principal U-admissible right ideal containing a, denoted by RY (a), 
to be the intersection of all U-admissible right ideals containing a. Similarly, we may give the 
definitions of an U-admissible left ideal and the principal U-admissible left ideal. 

Let S be a semigroup and x, y € S. We say that RY < Re if RU (x) C RY (y). A partial 
order on the £¥-classes can be defined in the usual left-right dual way. 

Lemma 2.1. RY, < RY, for any elements a and x of S. 

Proof. Clearly, the product az lies in aS, which is the smallest right ideal containing a. 
Since RY (a) is a right ideal containing a, we have aS C RY (a). Thus ax € RY (a). 

It follows immediately that 

RY (ax) C RY (a). 


Lemma 2.2. Let U C E(S) ande, f €U. Then RY < Re if and only if Re < Ry. 

Proof. Suppose first that RY < Re. Then we immediately have RY(e) C RY(f). We 
claim that eS is an U-admissible right ideal. 

In fact, for each a € eS, a = ea and so, for any b € RY, we have b = eb € eS. But RY (e) 
is a right ideal and e € U, so that eS C RY (e). 

Since eS is an U-admissible right ideal, we have that RY(c) = eS. Similarly RY(f) = fS. 
It follows that eS C fS, that is, R(e) C R(f). Hence Re < Ry. 

Conversely, suppose that Re. < Ry. Then eS C fS and so e = fax for some x in S'. Thus, 
by Lemma 2.1, we have RU = Ry < ee 

Corollary 2.3. The following statements hold on an U-semiabundant semigroup S(U) 
for any e, f € U: 
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(i) (e, f) € LY if and only if (e, f) € L; 
(ii) (e, f) € RY if and only if (e, f) ER. 


Theorem 2.4. Let S(U) be an U-semiabundant semigroup such that U is closed under 
basic products. Define two relations on S(U) as follows: 

For any x and y of S(U), 

r<,y if and only if RY < Au and there exists an idempotent rt € RU OU such that 
r=aty; 

xX y if and only if LY < LY and there exists an idempotent «* € LY MU such that 
x= yx. 

Then <,. and <, are respectively two partial orders on S(U) which coincide with w on U. 

Proof. We only need to prove that <,. is a partial order on S(U) which coincides with w 
on U since the proof of <j is similar. 

Reflexivity follows from the fact that S(U) is U-semiabundant. Now suppose that «<,.y 
and y<,a. Then RY = Re and there exist at € RYNU and yt € Re NU such that « = «ty 
and y = y*x. By Corollary 2.3, we have x = (atyt)x = yta = y. Next, suppose that r<,y 
and y<,.z. 

It follows that RY < A < RY and there exist et € RY NU and y* € R, NU such that 
x=atyandy=ytz. Thus x = (xt y*)z and Re, =RU< RO — Ru which gives R,+ < Ry+ 
by Lemma 2.2. 

It follows that x*S(U) C ytS(U) and so xt = ytax*. Since U is closed under basic 
products and (x*, y*) € w” C Dr, we deduce that aty* € U. Clearly, (ct,atyt) € R and so 
(x, xty+) € RY by Corollary 2.3. This leads to v<,.z. 

In fact, we have already shown that <,. is a partial order on an U-semiabundant semigroup 
S(U). It is easy to verify that <, coincides with the order w on U. 

Now the natural partial order < on an U-semiabundant semigroup S(U) is defined by 
<= — N< |. We first give an alternative description of the natural partial order < in terms of 
idempotents. 

Theorem 2.5. Let S(U) be an U-semiabundant semigroup such that U is closed under 
basic products and x,y € S(U). Then x<y if and only if there exist idempotents e and f in U 
such that « = ey = yf. 

Proof. We first prove the sufficiency part of Theorem 2.5. Suppose that « = ey = yf. 
From « = yf and Lemma 2.1 we have RY < A Choosing an idempotent xt in RU NU, we 
obtain that ¢= 2ta = (a e)y. 

Since ex = x and (a,2*) € RY, we have ext = xt. This implies (x+,e) € w” C Dg. By 
assumption, ete € U. Certainly, (x+,xt+e) € R and so (xt, rte) € RY by Corollary 2.3. Thus 
(x,xte) © RY. Hence x<,y. A similar argument shows that «<yy. 

The necessity part of Theorem 2.5 is straightforward from Theorem 2.4. 

Theorem 2.6. Let S(U) be an U-semiabundant semigroup in which U is closed under 
basic products and x,y € S(U). Then x<,y if and only if for each idempotent y+ € RO ANU 
there exists an idempotent «+ € RY QU such that xtwyt and x = x+y. The dual result holds 
for <j. 
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Proof. Suppose that TK,Y. Then RU < RY and x = ey for some idempotent e € RU NU 
by Theorem 2.4. 

Let f be an idempotent in RE NU. Then RY = RU < RU ms Re and so, by Lemma 2.2, 
Re < Ry. This leads to eS(U) C fS(U) and so e = fe giving ef € U by hypothesis. Clearly, 
(e,ef) € R which gives efRY eR” x by Corollary 2.3. Hence efwf and « = ey = (ef)y, where 
ef € RU NU. 

Conversely, suppose that for each idempotent yt € RU OU there exists an idempotent 
at € RU QU such that etwyt and « = aty. Then « = ytaty and so RY = RU ay < RY, - 
ee By Theorem 2.4, 2<,y. The proof is completed. 


§3. Locally V-semiadequate semigroups 


In this section we want to find the conditions on an U-semiabundant semigroup S(U) which 
make that the natural partial order < is compatible with multiplication of S(U). 

Recall in [2] that an U-semiabundant semigroup S(U) is called reduced if w" = w! on U. 
A reduced U-semiabundant semigroup S(U) is idempotent connected(IC) if it satisfies the two 
equations 

IC,: For any f € w(a2*) NU, «f = (af)t2; 

IC;,: For any e € w(x) NU, ex = x(ex)*. 

Lemma 3.1. Let S(U) be a reduced U-semiabundant semigroup then 


(i) If IC; holds then <, € <,; 
(ii) If IC, holds then <,.  <). 

Proof. We only need to prove (i) because the proof of (ii) is similar. If a<;y then «*wy* 
and x = yx* by the dual result of Theorem 2.6. Thus, by applying the condition [C;, we can 
obtain x = yx* = (yx*)Ty = at y. 

Certainly, y*(yx*) = yx* and so yt(ya*)* = (yx*)t. Since S(U) is a reduced U- 
semiabundant semigroup, we can easily see that 2* = (ya*)*wy*. It follows from Theorem 2.6 
that LXpY. 

Lemma 3.2. Let S(U) be an U-semiabundant semigroup in which U is closed under basic 
products and e € U. Then eS(U)e is a V-semiabundant semigroup, where V = UN E(eS(U)e). 

Proof. Let a be an clement of eS(U)e and let f be an clement of U with (f,a) € LY. 
Certainly, ae = aso that fe = f, that is, (e, f) € (w!)-1 C Dp. 

Since U is closed under basic products, the element ef € V. Clearly, (ef, f) € £ so that 
(ef, f) € LY by Corollary 2.3. It is easy to verify that (ef,a) € LY (eS(U)e). This implies that 
each element of eS(U)e is £Y-related in eS(U)e to an idempotent belonging to V. 

A similar result for RY gives us the required V-semiabundancy. 

An U-semiabundant semigroup S(U) is said to be LY-unipotent if U forms a right regular 
band. S(U) is called U-semiadequate if U forms a semilattice. 

For any e € U, we call eS(U)e a local submonoid of S(U). We shall say that S(U) 
is locally LY -unipotent(locally V-semiadequate) if every local submonoid is LY -unipotent(V- 
semiadequate). 
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A subset A of a poset (X, <) is said to be an order ideal if for each a € A and for any x € X 
with a<a then x € A(see [1]). An U-semiabundant semigroup satisfies the congruence condition 
if LY and RY are right and left congruences on an U-semiabundant semigroup, respectively 
(see [4]). 

Now we arrive at the main result of this section. 

Theorem 3.3. Let S(U) be an IC reduced U-semiabundant semigroup, in which U is 
closed under basic products, satisfying the two conditions: 


(C1) For any e € U, UN eS(U)e is an order ideal of EN eS(U)e; 
(C2) The congruence condition holds. 


Then the natural partial order < is right compatible with the multiplication if and only if S(U) 
is locally LY -unipotent, where V =UNeS(U)e. 

Proof. Suppose first that the natural partial order < is right compatible and x,y € V. 
Then x<e and so xy<ey = y. 

Thus, by Theorem 2.6, there exists f € U such that zy = yf = y(yf) = yay. It follows 
that (zy)(zy) = x(yzy) = «(xy) = xy and so that zy € ENeS(U)e. According to (Cl), 
ay € V. We have shown that V forms a right regular band. But, by Lemma 3.2, the local 
submonoid eS(U)e is V-semiabundant. Hence S(U) is locally £V -unipotent. 

Conversely, suppose that S(U) is locally £V-unipotent, that is, for any e € U, V = 
U meS(U)e forms a right regular band and a,b,c € S(U) with a<b. Then a<,b and so for 
each idempotentb+ € RY OU there exists an idempotent at € RY OU such that atwbt and 
a=artb. 

Since (be, (be)*+) € RY and b*(be) = be, we have b+(bc)+ = (bc)+. By the hypothesis that 
U is closed under basic products, (bc)tbt € UN b* S(U)bt. Certainly, ((bc) bt, (be)t) € R so 
that (bc) tb+ RY (bc)+ RY be by Corollary 2.3. 

According to (C2), (at(bc)+b+,ac) = (at(be)tbt,atbe) € RY. Since UN b+ S(U)bt 
is a right regular band and atwb*, we have at(bc)tbt € UN bt S(U)b* and at(be)tb* = 
(bc)t+btat(be)+bt. Again, ac = atbc = [at(bc)*b*]bc, where at(bc)tb+ € RY.NU. 

Thus 


cs = Riat (be)+ b+ ]bc < Rut (beyto+ = Robe)+b+at (be) t+ < Robey+ _ Bios 


It follows from Theorem 2.4 that ac<,,be. 

By Lemma 3.1, we also have <, = <;. Hence ac<be, as required. 

Combining Theorem 3.3 with its dual, we may obtain 

Corollary 3.4. Let S(U) be an IC reduced U-semiabundant semigroup in which U is 
closed under basic products. If 


(C1) For any e € U, UN eS(U)e is an order ideal of EN eS(U)e, 
(C2) S(U) satisfies the congruence condition, 


then the natural partial order < is compatible with the multiplication if and only if S(U) is 
locally V-semiadequate, where V = UN eS(U)e. 
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SDF*(n) is defined as the smallest positive integer m, such that m!! is divisible by n, where 
the double factorial m!! = 1-3-5---m, if m is odd; and m!! = 2-4-6---m, ifm is even. The 
main purpose of this paper is using the elementary and analytic methods to study the value 


distribution properties of SDF (n), and give an interesting mean value formula for it. 
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81. Introduction and results 


For any positive integer n, the famous Smarandache double-factorial function SDF (n) is 
defined as the smallest positive integer m, such that m!! is divisible by n, where the double 


factorial 


ne if m is odd ; 
mi! = 


2-4-6---m, if mis even. 


For example, the first few values of SDf(n) are: 


=1, SDF(2) =2, SDF(3) =3, SDF(4) =4, SDF(5) =5, SDF(6) =6, 
=7, SDF(8) =4,SDF(9) =9, SDF(10) =10, SDF(11) =11, SDF(12) =6, 
SDF(13) = 13, SDF(14) =14, SDF(15) =5, SDF(16) =6 ----:- , 


In reference [1] and [2], F.Smarancdache asked us to study the properties of SDF(n). 
About this problem, some authors had studied it, and obtained some interesting results, see 
reference [3]. In an unpublished paper, Zhu Minhui proved that for any real number x > 1 and 


fixed positive integer k, we have the asymptotic formula 


5a2 a? R ay x7 x? 
SDF = —_ - — O 
yy (n) ee ae a » + (<a -) ’ 


i 
nee =m ne 


where a; are computable constants. 
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The other contents related to the Smarandache double-factorial function can also be found 
in references [4], [5], [6] and [7]. For example, Dr. Xu Zhefeng [4] studied the value distribution 
problem of the F.Smarandache function $(n), and proved the following conclusion: 


Let P(n) denotes the largest prime factor of n, then for any real number x > 1, we have 
the asymptotic formula 


where ¢(s) denotes the Riemann zeta-function. 


The main purpose of this paper is using the elementary and analytic methods to study 
the value distribution problem of the double-factorial function SDF(n), and give an interesting 
asymptotic formula it. That is, we shall prove the following conclusion: 


Theorem 1. For any real number x > 1 and any fixed positive integer k, we have the 
asymptotic formula 


x és CG 0 x 
Y= (SDF(n) — P(n))? = 8) er +O (=a) ; 


a 
n<ux i=2 In’ 


where P(n) denotes the largest prime divisor of n, and all c; are computable constants. 


Now we define another function S(n) as follows: Let S(n) denotes the smallest positive 
integer m such that n | m!. That is, S(n) = min{m: n|m!}. It is called the F.Smarandache 
function. For this function, using the method of proving Theorem 1 we can also get the following: 


Theorem 2. For any real number x > | any fixed positive integer k, we have the asymptotic 


formula 


2 ¢(3) 2° Eg 23 x 
S> (SDF (n) — S(n))? = rr rae +0(—ar). 


§2. Proof of the theorems 


In this section, we shall prove our theorems directly. First we prove Theorem 1. We 
separate all integers n in the interval [1, 2] into two subsets A and B as follows: A= {n: 1< 
n<a, P(n) > Vn}; B= {n: 1<n<ax2, n ¢ A}, where P(n) denotes the largest prime 
divisor of n. If n € A, then n = m- P(n) and P(m) < P(n). So from the definition of A we 
have SDF (2) = 2. For any positive integer n > 2 and n € A, SDF(n) = P(n), if 2 7 n. 


Vol. 3 On the value distribution properties of the Smarandache double-factorial function 113 


SDF (n) = 2P(n), if 2 | n. From this properties we have 


do (SDF(n) — P(n))? 


n<u 
neA 
= © (SDF(2n)- P(2n))?+ S> (SDF(2n-1)- P(2n-1))’ 
2n<ax 2n-1l<a 
2nEeA 2n-1EA 
= = (SDF(2n)- P(2n))?= S> (2P(2n) - P(2n))’ 
nis l<n<z 
2nEA 2nEA 
= 2 6 @o= > 9 = 2). 25 P- (1) 
1<n<$ MPSS n<ZE 2N<PS Fy 
2neA p>2n ve 


By the Abel’s summation formula (See Theorem 4.2 of [8]) and the Prime Theorem (See The- 
orem 3.2 of [9]): 


k 

aj: 2 x 
= ) —+0O ‘ 
me) “¢ (<a) 


i= ln 
where a; (t= 1, 2,--- ,&) are constants and a, = 1. 
We have 
2 x ag 2 an 
SP = a-t(S)- On)? rn)-2 | "y-ray 
(2n) 2n on 
2n<pSsp 
2 (2) 

24n3 In x <4 lee n3 -In*tl gy 


where we have used the estimate 2n < \/z, and all b; are computable constants. 


=a 
Note that S- a ¢(3), from (1) and (2) we have 


n=1 


oF gaa x 
» (SDF(n) — P(n))’ ~ = Ing | Ss ne (a :) 8) 


nsx 1=2 
neA 


where all c; are computable constants. 
For any positive integer n with n € B, it is clear that SDF(n) < /n-Innand P(n) « /n. 


So we have the estimate 


s (SDF(n) — P(n))? « Son-nn< az? «In? g. (4) 
n<ax n<x 
neB 


Combining (3) and (4) we have 


S> (SDF(n) — P(n))? S> (SDF(n) — P(n))* + $> (SDF(n) — P(n))? 


n<a na n<a 
neA neB 


_ (3) 2 rae -O 8 
~ 24 Ina | : In®t1y }’ 
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where all c; are computable constants. This proves Theorem 1. 
Now we prove Theorem 2. Note that S(n) — P(n) = 0, if n € A; and |S(n) — P(n)| « Vn, 


ifne 


B. So from the result of the reference [4] and the proving method of Theorem 1 we have 


dU (SDF(n) — S(n))” = SY (SDF(n) — P(n))’ + D7 (S(n) — P(n))’ 


n<ux n<u n<ux 


—2) 7 (S(n) — P(n))  (SDF(n) — S(n)) 


n<x 


k 
¢(3) 2°? c+ x? ig 
=> > O a ao ae . 
24 mee In’ x x In*t1 


This completes the proof of Theorem 2. 
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81. Introduction and results 


As usual, n! is defined by: 


It plays a very important role in the studies of theory and application of mathematics on 
n!. Therefore the various properties of n! were investigated by many authors. For example, 


J.Stirling proved the classical Stirling formulae [1]: 


On 
n= 2rn(—)" exp(55_) (0 < On <1). 


It is well-known, any integers m(0 < m < n! — 1) can be uniquely written as[2]: 


n-1 
m= > aw, (1.1) 
i=1 


where 0 < aj <i(1 <i<n-—1). 
In this paper, we are interested in the sequence a(m) which counts the number of digit 
code in the Factorial Base. More precisely, taking m € {0,1,2,--- ,n! — 1}, we set 


n-1 
a(m) = x Qi, 
i=1 


and put 


M<x 


!This work is supported by the N.S.F. of China (Grant No. 60573040). 
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About (1.2), the authors disscused this problem in [3]. Now we shall use the different 
methods to give an exact calculating formula to (1.2) for r = 1,2. That is, we prove the 
following Theorems: 


Theorem 1. For any positive integer k, we have 


Ai(k!) = ohh ~1)-kl. 


Theorem 2. For any positive integer k, we have 


1 
Ao(k!) = take — 1)(9k? — 5k + 10) - kl. 


Theorem 3. For any positive integer N, let N = a,-ky!+aq-ko!+---+a,-ks! with 
ky < ko <+++ <ks,0 < a; < kj (1 <i < 8) under the Factorial base. Then we have 


(i) Ay(N) = a,Ai(ke!) + aly —1)-k,! + Ay(N —a,-k!) +a,(N — a, -kg!). 


(ii) Ag(N) = a, Ao(kg!) + 2(as — 1) Ar(ks!) + a(t —1)(2a, —1)-k,! + Ao(N — a, - k!) 
+2a,A1(N —a,-k!) + (N —-a,- k)a2. 


§2. Proof of the theorems 


In this section, we complete the proof of Theorems. First we shall prove Theorem 1. 
Prove of Theorem 1. For 0 <m< (k—1)!, we have 


a(m + t(k —1)!) =a(m) +t, (0<t<k-1). 


So 


m<k! 
= S- a(m) + > a(m) adil se a(m) 
m<(k—1)! (k-1)!<m<2(k—1)! (k-1)(k-1)!<m<k! 
= kA((k—1)!) hk 1)+(k=1)! 
= kAi((k-1)) + 508-1) - Bb (1) 
(2) 


Divides by k! both sides in (1), we get 


Ai(k!) — Ai((K-1)!) 1 
kis (K-11)! 2 


Ai (Rl) 


, from (2) we obtain 


Sk = Spit stk = 1), (4) 
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Note that 
Ax(1!) = S- a(m) = a(0) = 0, 


m<i! 


Ax (1! 
thus 5) = A) = 0, combining (3) we have 


1 


That is l 
A,(k!) = ake —1)-kl. 


This proves Theorem 1. 
Prove of Theorem 2. For 0 <m< (k—1)!, we have 


a?(m + t(k — 1)!) = (a(m) + #)? = a?(m) + 2ta(m) +, (0<t<k-1). 
So from Theorem 1, we have 


A2(k!) = S° a?(m) 


m<k! 
= > @m)+ S- a?(m) +e. + S- a*(m) 
m<(k—1)! (k-1)!<m<2(k—-1)! (k-1)(k-1)!<m<k! 
= kAo((k—1)!) +2A1((k-—1))(1+24+---+(k-1)) +(k-1)!- 
(17 +2?4.---+(k—1)?) 
= kAg((k —1)!) + B(&— 1) Ai((k — 1!) + E(h— 1) kh — 1) -(k-2) 
= kAg((kK-1)!)+k(k 1) Fk 1)(k — 2) -(k I+ 5 1)k(2k — 1) -(k—- 
= kAo((k I) + Fk 1)°(k 2)-R+ 2 (k 1)(2k —1)- ! 
= kAo((k—1)!) kl (hs 1)(3k? — 5k +4). 


Divides by k! both sides in (4), we get 
Aa(k!) — Ag((k-1)) 1 


2 
il (k—-i! 19k 1)(3k" — 5k +4). 
A2(k! 
Let T, = a(t ) from (5) we obtain 
1 
T= Ty + 19 (BR 8k? + 9k — 4). 
Since 
Ad(1!) = S~ a?(m) = a?(0) =0, 
m<l! 
| 
thus T, = ae = 0, combining (6) we have 


k(k — 1)(9k? — 5k + 10). 


1)! 
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No. 


That is 


Ao(k!) = Ea’ —1)(9k? — 5k + 10) - kl. 


144 


This completes the proof of Theorem 2. 


Prove of Theorem 3. Note that N = a,-k1!+ag-ko!+---+a,-k,!, and 


Ay(t- kl) = tAy(k!) + sit 1) -kI,(1<t<k). 


Then we have 


Ai(N) 


T 
M 
= 
2 


m<as:ks! 


= asAy (ks!) Ali s(t = 1) ks! + 


= asAy (ks!) ar sels = 1) ks! + 


= asAi(ks!) + site(s —1)-k,! + Ai(N —a,-k!) +a,(N — as: ks!). 


This proves (7) of Theorem 3. 


(ii) For any positive integer t(1 


Ad(tk!) = S° a?(m) 


a(m) 
as:ks!<m<N 


<t<k), we have 


L S- a*(m) +++ 


ki<m<2k! 


- S- a?(m+k!) +- 


m<k! 


m<k! 


= tAo(k!) + 2(t— 1)Ax(K!) + 


tS) (a(m) +1)? + 


m<k! 


S- a(m + as: ks!) 


0<m<N-as-ks! 


d 5 (alm) + a5) 


0<m<N-as-k! 


+ So alm) 


(t—1)k!<m<tk! 


oo S- a2(m + (t — 1)k!) 


m<k! 


+ DD (alm) + (t=)? 


m<k! 


- rt(t ~1)(2t — 1)kl. 
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So 
AN) = > a?(m) 
m<N 
= Fem Yeo 
m<as:ks! as:ks!<m<N 
= Ao(as-ks!) + S- a?(m + ag - ks!) 
m<N—as-ks! 
= Ao(as . kg!) Wa S- (a(m) + Gs)" 
m<N—as,:k! 
= Ag(as-ks!)+  S>  (a?(m) + 2a,a(m) + a2) 
m<N—as:k! 


1 
= asAo(ks!) + 2(a, — 1) Ai (ks!) 4 gules 1)(2a, — 1)k! 
+2a,A1(N —a,-k!) +(N —a,-k!)a?. 


This proves (77) of Theorem 3. 


Remark. For any positive integer r > 3, using our methods we can also give an exact 


calculating formula for A,(k!). But in these cases, the computations are more complex. 


References 


(1] I. Tomescu and S. Rudeanu, Introduction of Combinatorics, Collet’s Press, England, 


1975. 


[2] Lu Kaicheng and Lu Kaiming, Combinatorics, Beijing, Tsinghua Unversity Press, 2002. 


[3] Liang Fangchi and Jing Aiwen, Computation of k—th powers of digital sums in the 


factorial base, Joural of Air Force Engineering University(in press, in Chinese). 


ANY iy 
we) 


SCIENTIA MAGNA 


An international journal 


ISSN 1556— 6706 


